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Decays of the B meson into light and energetic particles are discussed. The calcula- 
tion of the corresponding decay amplitudes is non-trivial in many respects. Strong 
interaction effects are always present and cannot be computed reliably using an- 
alytic techniques. However, besides the intrinsic energy scale Aqcd of Quantum 
Chromo Dynamics, there also exists a much larger scale, the 6-quark mass nib, 
at which perturbation theory can be applied. QCD-factorization is the idea of 
separating the contributions that arise at these different scales and performing 
a systematic expansion in the ratio Aqcd/""^?)- Hard processes at the large scale 
can be computed using perturbation theory, while soft processes are encoded in 
non-perturbative structure functions. The observation that for many decay modes 
the same structure functions are needed make this a useful approach. However, 
factorization theorems need to be proved for every single decay, since some ampli- 
tudes do not factorize. The intent of this thesis is to study examples of factorizable 
and non-factorizable amplitudes in a systematic framework, by using effective field 
theory techniques. 

The advancing precision of experimental measurements of -B-decays make it 
necessary to improve the accuracy of theoretical predictions. To achieve this, it is 



necessary to perform a resummation of Sudakov logarithms, which enter at every 
non-trivial order in perturbation theory. 

In this thesis we present an introduction to Soft-CoUinear Effective Theory, 
which can be used to prove (or disprove) factorization theorems to all orders in 
the strong coupling constant for some B decays into light and energetic particles. 
Specifically, the factorizable amplitudes for inclusive B — > I'v and exclusive B~ 
— s> 7/"// are calculated in renormalization-group improved perturbation theory to 
first non-trivial order. Form factors encoding the exclusive decay amplitudes for 
B Pl^v and B Vl^v {P = light pseudoscalar meson, V = light vector 
meson) are studied and proved to be dominated by the non-factorizable Feynman 
mechanism. 
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CHAPTER 1 
INTRODUCTION 



1.1 Preface 

The conscious and active interplay between the observation of Nature and its math- 
ematical description is at the core of the physical sciences since the 16*^ century. In 
the quest of finding the ever-more fundamental processes, this interplay has led to 
emergence of the "Standard Model" in the second half of the 20**^ century. Our cur- 
rent theoretical understanding of elementary processes at energies testable through 
terrestrial experiments is based on the description within quantum field theories 
in combination with the gauge principle. The effects of electromagnetic, weak 
and strong nuclear forces through which the matter, quarks and leptons, interacts, 
is captured in a quantum field theory which is invariant under the gauge group 
SU (3)c X SU {2)l X U (l)y . The electromagnetic and weak forces are unified in the 
gauge group SU{2)i x U{1)y around the energy scale of roughly 100 GeV. Below 
the electroweak scale this gauge group is spontaneously broken down to f/(l)em5 
which is the group of Quantum Electro Dynamics, and the and Z bosons me- 
diating the weak interactions become massive. In the most simple mechanism that 
breaks the electroweak symmetry, the Higgs mechanism or Glashow-Weinberg- 
Salam Model, masses are also given to quarks and leptons. Fermion masses are 
generated through Yukawa interactions with the Higgs field. Using global unitary 
transformations in flavor space, the Yukawa interactions are diagonalized to ob- 
tain the physical mass eigenstates. It can be chosen such that all up-type quarks 
(up, charm, top) are not rotated in flavor space, whereas the weak eigenstates of 
the down-type quarks (down, strange, bottom) are related to the mass eigenstates 
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through the unitary matrix Vckm- Numerically, the Cabbibo-Kobayashi-Maskawa 
matrix Vckm is found to be almost diagonal, as the non-diagonal entries, that allow 
for transitions between quarks of different generations, are small. Such transitions 
are mediated by the weak currents. Since the and Z bosons have masses very 
close to the electroweak scale, it is possible to construct an effective theory for 
weak processes far below this scale. This theory is called the Fermi theory of weak 
interactions, or the weak effective Hamiltonian. 

The theory of strong interactions is Quantum Chromo Dynamics (QCD), the 
non-abelian gauge theory of SU{3)c- The condition to be renormalizable deter- 
mines the Lagrangian (almost entirely) to be 



where iD^ = id^ +gsA^ is the covariant derivative. We have suppressed any spinor 
indices and the color index a on the gluon field = A^j^T"-, which reflects the fact 
that it is a color-octet. The gluon field strength is = d^Al—d,yAl^+gsf"'''^A^^Al. 
Finally, the generators of SU {N) are denoted by with a = 1, . . . , A^^ — 1, f"-^'^ are 
the structure constants, and gs is the strong coupling constant. Every procedure to 
renormalize the theory necessitates the introduction of a dimensionful parameter 
/i. It is generally thought of as the energy scale at which the theory is defined. 
The dependence of various quantities in the theory on this scale is governed by 
renormalization-group equations (RGEs), which can be calculated in perturba- 
tively as long as the coupling is reasonably small compared to unity. Qualitatively, 
the strong coupling "constant" is large for values of fi around the intrinsic 

QCD scale Aqcd ~ 0.5 GeV, and tends asymptotically to zero for increasing fi. 
This effect is in contrast to the other sectors of the Standard Model and reflects 
the fact that quarks are confined in hadrons for low scales, but are "asymptotically 






•ip=u,d,... 
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free" for highly energetic processes such as deep inelastic scattering at scales much 
larger than Aqcd- 

All of the above is standard class-work and text-book material, which is why 
we do not feel the need to cite [1] and discuss the original works in great detail 
here. 

Corrections to the Standard Model, generally referred to as "New Physics", 
can be searched for in two complementary ways: in direct searches where new 
particles are produced and observed, or in indirect searches in which new particles 
enter as quantum corrections. The former method is obviously a very effective way 
to find New Physics, albeit a formidable financial endeavor. Indirect searches, on 
the other hand, require both precision measurements on the experimental side, as 
well as a level of theoretical Standard Model predictions that compares well with 
the experimental uncertainties. Only if both sides of the scientific medallion, the 
observation of Nature through High- Energy Experiments and its mathematical de- 
scription within the Standard Model, are met with equal precision can we conclude 
on the validity of the underlying theory. 

Needless to say, the knowledge of the parameter values entering the Standard 
Model is crucial for theoretical predictions to be reliable. The flavor sector is in this 
respect very challenging. In particular, the smallest of the CKM matrix elements, 
Vub and Vtd, are least well known, but are responsible for a wealth of phenomena, 
such as CV violation. In general, the flavor sector of the Standard Model displays 
a broad spectrum of phenomena and mechanisms that need to be understood in 
order to engage in indirect searches for New Physics. Primarily one is interested in 
the study of flavor-changing heavy quark decays, which are governed by the weak 
interactions. However, it is the strong force that is responsible for the formation of 
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the hadrons that can be observed in experiments. The B meson is a unique system 
in that it is the simplest hadron containing a heavy quark. Simple, because it is 
a bound state of only one heavy and one light (anti-) quark. A heavy quark is 
favorable for New Physics searches, because flavor-changing processes from heavy 
particles of the third family to light particles of the first two families are very 
rare. If the Standard Model predicts a quantity to be small or even vanishing, 
the chances of observing New Physics are more pronounced. Even though top 
quarks are by far the heaviest, they are less suitable because they decay before 
hadronization can occur. 

The broader aim of the "S-physics community" is the study of weak decays 
of B mesons and the interpretation of the data collected at dedicated B factories, 
such as the BaBar, Belle and CLEO experiments. The major goals are to test the 
Standard Model, extract its fundamental parameters, explore the phenomenon of 
CV violation, and to search for New Physics at higher energies. This quest will 
continue in the future with the help of the Large Hadron Collider and the LHCb 
detector in particular, as well as BTeV and possible High-Luminosity B factories. 
The theoretical challenge is the development and application of tools that enable 
us to control the effects from strong interactions. These theoretical tools include 
effective field theories, factorization theorems, symmetries, and heavy-quark ex- 
pansions. It is our belief that the present work contributes a significant step 
toward meeting this challenge. 

1.2 Methodology 

The calculation of hadronic processes is extremely difficult because of our lack 
of understanding of the connection between quark and hadron properties at low- 
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energy scales. While it is possible to apply the perturbative technique to "hard" 
processes (for example asifrib) ^ 0.22), soft processes make it difficult for analytical 
methods to make reliable predictions, since the expansion in the coupling constant 
o^sif^) = Qsif^)/^'^ longer converges. This puts some grave limitations on the 
amount of information about the quark level, that can be extracted from the 
hadronic level, and vice versa. Often times it is only possible to access such 
information if certain symmetries are present, such as the chiral symmetry for 
massless quarks, or the heavy-quark symmetry in the limit of infinite quark masses 
[2]. 

In the present work, we will concentrate on the subgroup of B decays into 
light particles only. We are thus facing the problem of having neither of the two 
symmetries mentioned above, but rather a system that involves both heavy and 
light quarks. As we shall see below, in some cases a certain "universality" emerges, 
that all soft processes are captured in a set of structure functions entering the cal- 
culation of factorizable decay amplitudes, regardless of the specific decay channel. 
To achieve this it is necessary to disentangle the short-distance effects associated 
with the large scale m^, from the long-distance effects at the low scale Aqcd- This 
is the idea of QCD-factorization, as put forward by Beneke, Buchalla, Neubert, 
and Sachrajda [3, 4, 5, 6]. It states that factorizable amplitudes can be expressed 
in the heavy-quark limit and to all orders in perturbation theory as convolution 
integrals of a perturbatively calculable "hard-scattering kernel" with some of the 
universal structure functions, which are treated as input parameters, i.e. quantities 
that need to be extracted by other non-perturbative techniques or directly from 
experiment. For practical purposes they may simply be replaced by a model. We 
shall stress, however, that this would obviously introduce a model dependence to 
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the prediction, whereas QCD-factorization itself is a model independent property 
of QCD. We were careful to add the pronoun "factorizable" to the amplitudes in 
question, because not all processes share this feature. Factorization needs to be 
proved for every individual decay channel. 

QCD-factorization combines the disentangling of physics effects from different 
energy scales with the idea of "naive factorization" of matrix elements. For a simple 
example of naive factorization consider the weak semileptonic decay amplitude of 
a neutral kaon into a charged pion and an electron-neutrino pair. The parton-level 
tree diagram mediating this decay is as follows: the strange quark couples to an up 
quark and an highly off-shell W boson, which in turn decays into a lepton-neutrino 
pair. After integrating out the W boson, the decay is mediated by GpVus times 
a four-fermion operator. The matrix element of this operator can be factorized to 
all orders in the strong coupling constant into 

(e^z>evr+| (rz/)v-A(ws)v-A IK'') = (e^z>e|(^^^)v-A|0) (7r+|(Ms)v-A|i?°) , (1.2) 

since gluons cannot couple to the lepton-neutrino pair. Note that this argument 
holds for any lepton-neutrino pair that is kinematically allowed, e.g. also for a 
muon- neutrino pair. The purely hadronic matrix element on the right-hand side 
is identified with a structure function called the transition form factor. Therefore 
equation (1.2) serves as a simple example of the "universality" mentioned above. 
Note also that an attempt to further factorize this hadronic matrix element must 
fail because of unsuppressed soft gluon exchange. The question of whether this 
is nevertheless possible for a 5 ^ tt transition form factor (exchange s b 
quarks and K ^ B states) in the limit of mi, ^ oo is a very non-trivial one and 
will be addressed in this thesis. (The answer is that the form factor does not 
factorize. However, the study reveals many surprising properties of this quantity.) 
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To see how this might possibly work consider a different example, the non-leptonic 
B Dtt decays in which the D meson picks up the spectator quark of the B 
meson, i. e. B ^ D^n' and B^ D^n^ [7]. The spectator quark and the other 
light degrees of freedom inside the B meson can easily form a D meson after the 
weak b ^ c transition. The remaining two light quarks are very energetic. To form 
a pion they must carry momenta collinear to the pion momentum and must form 
a color-singlet state. Such an energetic "color-transparent" compact object can 
leave the decay region without interfering with the formation of the D meson. The 
decoupling of soft degrees of freedom from collinear ones therefore lies at the core 
of the color-transparency argument, and likewise of soft-coUinear factorization. 

Let us now turn to the question of energy scale separation. Perturbative ef- 
fects are included in the hard-scattering kernel. In many processes one has to 
deal with multiple scales that enter the calculation. Typically, aside from a large 
scale E and the soft scale Aqcd there exists also an intermediate scale of order 
•^E'Aqcd- For the sake of the argument, let us simply state that two perturbative 
scales Ml and M2 enter the calculation of the hard-scattering kernel. Whereas one 
would naturally expand the kernel in the strong coupling constant ^^(/i) if /i can 
be chosen close to either of the two scales, it could also happen that as(yu) is mul- 
tiplied by powers of the large logarithm lnMi/M2, thus upsetting the expansion 
procedure. The effects of a single power of such a logarithm at one-loop order is 
well understood. However, double (Sudakov) logarithms appearing in the series 
Os In^ +a'^ In^ -|- . . . are more troublesome. A clean separation of scales requires the 
resummation of large logarithms to all orders in perturbation theory. 

Both issues, the soft-coUinear factorization and the separation of energy scales 
including Sudakov resummation, can be most elegantly addressed using effective 
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field theory technology. By definition the use of effective field theories achieves 
the separation of scales in that infra-red physics are reproduced in the matrix 
elements of effective operators, while ultra-violet physics effects give rise to Wilson 
coefficients. A resummation of large Sudakov logarithms is performed by solving 
the RGEs of effective operators and running Wilson coefficients from a high scale 
Hh ~ Ml down to a low scale fii ~ M2. 

The theories we have in mind are the familiar heavy-quark effective theory 
(HQET) and two versions of soft-coUinear effective theory (SCET), called SCET/ 
and SCET//. These theories have been proposed and pioneered in the year 2000 
by Bauer, Pirjol, Stewart et. al. [8, 9, 10, 11], with further development by Beneke, 
Feldman, et. al. [12, 13] and Becher, Hill, Neubert, et. al. [14, 15, 16] as of today. 
(Here we refer to the major developments only. Many more researchers contributed 
as well and were not mentioned, for which we apologize.) These theories can be 
applied to any process in which soft particles interact with light but highly energetic 
particles such as heavy-to-light decays, deep inelastic scattering and jet physics. 
In the present thesis SCET will be used to prove (or disprove) QCD-factorization 
for some exemplary inclusive and exclusive B meson decays. 

1.3 Structure of Thesis 

The thesis is structured as follows: 

We begin with a pedagogical introduction to SCET in Chapter 2. Since the 
theory is a relatively recent development within the Standard Model (its proposal 
dates only four years back), a comprehensive introduction on a basic level seems 
worthwhile. It should be understood that the presentation does not follow the 
chronology of development, but rather follows the logical steps in the construction 
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as understood in retrospect . As mentioned above, two different theories, called 
SCETj and SCET//, are needed in later Chapters. They differ in the field con- 
tent and will be discussed separately. Conceptually, SCET/ is less complicated in 
that it contains only two distinct sectors of momentum modes, whereas SCET// 
contains three and therefore displays a richer gauge structure. Before we address 
the two theories separately and in detail, a brief motivation and an outline of the 
general strategy is given. As a first application, we consider SCET currents and 
general four-quark operators, including their renormalization. These quantities are 
of special interest to any heavy-to-light decay process, in which a local operator 
product expansion (OPE) fails due to the appearance of soft and collinear singu- 
larities. We close this Chapter with a discussion on reparameterization invariance, 
which provides useful information on the structure of SCET operators and the 
large-scale dependence of their Wilson coefficients. 

Hadronic matrix elements of low-energy effective operators cannot be calcu- 
lated perturbatively. They define long-distance structure functions, whose precise 
functional form must be extracted by other means, for example Lattice QCD or di- 
rectly from experiment. However, much can be learned apart from the functional 
form using perturbative techniques. For two structure functions, the 5-meson 
light-cone distribution amplitude (LCDA) and the so-called shape function, such 
calculations are performed in Chapter 3. While LCDAs enter the factorization 
theorems of exclusive B decays, the shape function encodes the "Fermi motion" of 

the b quarks inside the B meson and is needed to factorize inclusive B decay ampli- 

^In particular we will not formulate SCET in the "label formalism" which was 
used in early papers by Bauer et. al. The formulation in position space, first used 
by Beneke et. al. and then by Neubert et. al., is equally intuitive but easier on the 
technical level. We will use the latter framework in this thesis. 
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tudes. Besides a study of renormalization properties of these structure functions, 
several novel constraints emerging from a moment analysis of the shape function 
are derived. 

As a first application we factorize the decay amplitude for inclusive B — > I'v 
in Chapter 4, using a two-step matching procedure QCD SCET/ HQET at 
next-to-leading order in perturbation theory and leading power. Such a calculation 
is appropriate in the phase-space region of large energy and moderate invariant 
mass of the Xu system. (It would be possible to avoid this two-step matching 
procedure and apply a local operator product expansion instead, only if the Xu 
invariant mass is also large. Such a kinematic setup is, however, not suitable due 
to large backgrounds, see below.) The inclusive decay is, in a sense, simpler than 
exclusive ones because the more complicated theory SCET// does not enter the 
computation. However, the matching calculation and subsequent resummation of 
Sudakov logarithms by analytically solving renormalization-group equations is far 
from trivial. To discriminate the background B X^l^v decays it is necessary 
to restrict the kinematics to a subspace in phase-space by means of certain ex- 
perimental cuts. We present results for event fractions that pass such cuts. To 
give a final numerical answer the functional form of the shape function is required. 
Here we adopt a model that is consistent with all constraints derived in the pre- 
vious Chapter 3. One of the discussed methods of kinematic cuts, the P+ cut, is 
theoretically favored and deserves a closer study. This is done at the end of this 
Section, including a thorough estimate of all theoretical uncertainties entering a 
\Vub\ determination using this method. 

In the remainder of this thesis we turn our attention to exclusive decays. In 
Chapter 5, the outhned methodology is applied to B ^ 7^"^, which is one of 
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the simplest exclusive decay modes in that there are no hadrons in the final state. 
Still there is sensitivity to the light-cone structure of the B meson, because of the 
coupling of the high-energy photon to the soft spectator quark inside the heavy 
meson. Many techniques can be learned in this environment. With minor modi- 
fications they can then be applied to the hard-scattering term of such important 
decay amplitudes as the ones for B K*'~f oi B ^ vrvr, etc. One of the important 
remaining contributions to the corresponding amplitudes involves B ^ M form 
factors, where M is some light meson. While it is not the goal of this thesis to 
study these difficult decay modes in all completeness, we wish to contribute by 
analyzing some of their ingredients in a clean environment. 

For that reason we devote the Chapter 6 to analyse B ^ M form factors in 
the high recoil region, meaning that the light meson M is highly energetic. The 
underlying processes are B nl^u and B pl^v. The calculation of these 
form factors using a two-step matching procedure QCD — SCET/ SCET// 
reveals some almost expected outcomes (such as the fact that form factors do not 
factorize, and that form factors at large recoil are power suppressed quantities) 
and also some surprises. The most important one might be that, although going 
through SCET/ enables us to resum large perturbative logarithms, there exists 
a purely long-distance contribution for which the intermediate theory is without 
any physical significance. In other words, the non-perturbative "soft overlap" or 
Feynman mechanism (of lots of small momentum kicks to bring the spectator quark 
up to speed) exists and dominates over the hard-scattering picture. Furthermore 
we find that there are unsuppressed contributions from higher Fock states of both 
the B and M meson. More precisely, three-particle configurations (two quarks plus 
one gluon) contribute at the same power as two-particle ones. This is in contrast to 
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the previous understanding as found by e. g. light-cone sum rules. It is tempting to 
try to estimate the contribution of this effect; however, this would be accompanied 
by large uncertainties. Note also that we only perform the matching at tree-level, 
due to the complexity of the problem. A numerical analysis is therefore left for 
future work and will not be given here. Finally we conclude in Chapter 7. 



CHAPTER 2 
SOFT-COLLINEAR EFFECTIVE THEORY 



2.1 General Considerations 

Many processes in 5-decay physics involve hadrons or jets of hadrons with energies 
much larger than their masses. The assumption, that underlies the construction 
of an effective field theory describing the dynamics of the process, is that the 
constituents of these hadrons or jets carry momenta coUinear to the hadronic mo- 
menta. Consider a light but energetic hadron with momentum in the ^-direction 
and invariant mass = m'jj much smaller than its energy E. We may decom- 
pose P^ = {E — mjj/AE + . . + {m'jj/AE + . . .)n'^ with the light-like vectors 
= (1, 0, 0, 1) and = (1, 0, 0, —1) satisfying = = and n - n = 2. Up to 
small corrections of order / E^ the coefficient of is given by the large energy 
E of the hadron. A natural description of these kind of momenta is given in terms 
of the light-cone decomposition 

P^ = in-P)- + {n.P)'^ + P^. (2.1) 

Four-vectors are characterized by the scaling of the individual light-cone compo- 
nents. The three terms in (2.1) are referred to as P+, P_ and Pi. In other 
words, P+ is a four-vector P'^ = {n ■ P) n''/2, etc. In the example above, we have 
(P^, P^, P'l) = (^ f^, 2E ^, 0) + . . . as the leading contributions. We assume 
that a parton inside this hadron carries momentum p^, which also points essentially 
in the n'^ direction, but it may also have non-vanishing momentum in the n and 
perpendicular direction. Those are, however, dynamically generated and typically 
much smaller than the large minus component. By introducing a dimensionless 
scaling parameter A we may characterize the order of the light-cone components 

13 
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~ (A^,1,A). The scaling of the "plus component" p+ ~ follows from the 
assumptions that the invariant mass of the constituent p"^ = 2p^ ■ p^ + p\ is not 
larger than O(A^). All the above assumptions are well justified for the processes 
described below. For example, consider a decay B tt in the kinematic situ- 
ation where the B meson is at rest and the pion takes away almost half of the 
5-meson mass as kinetic energy E ^ Mb/2. The quarks and gluons inside the 
pion are massless, but off-shell by parametrically ^qcd- Throughout this work, 
Aqcd denotes the soft QCD scale of non-perturbative physics. We would define 
A = Aqcd / E and assume n ■ p = 2xE where x is a positive dimensionless variable 
less than (but of order) 1. Generally, the precise definition of the scaling parame- 
ter A depends on the particular process and should be thought of as the ratio of 
two scales whenever we are facing a multi-scale problem. The effective theory is 
constructed around the systematic expansion in inverse powers of the larger of the 
two scales. We shall distinguish between labels and dynamical components. While 
large momentum components appear as labels and are not altered by soft inter- 
actions (only by collinear interactions), the remaining ones are dynamical. This 
is similar to the notion of the velocity label v and "residual" momentum k in the 
HQET expansion pb = rribV + k [2]. For the sake of simplicity, we will not write 
out explicitly the label v on heavy quark fields or the label En on collinear fields. 
In fact, we will not use the label formalism of earlier papers, and keep the large 
momentum dependence explicit. 

The relevant degrees of freedom of SCET are such that for any given perturba- 
tive process in QCD involving the kind of momenta described above, each infra-red 
contribution is reproduced by the effective fields. In this spirit, SCET is founded 
on the method of regions [17]. However, the approach of going through the effec- 
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tive theory allows for a systematic power counting already at the beginning of a 
calculation. 

There is a certain spin-symmetry realized in the large energy limit. This can 
be understood when decomposing the collinear QCD spinor ip^ into two separate 
fields ^ and rj, which are both two-component spinors and subject to the constraints 

= and ^rj = 0. One defines 

e = |^,. „ = (2.2) 

The power counting associated with these two-component spinor fields is chosen 
such that propagation of the fields do not enter the power counting of any given 
diagram, i.e. count as unity. Denoting the time-ordered product by T, the propa- 
gator reads 

{0\T{M.UM} 10) = / • (2-3) 

When assuming collinear momentum scaling ~ (A^, 1, A), the integration measure 
d'^p = ^dp^dp-(Pp± scales like 0{X^), and ~ C){X^). It follows that the indi- 
vidual fields count as ^ ~ 0{X) and rj ~ O(A^). In the limit of large energy, the 
4-component QCD spinors ipc reduce to a two-component spinor ^, thus giving rise 
to a reduced Dirac basis. This is the manifestation of large-energy spin-symmetry 
[18]. Similarly, analyzing the gluon propagator in covariant gauge 



{0\T{A^,{x),A:{y)}\0) 



dp 



^-ip-ix-y) ^2.4) 



(27r)4 p2 + 

reveals that A'^ counts in its light-cone components like the momentum itself. 

After the identification of the small component ?7, the SCET-Lagrangian, which 
contains only the leading collinear field is constructed by integrating out rj. Let 
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'i/'c be the full QCD quark field which carries coUinear momentum. 

= 4^^^, = {^ + r])(^zn-D^ + zn-D^ + ip^ + r,) (2.5) 

The small component rj is integrated out by solving its equation of motion 
6Cf^^/6f] = 0. The solution is 

V = -- ^ ip±^- (2.6) 

' 2 in - D + ie ^ ' 

The SCET Lagrangian is therefore identical to the QCD Lagrangian, but boosted 

to a reference frame in which the fields carry collinear momenta. As a result, the 

coUinear Lagrangian is exact to all orders in A and not renormalized [12]. At this 

stage we have 

C^ = ^l(tn-D + ip^—^—-tPj) e • (2.7) 
2 \ in ■ U + le J 

Note that the Lagrangian itself scales like A^, which is consistent with the action 
S = J d*x Cc ~ since the integration measure df^x ~ 0(A~^) for collinear 

momenta. The factor n - D m the denominator is symbolic. After the introduction 
of Wilson lines, it leads to an integration of the fields along the n direction. This 
point will be explained later on in more detail. We have also inserted a +ie pre- 
scription, which is an arbitrary choice and not dictated by the QCD Lagrangian 
[12]. This regulator is necessary for practical purposes but bears no physical im- 
plications. From now on we will omit this prescription in the notation and only 
refer to it explicitly when needed. 

The next step toward a consistent power counting is to achieve a "homoge- 
neous" scaling in A, i.e. an expansion. In the above formulation, D denotes the 
full covariant derivative, which means that the gluon fields can carry any momen- 
tum. As mentioned before, the strategy for the construction of SCET is to assign 
an independent effective field for each momentum configuration that contributes 
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to a given process at leading power. Depending on the process, there can be many. 
Typically, we work in a reference frame in which there are "soft" momenta present, 
e.g. fields that carry momenta with small light-cone components. The largest soft 
mode that is kinematically allowed to couple to a collinear field without changing 
its momentum scaling is an "ultrasoft" field with momentum ~ (A^,A^,A^). Let 
us assume for the moment, and for the sake of simplicity, that only these two 
momentum configurations contribute. (Such a theory is called SCET/.) It is then 
necessary to split up iD^ = id^ + + M^^. The subscripts c and us label the 
collinear and ultrasoft fields, respectively. Ultrasoft gluons scale like ultrasoft mo- 
menta, and it is therefore necessary to expand the Lagrangian in powers of A for 
a homogeneous scaling. This, and also the so-called "multipole expansion" below, 
provide a rigorous definition of leading and subleading Lagrangians in terms of 
inverse powers of the large scale E. 

Many ingredients of SCET, including the definition of the expansion parameter 
A in terms of physical quantities, and the number and nature of degrees of freedom 
in the effective theory, depend on the particular process that is analyzed. In the 
next sections we will construct two different theories called SCET/ and SCET//, 
suitable for inclusive and exclusive heavy-to-light decays. The heavy system pro- 
vides for a natural reference frame in which to carry out the calculations, namely 
its rest frame. The constituents of the heavy system contain both heavy and light 
soft fields. In the effective theory, the mass of the heavy system is a large scale and 
should not enter the low-energy description. This is provided by the heavy-quark 
effective theory, which will be part of SCET. In the decay applications considered 
here, the fe-quark mass mf, and the large energy E of the collinear fields will be of 
the same order. 
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Table 2.1: Nomenclature for various momentum modes. Quark fields car- 
rying these momenta scale according to the last column, whereas 
gluon fields scale in their light-cone components like the corre- 
sponding momenta. The effective field theories will contain a sub- 
set of these fields, depending on the application. 
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Before going into detail, let us set the nomenclature for the various fields and 
their scaling. The names are tied to the scaling in A, regardless of the physical 
definition of A. They are listed in Table 2.1. Not all fields listed here will be 
included in the final low-energy theory. Nevertheless, it is convenient to define a 
set of momentum modes which allow to address contributions to processes as they 
are found by the method of regions. This can be demonstrated in the following 
example [15]: Consider a scalar triangle graph in the kinematic setup where the 
external momenta are /'^ ~ (A, A, A) soft, ~ (A^,1,A) coUinear, and = (/ — 
p)^ ~ (A, 1, A) hard-collinear, as shown in Fig. 2.1. [q^ is called a hard-collinear 
momentum even though the perpendicular component is smaller than \/~X. Only 
the first two components, "plus" and "minus" are of importance.) One defines the 
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Figure 2.1: An Example: Scalar triangle graph with external momenta I 
(soft) and p (collinear). The loop momentum is denoted k. 



loop integral 

^ J {k^ + zo)[{k + iy + io][ik+py + zo] 

in d = 4 — 2e space-time dimensions and analyzes it for arbitrary external momenta 
obeying the above scaling relations. It will be convenient to define the invariants 

L'^ = -P-iO, P^ = -p'^-iO, Q"^ = -{l-py-i0 = 2l+-p_-i0 + ... , 

(2.9) 

which scale like ~ A^, ~ A^, and Q'^ ~ A. (This P has nothing to do 
the momentum defined in (2.1), where P served as an example for the light-cone 
decomposition of a 4- vector.) The exact result is 



1 



1 Q\ ^V^Vna^ 



(2.10) 



The next step is to reproduce the result using the method of regions [17], which 
means that one assumes a certain scaling for the loop momentum k and expands 
the integrand in powers of A before performing the integration. Going through 
the list of modes in Table 2.1 we find leading power contributions for k being 
hard-coUinear, collinear, soft, and soft-coUinear. The individual regions give the 
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following contributions [15]: 



'hc = zvr ' /i 



2 1 2 9 

In — + - In^ 

1 



'^'^^^^^ + -ln4 + iln^4-;!V0(e), (2.11) 



r(l + e) / 1 1 , /i^ 1 , 9 /i^ 7r2\ ^ ^ ^ 

^ In^-oln'^ + ir +0(6). (2.12) 



(P + zO) (2A;_ ■ /+ + zO) [(A; + p)2 + zQ] 
ilni^-iln^i^ 

g2 \^ e2 g p2 2 P2 

'^/^u^-'^ / d^k ^ 



(P + iO) [(A; + ly + iO] (2A;+ ■ p. + iO) 



r(l + e) / 1 1 , /i^ 1 , 9 /i^ ^ ^ ^ 

^ ^' ln^--ln2^ + - +0(e). (2.13) 



g2 y £2 g ^2 2 L2 ' 6 



(P + iO) (2P ■ /+ + /2 + iO) (2fc+ • p_ + p2 + io) 



r(l + e) / 1 1 /i2Q2 1 ^ ^2g2 



TT 



A couple of comments are in order. 

(i) It is a simple check to add all the above contributions and find that they 
indeed reproduce the full result / = /hc + -^c + -^s + -^sc- 

(ii) Each individual integral is, of course, Lorentz invariant. In this sense 
Lorentz invariance was never broken by picking a particular frame. Had we picked 
a different reference frame, for example the Breit frame, the results would not 
change. However, we would have found different momentum modes that con- 
tribute. Clearly, we could perform a boost to this frame before evaluating the inte- 
grals, thus finding a different set of momentum modes. In the example of the Breit 
frame, the formerly soft leg / becomes coUinear in the n direction (anti-collinear), 
and the coUinear momentum p stays coUinear in the n direction. Leading contribu- 
tions would come from hard loop momenta (formerly hard-collinear), anticollinear 
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Figure 2.2: Corresponding contributions in the effective low-energy theory. 

5C denotes a Wilson coefficient, the solid lines are soft fields, the 
dashed lines coUinear, and the dotted line soft-coUinear. 

(formerly soft), coUinear (formerly collinear) and ultrasoft (formerly soft-coUinear). 
It is therefore not only the physical process that defines which modes must be in- 
cluded in the effective theory, but also the frame in which we choose to work. 
One further point can be made by boosting the frame even further, such that the 
formerly collinear momenta become soft. Clearly, the formerly soft momenta then 
become collinear in the n direction and we find the following "symmetry" (this is, 
of course, nothing else but Lorentz invariance): The contribution from soft and 
collinear loop momenta are related to each other under the simultaneous exchange 
of soft and collinear external momenta and the interchange of n and n. However, if 
a heavy quark is present there is a preferred Lorentz frame, namely the rest frame 
of the heavy quark, and the above symmetry is broken. 

(iii) We went through a lengthy discussion of the various contributions to the 
full theory diagram. The benefit is the separation of short- and long-distance 
physics. The strategy for constructing a low-energy effective theory is to let 
the collinear, soft, and soft-collinear contributions be reproduced by exchange of 
collinear, soft, and soft-collinear fields, while the hard-coUinear contribution in this 
case will be absorbed into a Wilson coefficient. In fact, the sum of the collinear, 
soft, and soft-collinear integrals is free of infra-red singularities, which are regular- 
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ized by keeping the external legs off-shell (e.g. the sum is free of ln(/i^/L^) and 
ln(/x^/P^) ). The view of associating the infra-red contributions with diagrams in 
the effective theory is depicted in Fig. 2.2. Note that when a soft gluon connects 
to a coUinear line it will lead to a (hard-collinear) off-shell mode, which will be 
integrated out. Therefore the soft gluon emerges out of the vertex in the figure. 
A similar argument applies to the collinear gluon. A deeper investigation of this 
point will follow in later sections and lead to Wilson lines. Here we wish to stress 
the close connection between the method of regions and the idea of constructing 
SCET clS db low-energy effective theory. 

Before beginning the discussion of SCET/ in greater detail, note that the 
same (physical) momentum modes may have different names, depending on the 
(physical) definition of the expansion parameter A. Consider, for example, a 
momentum p ~ (Aqcd, E, a/ -EAqcd) with some scales E ^ Aqcd- This mode 
can be called collinear when defining A = a/AqcdZ-E. It can, however, also 
be called hard-collinear when defining A = Aqcd/E. Similarly, a momentum 
/ ~ (Aqcd, Aqcd, Aqcd) might be called ultrasoft or soft for A = a/Aqcd/-E or 
A = Aqcb/E, respectively. 

2.2 Renormalization 

For consistency we review the main strategies of renormalization using dimensional 
regularization [19] in 4 — 2e space-time dimensions. This scheme has already been 
adopted in the previous section and will be used throughout this work. As apparent 
in (2.11), for example, loop integrals which are logarithmically divergent in four 
dimensions become finite for e 7^ 0. The general strategy is to absorb terms that 
diverge in the limit e — > into factors that relate bare quantities to renormalized 
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ones. In the end, renormalized quantities are finite for e ^ and liave meaningful 
relations to physical observables. 

Since the action J d'^~'^^x C{x) is dimensionless, it follows that the bare QCD 
coupling constant 5'^^''^ has a mass-dimension of e. Defining {g^^'^'^Y = Ana^^^^, the 
relation to the dimensionless renormalized coupling ^^(/i) can be written as [20] 

«^'-^ = /i2^Z,(«,(/i))«,(/i) , (2.15) 

where /i is an arbitrary mass parameter, called the renormalization scale. In a 
minimal subtraction scheme, the renormalization factor Za = 1 + ^~ 

Za is 

k=l 

expanded in inverse powers of e, and the evolution of as{f^) is encoded in the "beta 
function" 

■a,(yu) = /5(a,(/i),e) . (2.16) 



dln/i 

We will consistently work in the modified minimal subtraction scheme MS, in 
which all terms multiplying powers of 1/e — 7^; + ln47r are absorbed into the 
renormalization factors. Since the beta function (3{as, e) is finite in the limit e — 
(in a renormalizable theory), it can be Taylor expanded around e = 0. It follows 
[20] that only the first term in this expansion is non-zero and the exact relation 
reads 

l3{as,e) = P{as) -2eas , 

d ^2.17) 
/3(a.)=2a^ — Z«(a.). 

Throughout this work we will use the following expansion of the beta function in 
terms of a^: 



with [21] 



/?(°) = ll-^n;, (2.19) 
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/5(^^ = 102-f n,, 

(2.20) 

2857 5033 325 o 
^ 2 18 ^ 54 ' 

where Uf is the number of flavors that contribute in the loop calculation. 

We will also need the wave-function renormalization of the quark fields. In 
QCD, the (massless or massive) bare and renormalized quark fields are related 
by = with Z^, = 1 - a,/(37re). The only quark field in SCET 

that renormalizes differently is the heavy quark field. Its wave-function renor- 
malization is obviously identical to the heavy quark field in HQET and reads [2] 
Zh = l + 2aJ{37Te). 

Composite operators Qi require renormalization beyond that of its field and 
coupling components, and are renormalized in the same manner. Under renormal- 
ization they can mix among each other if they share the same quantum numbers 
and dimensions. We will only consider operators that are gauge invariant and do 
not vanish under the application of equation of motions. Therefore one relates 

g^-<= = z,,(/i)g,(/i) (2.21) 



with Zij = 6ij +^ e ^Z^'^'^ in the MS scheme. The yU. dependence of an operator Qi 

k=l 

is encoded in the renormalization group equation (RGE). Since the set of all allowed 
operators {Qi, Q21 • • •} is a complete basis, one can express the operator dQi/d\\i fi 
as a linear combination of Qi, Q2, ■ ■ ■■ Defining the anomalous dimension matrix 
7 through the RGE 

^ Qj = -ljkQk, (2.22) 



dln/i 

it follows from the fact that the bare operator Q^^^° is independent of /i, that 

d „ ^ d ^(1) 



25 

Here we have used (2.16) and (2.17). In words, the anomalous dimension is, 
apart from a minus sign, twice the coefficient of the 1/e term, when all ultra-violet 
divergences are regularized dimensionally and all infra-red effects are regularized 

by other means. For any anomalous dimension 7(0^) we will adopt the expansion 

00 , 

k=0 

unless otherwise stated. 

Physical amplitudes are given as products of operator matrix elements and 
Wilson coefficients Ci{fi) {Qi{fi)) and are independent of the renormalization scale 
H. This implies that the Wilson coefficients must obey the RGE 

(2.25) 

where T denotes the transpose of the matrix 7. The above equation will be used 
below to resum large logarithms to all orders. 

2.3 SCETj 

This theory is applicable to processes in which there are only two different mo- 
mentum modes present, namely hard-collinear and soft ones. An example for such 
a process is the inclusive decay of a heavy meson such as the B meson. In the rest 
frame of the B meson the partons carry, with the exception of the massive b quark, 
momenta that scale like (Aqcd, Aqcd, Aqcd), with Aqcd being the soft QCD scale. 
The b quark momentum can be split up into the static component mi,v and the dy- 
namic residual momentum k, which also scales like (Aqcd, Aqcd, Aqcd)- A flavor- 
changing current turns the b quark into some light but energetic quark. This quark 
could be a parton of a particular hadron or, more generally, of a hadronic jet of 
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invariant mass much larger than Aqcd- We are then faced with the kinematic 
situation in which a hierarchy of scales emerged: Aqcd <^ M <^ E', where E is the 
energy of the jet (close to half of the 6-quark mass) and M is the jet invariant mass. 
Specifically, we assume that M is of the "intermediate"^ order M ~ -^/^Aqcd- In 
this kinematic setup the two momentum modes are therefore (Aqcd, Aqcd, Aqcd) 
and (Aqcd,-^', a/^Aqcd)- When choosing the expansion parameter A = Aqcd/-^ 
we can refer to them as soft and hard-coUinear, respectively. (Alternatively, they 
are also referred to as ultrasoft and coUinear when choosing A = a/AqcdZ-E.) 

2.3.1 Power Counting 

The full massless QCD fields are separated into hard-coUinear and soft momentum 
modes, e.g. the quark field splits up into the soft quark field q and the hard- 
coUinear quark field ip^c = ^hc + Vhc- Analogously, the QCD gluon field A is 
decomposed, yielding 

4^ = ^hc + r]hc + q, ^D^' = + ^7^^, + gA^^ . (2.26) 

Let us study the effect on the kinetic term in the Lagrangian C = ip ip ip + . . .. 

Since the spinors ^hc and rjhc are constrained by i/iC^hc = "^Vhc = 0, we decompose 

the Dirac matrix in light-cone coordinates to n'^yi /2 + n^i^ /2 + 7^ and expand 

the Lagrangian. 

— — ^ 7^ _ 

ipipip = ^h^-in- D ^hc + f]hc-in- D rjhc + Vhc %pi.ihc + ihc ipiJ]hc 

^ ^ (2.27) 

^Ihc g4^hc q + v gAhc q + qg4^hc ihc + qg^^hc Vhc + q{i^ + g4s) q 

^For instance, in inclusive B —>■ decays we need to restrict the jet invariant 
mass to be less than the charm mass for discrimination purposes. Numerically, rric 
compares very well with y^m^^AQco- 
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Note that we do not allow for interaction terms that are forbidden by kinematics, 
for example two soft fields cannot couple to a single hard-coUinear field. The small 
component rjhc is then integrated out by solving its equation of motion. This yields 

Vhc = ~ {^P± ihc + 94^hc q) , (2.28) 

which must be plugged back into (2.27). So far we have merely rewritten the 
Lagrangian into a more complicated but exact form. In a next step one needs 
to perform a systematic expansion in A of this Lagrangian. To do this, it is 
important to take the scaling of the integration measure in S* = J d^x L into 
account. (So far, we omitted the explicit dependence of the various fields on 
the position x. We will continue to do so, unless otherwise stated.) For terms 
involving only soft fields, d^x scales like A~^. To see this consider a soft field 
at position y. Translational invariance relates this field to one at position x + y 
through an exponential exp(?/ ■ x), where / is the soft momentum, and I ■ x = 
Z+ ■ a:_ + /_ ■ x+ + /_L ■ a;_L. Because of the scaling properties of the soft momentum, 
leading contributions to the action are obtained when x scales like (A~^, A^^, A^^). 
Similarly, for terms involving hard-coUinear fields, x scales like (1, A^^, A~^/^) and 
the integration measure d^x ~ A~^. 

Inverse covariant derivative operators are also expanded. For brevity, let us 
introduce the intuitive notation Dhc = d — igAhc (and Dg = d ~ igAg will also be 
used below). When {in ■ D)^^ acts on hard-coUinear fields, we must expand 

= - gn-Ag .-i— + O(A^) . (2.29) 

m- D m- Dhc ■ Dhc in ■ Dhc 

Every term in the effective theory should have a single and homogeneous scaling 

behaviour. This requirement leads to a subtlety in interaction terms involving soft 

and hard-coUinear fields. In the above notation we have omitted the dependence 
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on the position and treated all terms as local in that all fields are evaluated at the 
same position x. According to the above argument, the light-cone components of 
X in interactions of hard-coUinear and soft fields scale like (1, A~^, A~^/^), because 
the vertex carries hard-collinear momentum (since two hard-collinear momenta 
can not add up to a soft momentum). One must therefore perform a "light-front 
multipole expansion" [13] of the translation operator with —id ~ (A, A, A) 

(j)s{x) = exp(x ■ (9) 0s(O) 

= exp(x+ ■ + x± ■ d±) exp(x_ ■ d^) 0^(0) 
= 4>s{x-) + [x± ■ d±(f)s] (x-) 

+ ^x+ ■ [d_<Ps] {x-) + ^ [x'ix'id.dM {x-) + ... (2.30) 

for any soft field 0^ in such interaction terms. This means that all soft fields must 
be evaluated at a position x^ = {n ■ x) n^/2 on the light cone, whereas the hard- 
collinear fields remain at position x. This non-locality is a refiection of the fact that 
at leading power in A the total momentum in such interactions is not conserved. 
However, to put this shocking observation into the right context, we emphasize 
that translational invariance will be restored order by order in the expansion in 
A. At leading power only the "plus component" of the soft momentum is kept 
when adding it to a hard-collinear momentum: (A, 1, A^^^) + (A, A, A) = (A, 1, A^^^). 
The expanded Lagrangian is then homogeneous in the power counting term by 
term. This concludes the systematic expansion in powers of A. We will address 
the important question of gauge invariance in a separate section below. 

Let us turn our attention to the treatment of the h quark in SCET. Heavy 
quarks are described by the familiar HQET Lagrangian [2]. In the low energy 
effective theory, the heavy quarks do not couple to hard-collinear fields because 
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such an interaction would put the quark off-shelL The HQET Lagrangian reads 



Ch = hiv ■ Dg h + 



h{iD,fh + C^ag{^i)h^-a^,G^-h +0{l/ml). (2.31) 



2m;, 

Finally the pure glue (Yang-Mills) Lagrangian £giuc can be split into soft plus 
hard-coUinear plus interaction terms. The soft and hard-collinear sectors take the 
same form as in full QCD, including gauge-fixing and ghost terms. The interac- 
tion terms are restricted to those that are kinematically allowed, using the same 
arguments as above. 

2.3.2 Wilson lines 

The introduction of Wilson lines is a well known concept in HQET. Consider the 
effect of coupling a gluon with momentum of order O(Aqcd) (soft) to a heavy 
quark field with momentum mi,v. This will lead to a propagator that is off-shell 
by an amount {mi,v + — rnf ^ 2mbV ■ k ~ O(mfeAQCD)- If the effective theory 
describes physics of energies below this scale, such propagators should be integrated 
out. We can repeat this argument consecutively and find that attaching n gluons 
with momenta ki, i = 1 .. .n, to a heavy quark line leads to 



V ■ {ki + . . . + kn) " ' V ■ {ki + k2) V ■ ki 
These are precisely the Feynman rules for the object 



(2.32) 



Sy{x)=Pexp^g J ds V ■ As{x + sv)^ , (2.33) 

which is called a Wilson (or Schwinger) line extending from minus infinity to the 
interaction point x along the v direction. The path-ordering symbol P denotes 
that the gluon fields are ordered from left to right in order of decreasing s. In the 
spirit of effective theories, we integrate the off-shell propagators out by performing 
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the field redefinition h{x) Sv{x) h^^\x). The Wilson line Sy has some very 



useful properties, SlS^ = S^Sl = 1 and, most importantly, Sliv ■ Dg = iv ■ d. 
Performing the field redefinition leads to a free-particle Lagrangian at leading 
power, Ch = h^^^ iv ■ dh^'^^ + . . ., so that the "sterile" field h^^^ no longer couples to 
soft gluons. In practice we will use the Wilson lines 5"^ only when discussing the 
renormalization properties of operators that include heavy fields. For explicit loop 
calculations one can choose to work with either h or Svh^^\ since they lead to the 
same Feynman rules. 

The purpose of this recapitulation is to draw a close analogy to the treatment 
of hard-coUinear gluon attachments, which can be dealt with in the same manner. 
In full QCD the attachment of a gluon with hard-coUinear momentum En + ... to 
a massive quark line carrying momentum rribV will provide for a Dirac structure 



where the dots represent power corrections. An infinite succession of such attach- 
ments results in the Wilson line Whc [10], which can be written in position space 
as 



We should point out that the above (somewhat historic) discussion does not sug- 
gest that the massive b field matches onto the product (Whch). In fact, it is the 
combination with a hard-coUinear field, {^hc^hc), that appears in explicit match- 
ing calculations, e. g. (ipTb) — > {^^c Whc T h). As discussed below, this combination 
is also dictated by hard-coUinear gauge invariance. 

We also need the conjugate operator Wl^, which has the opposite path ordering 
and the factor ig replaced by —ig. The most important operator properties are 




(2.34) 




(2.35) 
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that WljVhc = WhcWl^ = 1 and in ■ Dhc Whc = in - d, following from the fact 
that {in- Dhc Whc) = by definition. As a consequence, one obtains the very useful 
identity [11, 12, 14] 

1 1 



(x) = Whcix) Wlix) 0(x) = -t I ds Whcix 



in - Dhr in - d 





'x + sn) 



(2.36) 

for any field (or product of fields) 0(x). This enables us to write the purely 
hard-collinear quark Lagrangian Che (the interaction terms with soft fields will be 
denoted by £mi) in the following form: 

Chc{,x) = ^hc{x) - in - Dhc ^hc{x) (2.37) 





i I ds 



ihclphclWhc (X)| WllphcUhc 



[x + sn) . 



This Lagrangian sums up an infinite number of leading-order couplings between 
hard-collinear quarks and hard-collinear gluons. For later purposes, note that 
Whc = 1 in "coUinear light-cone gauge" n - Ahc = 0. 

One can also introduce the soft Wilson line 5*^, analogous to the hard-collinear 
Wilson line, with n ^ n interchanged [11, 12]. Again, = 1 if we choose the 
n - As = Q gauge. Although Ss is not necessary to construct the Lagrangian 
for soft fields, it turns out to be very useful in understanding how soft-collinear 
factorization arises as a property of the effective theory at lowest order. This 
is similar to the decoupling of HQET fields h from soft fields, which leaves the 
heavy field /i*-"^ sterile to gluon interactions. It has been shown [11] that the field 
redefinitions 

ihc = Ss C,hc ; ^hc = Ss ^ii^ Sl , c= Ss c^^^ Sl , (2.38) 
remove all hard-collinear - soft interactions from the leading-order Lagrangian. 
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(Above, c denotes the ghost fields.) We will find (and explain in some more detail) 
an analogous mechanism for the decoupling transformations in SCET// below. 



2.3.3 Gauge Transformations 



Gluon fields have been split up into two components, hard-coUinear and soft. The 
gauge symmetry of QCD also decomposes into hard-collinear and soft gauge sym- 
metries. It is clear that soft fields do not transform under hard-collinear gauge 
transformations Uhc, since this would change their momentum scaling. Only hard- 
collinear fields can transform. On the other hand, soft gauge transformations Ug 
do not alter the kinematic scalings, and therefore all fields transform, including 
hard-collinear ones. For that matter, soft gluons can be viewed as slowly vary- 
ing background fields. Hard-collinear gluons transform then covariantly under soft 
gauge transformations, while they transform inhomogeneously in the hard-collinear 
sector. 



hard-coll.: 



soft: 



Ahc 
As 

Ahc 
As 



Uhc Ahc Ul + -Uhc 
9 

As , 

UsAhcUl , 



UsAsUl + -Us[d,Ul] 
9 



Ds , Ul 



he 



ihc — ^ Uhci. 

q^q , 

ihc^ Us ihc 

q-^Usq . 



(2.39) 



Note that the sum Ahc + Ag transforms in the usual way under both hard- 
collinear and soft transformations. We have also introduced the hard-collinear 
Wilson line Whc{x)i which extends from {—oo) to x, and therefore transforms like 
Whc{x) Uhc{x) Whc{x) Ul^{—oo). Let us agree that the gauge transformations 
Uhc{x) become irrelevant for x — — oo and set Uhc{—oo) = 1. This will obvi- 
ously not alter the transformation properties of finite Wilson lines Whc{0) Wl^{tn) 
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along the light cone n, because the transformations Uhc{—oo) cancel. It will, how- 
ever, simplify the discussion since now the Wilson lines transform like Whc{x) 
Uhc{x) Whc{x) under hard-coUinear gauge transformations. Note that, as a simple 
but important example, the combinations 



Xhc = Wl a. , Al = [Wl iDl Wnc\ , (2.40) 

are left invariant. In light-cone gauge n ■ A^c = 0, the calligraphic gluon field 
Ahc is identical to the QCD gluon field, which will make matching calculations 
particularly simple. We will sometimes refer to the calligraphic fields as "gauge- 
invariant building blocks" , because it is easy to construct gauge-invariant operators 
using these fields [14]. 

The effective Lagrangian before multipole expanding is invariant under the 
above gauge transformations. However, such a transformation reintroduces terms 
of different power counting in the Lagrangian. In the end, we wish to construct 
the theory in such a way that every term in the Lagrangian has a simple (ho- 
mogeneous) scaling behaviour. This leads to the notion of "homogeneous gauge 
transformations", and is discussed in detail in [13]. We have also explored that in 
interactions with hard-collinear fields, soft fields need to be positioned at x_ rather 
than X. A new Wilson line extending from x to a;_ is then necessary to restore 
gauge invariance. These Wilson lines need to be expanded in A, too, which leads 
to new interaction terms in the subleading Lagrangians. The above arguments will 
be expanded and discussed in more detail in the case of SCET//. 



2.3.4 Heavy-to-Collinear Currents 

An important application for SCET is the decay of the heavy b quark to a highly 
energetic light quark q. Such flavor- changing currents are the result of the weak in- 
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teractions and appear in the weak effective Hamiltonian description of the Standard 
Model at energy scales well below the weak scale ~ 100 GeV. The QCD current 
is qTb, where we treat F as an arbitrary Dirac structure. At leading power in A, 
this current matches onto the SCET current Xhdx) T h{x^). As noted before, the 
dynamical heavy field is treated as a soft particle and is therefore positioned at 
x_. The Wilson line in the definition of the gauge-invariant field X^c = ^Ic^hc 
correctly accounts for an arbitrary power of n ■ Ahc gluon insertions. Only the large 
components n ■ Ahc need to be exponentiated since all other components are para- 
metrically smaller and can be expanded. Soft gluon insertions enter at subleading 
power. Although these subleading currents are known in the literature [12, 22, 23], 
they are not needed for this work at present (although very important for future 
continuations of it). There are two categories (labeled as "A" and "B") which are 
separately invariant under reparameterization transformations (see section 2.5). 
The leading current Xhc{x) T h{x-) belongs to the "A" category and at order A^^^ 
one finds, for example, 

J(^i) = XhcTx^-D,h- Xhc i^hci. 

m - d ^ (2.41) 

Irrii) 

The calligraphic covariant derivative is defined as iD^c = id + Ahc- AH hard- 
coUinear fields are located at position x, and all soft fields at position x_. "A"- 
and "B"-type currents are also known at order A. For the remainder of this section 
we will consider the most simple case and discuss the matching at leading power 
only. 

The effective fields Xhc and h are two-component spinors, in contrast to the 
four-component spinors q and b. Thus, the Dirac basis {1, 75, 7^^, 7^75, ia'^'^, ia'^'^'-f^} 
will be matched onto a smaller basis, which is commonly chosen as {1, 75, 7^}. Let 
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us define the currents 

Ji = 'Xhch , J5 = Xhc 75 ^ , ■J± = y^hc 1± h , (2.42) 

and the following symmetric and anti-symmetric tensors in the transverse plane 

9T = 9'" - ""'""^ \ , 6^ = e^-'v^n, , (2.43) 

for convenience. We use eoi23 = 1 and 75 = i'^'^^^^'^x' ■ 

Omitting the dependence on the renormalization scale and denoting the Wilson 
coefficients by Cj, the matching of heavy-to-light currents onto the above SCET 
currents reads to all orders in the strong coupling constant and at leading power 

qh ^ CiJi, (2.44) 
^756 C2J5, 

q^^h ^ C3J^ + (C4n'^ + C5t;'^)Ji, 

q,a>^^b ^ Cg{n'^gl'-n^gf)Ji-C^o^eTJ5 

+Cn {W - v-'nn Ji + C12 iv'^g"/ - v'gf) Ji , 

+C10 {v^'rf - v^n^") J5 - C12 {v^'ief - v^^ief) Ji . 

This follows from the reduction of the Dirac basis. Furthermore, it has been 
shown [9] that in the naive dimensional regularization scheme of anticommuting 
75 in 4 — 2e dimensions, there are relations Ci = C2, C3 = Cq, = C7, C5 = Cs, 
Cio = Cii, C12 = 0, among the coefficient functions that hold true to all orders for 
massless light quarks. At tree-level, the coefficient functions are 

Ci, 2,3,4,6,7,9,10,11 = 1 7 ^5,8,12 = . (2-45) 
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One-loop corrections to these Wilson coefficients have been computed and can be 
found in [9]. Apart from the renormalization scale /i they depend on the large scales 
n ■ p and mb. This matching calculation achieves the ffist step in scale separation, 
and we will typically need some combination of these Wilson coefficients in later 
applications (e.g. factorization of the hadronic tensor in inclusive B Xul^v 
decays) . 

In this short section we do not present the radiative corrections, and therefore 
we will also delay the discussion of renormalizat ion-group flow of the coefficient 
functions to later chapters, i.e. when we need it. Let us just mention without 
proof that all Wilson coefficient functions Cj share the same anomalous dimension 
[9], which will be the same (up to a minus sign) as the anomalous dimension of 
the SCET^ current "Xhc T h. 

2.4 SCETjj 

The goal is to construct an effective theory in which momentum fluctuations are 
at most of order the soft QCD scale Aqcd- The main application for such a theory 
is the calculation of exclusive Heavy-to-Light decay amplitudes. Introducing the 
expansion parameter as A = AqcdZ-E'? the constituents of the B meson have a 
soft momentum scaling (A, A, A). If the decay product is a light energetic hadron, 
we assign coUinear scaling (A^, 1, A) to the partons inside this hadron. These are 

the only momenta that appear on external legs. However, as we have seen in the 

^We will see later that the anomalous dimension is independent of the hard- 
coUinear momentum square p^. This must be true, because serves as an infra-red 
regulator, which does not enter the ultra-violet behaviour encoded in the anoma- 
lous dimension. Therefore the anomalous dimension will be the same in both 
theories SCET/ ijp' ~ .EAqcd) and SCET// {p^ ~ K^qy))^ and we will discuss it in 
the latter theory. 
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introductory example in Fig. 2.1, there is an additional long-distance mode (soft- 
collinear) which scales like (A^, A, \^^^). Such a field communicates with either the 
soft or the coUinear sector but does not appear on external legs, and is therefore 
referred to as a "messenger" . Whether or not this mode will contribute in the end 
depends on choices such as infra-red regulators and quark masses. We shall work 
consistently in the limit of vanishing quark masses for the first family, and therefore 
include this mode in the effective theory. We will see later that the soft-coUinear 
sector of SCET// provides us with very elegant arguments to discuss important 
physics such as factorizability and endpoint singularities of decay amplitudes. 

Let us deal with the various sectors of SCET// one by one. The effective 
Lagrangian contains interactions between soft and collinear fields by either the 
exchange of messenger modes or induced interactions resulting from the exchange 
of hard-coUinear off-shell propagators. The latter are, however, kinematically for- 
bidden in the applications considered in this work [15], because they require ex- 
ceptional momentum configurations which are typically not generated at leading 
power in loop graphs. We may therefore neglect such terms. The Lagrangian can 
be split up as 

C = Cc + Cs + Ch + Csc + Ant . (2.46) 

The construction of the collinear, soft, and heavy Lagrangians proceed analogously 
to our previous discussions [12, 9]: 

A = dzn-D,^-^zp,^ I e ^ ^ q^p^q+Ct^ , (2.47) 

2 2 tn ■ Dc 

and the HQET Lagrangian as given in (2.31). Note that, in contrast to the con- 
struction of SCET/, the above Lagrangians are exact to all orders in A. The gluon 
Lagrangians in the three sectors retain the same form as in full QCD, but with 
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the gluon fields restricted to the corresponding subspaces of their soft, collinear, or 
soft-coUinear Fourier modes. The soft-colhnear Lagrangian resembles in its form 
the collinear Lagrangian. This is not surprising, since the light-cone components of 
a soft-coUinear momentum displays the same hierarchy of scales. For completeness 
(and at the same time reviewing the strategy for £c); let us repeat the derivation 
briefly. Using the projectors and the quark field Qsc = d + <^ splits up 

into small and large components. Analyzing the propagator reveals that ^ ~ is 
the large component, and a ~ A^/^ can be eliminated using the equation of motion 

(J = -t:-^^^Psc±e. (2.48) 
/ m ■ Use 

Inserting this result back into the Lagrangian yields the exact result 

Csc = e^in-Dsce-e ipsa. % ipsc^ • (2.49) 

/ 2 tn ■ Use 

The most important new ingredient of SCET// is the interaction sector between 
soft-colhnear fields with soft or collinear degrees of freedom. The formalism follows 
along the same lines as the discussion of hard-collinear - soft interactions in SCET/, 
but is richer (and more interesting). It has been developed and discussed in detail 
in [15] and will be outlined in the next section. 

In interactions with other fields, the soft-collinear fields (but not the soft and 
collinear fields) are multipole expanded as 

0sc(3^) = (t^sc^xJ) + x_L ■ 9_L 0^,c(x„) + . . . in collluear interactions , 

(2.50) 

0sc(a;) = 0sc(a^+) + x± ■ d±_ (psci^Xj^) + . . . in soft interactions . 

As before, the four- vectors x_ = [fi ■ x)nj1 and x+ = {n- x)n/2 are light-cone 
positions. The first correction terms in (2.50) are of 0(A^/^), and the omitted 
terms are of 0(A) and higher. 
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2.4.1 Gauge Transformations and Interaction Terms 

Soft-coUinear fields can couple to soft or collinear fields without altering their 
scaling properties. This motivates the treatment of the soft-coUinear gluon field 
as a background field. However, in order to preserve the scaling properties of the 
fields under gauge transformations one must expand the transformation laws in A. 
In close analogy to the discussion in SCET/, this leads to the following set of 
"homogeneous" gauge transformations for the quark fields: 

soft: Qsix) — > Us{x) qs{x) , coll. and soft-coll, fields invariant 
collinear: ^(x) — > Uc{x) C,{x) , soft and soft-coll, fields invariant 

soft-coUinear: g^(x) Uscix+) qs{x) , ^(x) ^ Usc{x-) ^{x) , 

Qscix) Usc{x) Qscix) . 

(2.51) 

The transformation laws for gluons are more complicated, but straightforward [15]. 
Treating the slowly varying soft-collinear gluons as background fields, we have the 
following transformation laws under soft and collinear gauge transformations before 
multipole expansion: 

Us Uj + Us [zD^,, Uj] , A^.^U, A^f f/] + U, [zD^^, t/t] , (2.52) 

while the collinear and soft-collinear fields remain invariant under Us, and soft and 
soft-collinear fields remain invariant under Uc- Finally, under soft-collinear gauge 
transformations Use the fields transform as 

A^^^UscA^^Ul, ^^Usc^, 

A^: ^ Use A>: Ul , qs ^ Use Qs , (2-53) 

^ Use Ule + Use [^^^ f/jj , g,, ^ Use qsc ■ 
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It can be seen from these relations that the combination (A^ + A^^) transforms 
in the usual way under both soft and soft-collinear gauge transformations, while 
(A^+A^^) transforms in the usual way under both collinear and soft-collinear gauge 
transformations. In order to preserve a homogeneous power counting, we expand 
the above rules and keep only the leading-order terms for consistency. This means 
that the transformation laws are replaced by the homogeneous transformations 

soft: n ■ As ^ UsU ■ AsUl + Us[in ■ d,Ul] , qs-^UsQs, 

n-As^Usn-AsUl + [m ■ Dsc{x+), Uj] , 

soft-collinear: Uscix+) A'^ Ul^{x+) , f/.c(a;+) Qs , 

(2.54) 

for soft gluon fields and 

collinear: n ■ A^ UcU ■ A^Ul + Uc[in ■ d,Ul] , ^^U^^, 

n-Ac^Ucn-AcUl + Uc [in ■ D^dx^), Ul] , 

soft-collinear: A^f ^ Usc{x^) Ul{x.) , ^^Usc{x-)^, 

(2.55) 

for collinear gluon fields. The soft-collinear gauge sector transforms in the usual 
form. Let us now return to the interaction Lagrangian. There are no interactions 
between soft and collinear fields, since such interactions would lead to off-shell 
momenta of 0{E Aqqj)). Soft and collinear fields can, however, couple separately 
to soft-collinear fields. At leading power we have 

^hli^) = Qsix) I gn ■ Ascix+) qs{x) + h{x) gn ■ A,c(x+) h{x) 

+ ^{x) — gn ■ Asc{x^) ^{x) + pure glue terms . (2.56) 
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Subleading terms have been calculated and can be found in [15]. 

Momentum conservation implies that soft-collinear fields can only couple to 
either soft or coUinear modes, but not both. More than one soft or collinear 
particle must be involved in such interactions. The gluon self-couplings can be 
derived by substituting ^ + |n'^n ■ Asc{x+) for the gluon field in the 
soft Yang-Mills Lagrangian and A'^ —>■ A'^ + ^n'^n ■ Asdx^) for the gluon field 
in the collinear Yang-Mills Lagrangian, and isolating terms containing the soft- 
collinear field. The precise form of these interactions will not be relevant to our 
discussion. Finally, let us note that none of the terms in the SCET Lagrangian 
(2.46) is renormalized beyond the usual renormalization of the strong coupling and 
the fields [12, 15]. 

The multipole expansion of the soft-collinear fields implies that momentum is 
not conserved at these vertices. When a soft (light or heavy) quark with momen- 
tum ps absorbs a soft-collinear gluon with momentum k, the outgoing soft quark 
carries momentum Ps + k_. Likewise, when a collinear quark with momentum 
Pc absorbs a soft-collinear gluon with momentum k, the outgoing collinear quark 
carries momentum pc + 

In order to match the quark and gluon fields of the full theory onto SCET fields 
obeying the homogeneous gauge transformations one first adopts specific gauges 
in the soft and collinear sectors, namely soft light-cone gauge n ■ As = (SLCG) 
and collinear light-cone gauge n ■ Ac = (CLCG). At leading order in A, one then 
introduces the corresponding SCET fields via the substitutions [13, 15] 

il^sls^^.^ ^ Rs Sl qs , b\^^^^ ^ Rs Sl h , ■0c|clcg ^ ^ • {'^■^7) 

The corresponding replacements for gluon fields can be found in [15]. They are a 
little more complicated, because the homogeneous gauge transformation mixes the 
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soft and the coUinear gluon fields with soft-collinear ones. For example, the two 
non-zero components of the soft gluon in SLCG are given by 



n ■ As\^^^^ ^ Rs [Slitn ■ D'^^ 5,) + Sl[n ■ ^c(a;+), Ss]]RI . 



(2.58) 



Similarly for coUinear gluons. The quantities 



Ss{x) = Pexp^g J dt n ■ As{x + tn)^ 
Wc{x) = P exp J dtn ■ Ac{x + tn)^ 



(2.59) 



are SCET Wilson lines in the soft and collinear sectors [12, 11], which effectively 
put the SCET fields into light-cone gauge. Small subscripts "s" and "c" are pro- 
vided to distinguish these Wilson lines from corresponding ones with the soft or 
collinear gluon fields replaced by soft-collinear ones. The objects Rg and Rc are 
short gauge strings of soft-collinear fields from x+ to x (for Rg) and a;_ to x (for R^). 
They differ from 1 by terms of order A^^^ and so must be Taylor expanded. Note 
that Ss transforms as Ss{x) Us{x) Ss{x) and Ss{x) — > Usc{x+) Ss{x) f/]^(x_|_) un- 
der soft and soft-collinear gauge transformations and is invariant under collinear 
gauge transformations. Likewise, Wc transforms as Wc{x) Uc{x) Wc{x) and 
Wc{x) Usc{x_)Wc{x)Ul^{xJ) under collinear and soft-collinear gauge trans- 
formations and is invariant under soft gauge transformations. The short strings 
only transform under soft-collinear gauge transformations, in such a way that 
Rs{x) Usc{,x)Rs{x)Ul^{x+) and R^^x) f/,c(x) i?c(x) t/],(x_). It follows that 
the expressions on the right-hand side of (2.57) are invariant under soft and 
collinear gauge transformations and transform as ordinary QCD quark fields under 
soft-collinear gauge transformations. 

Again it is sometimes useful to introduce the calligraphic fields that are gauge- 
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invariant building blocks 

Wl i = Ssc{x. ) X , SUs = WscM Qs , Slh = Wscix+) J{ , 

(2.60) 

where fields without argument live at position x. The quantities Wsc and Ssc are 
a new set of Wilson lines defined in analogy with Wc and Ss in (2.59), but, as 
mentioned above, with the gluon fields replaced by soft-collinear gluon fields in 
both cases. The calligraphic fields in (2.60) are invariant under soft, coUinear, 
and soft-collinear gauge transformations. This follows from the fact that the new 
Wilson lines transform as 

Wsc{x+) ^Usc{x+)Wsc{x+), S^x-) ^ Uscix.) Ssc{X-) (2.61) 
under soft-collinear gauge transformations. 

2.4.2 Currents 
Heavy-to-Collinear Currents 

Let us recall from the discussion of Wilson lines that Whc emerged in heavy-to- 
collinear currents in SCET/ by keeping leading-order attachments of collinear glu- 
ons to the heavy quark line. In SCET// we denote this Wilson line by Wc- Note 
that attachments of soft gluons to the collinear quark line will also lead to off-shell 
modes, and thus to the appearance of Sg. Naively, one might therefore expect 
to build up the current "S"! E Wc h, and in a non-abelian theory Sl and Wc do 
not commute. However, adding diagrams involving non-abelian gluon couplings 
reverses the order of the two Wilson lines [11], as dictated by soft and collinear 
gauge invariance. 
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The matching procedure (2.57) reproduces this exphcit finding in the following 
simple way. At tree-level and leading power the QCD current is matched onto the 
gauge-invariant object (omitting the large HQET phase e~^"^'''"'^) 

Mx)rbix) ^ [^WcRl]{x)r[RsSlh]{x) 

= [^Wc]{x+ + x^)T[Slh]{x.+x^) + 0{\). (2.62) 

Note that the expression in the first line is not homogeneous in A. In interactions of 
soft and coUinear fields, the soft fields must be multipole expanded about a;+ = 0, 
while the collinear fields must be multipole expanded about x_ = 0. Also, as 
mentioned above, the quantities Rs and Rc must be expanded and equal 1 to first 
order. This leads to the result shown in the second line. The terms of 0(A^/^) in the 
expansions of Rs and Rc cancel each other [15]. The leading-order SCET current 
in the final expression is gauge invariant even without the Rg and Rc factors, since 
soft fields at x+ = and collinear fields at x_ = both transform with Usc{0) 
under soft-coUinear gauge transformations. 

When radiative corrections are taken into account, the current mixes with 
analogous operators at different positions on the light cone, and for the case of 
the heavy-to-coUinear currents different Dirac structures can be induced by hard 
gluon exchange. The correct matching relation reads (setting x = for simplicity) 

V'c(0)r6(0)^^ j dsC,{s,fi)[^Wc]{sn)r, [Slh]{0) + OiX) 

' (2.63) 

= c^{n ■ P^ ^^) [iWc r. si h] (0) + 0(A) , 

i 

where translation invariance is used to rewrite the expression in a local form. As in 
the case of SCET/, there are only three independent Dirac structures Fj allowed. 
The coefficient functions Ci{n ■ P'^) = J dse^'^^'^''Ci{s) are the Fourier transforms 
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of the position-space Wilson coefficients Ci and depend on the total coUinear mo- 
mentum operator P'^ = P^^^. — P;^. In fact, by reparameterization invariance (see 
Section 2.5 below) the dependence of the Wilson coefficient functions on P'^ must 
be through the combination 2v ■ Pf = {v ■ n){n ■ P'^). Unless otherwise stated, we 
set t; ■ n = 1 for convenience, and work with the above expressions. 

Soft-to-Collinear Currents 

Arguments along the same line as in the previous section would suggest that the 
QCD current ijjc{^) T il)s{Q) matches simply onto terms of the form T \Sl q\ 

with the fields evaluated at different positions. This is, however, not quite correct. 
Note that the soft massless quark field is a full four-component QCD spinor, so that 
we cannot use the reduced Dirac basis Fj. A careful perturbative matching analysis 
shows that there is a second structure that appears at leading power, involving a 
perpendicularly polarized coUinear gluon [14]. At tree-level, the expression 

[eiye](o)r[5tg](o) - f dtm^ 

sums up an infinite set of leading-order interactions. Despite of the extra perpen- 
dicular derivative (which scales like A), this operator is not subleading because 
of the large non-localities associated with soft fields (integration in t over a large 
domain of order A~^). Radiative corrections mix these currents with similar oper- 
ators at different positions on the light-cone. Note that because of the factor of yi, 
the operator vanishes if F commutes with ^. The existence of this second operator 
therefore depends on the particular Dirac structure. For the present work we will 
only need Dirac structures that are of the form F^, so that the second term in 
(2.64) vanishes since = 0. The matching relation of interest therefore reads 



(0) [si q] (tn) 



(2.64) 
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^,(0)r^^,(0) ^ [ dsdtD{s,t,fi) [^W,]{sn)T^ [5] g,] (tn) + 0(A) 

^ J ^ (2.65) 

= D{P^ .Pl,f,) [^W^T^ Si qs] (0) + 0(A) . 
Again, reparameterization invariance dictates that the Wilson coefficient can only 
depend on the scalar product 2P^ ■ = {n ■ P^)(n ■ P^). 

Radiative Corrections and Anomalous Dimensions 

The momentum-space coefficient functions in (2.63) and (2.65) are renormalized 
multiplicatively and obey renormalization-group (RG) equations of the Sudakov 
type [9, 24], i.e. contain an explicit logarithmic renormalization scale dependence, 
d 



Ci{v-pc-,n) = -f^h{v ■ Pc-, IJ') Ci{v ■ Pc-, 1^) , 

^ D {Ps+ ■Pc-,^)= l^q {Ps+ ■Pc-,fJ')D{ps+- Pc- , /i) 



(2.66) 



dln/i 

where the anomalous dimensions take the form 

1 /i^ 
lih{v ■pc-,fi) = -- rcusp[a.(/i)] In y + T^h[as{fi)] 

1iq{Ps+-Pc-,f^) = -rcusp[«s(/i)] In- + r^g[asifi)] . 

Pc- 



(2.67) 



The coefficients of the logarithmic terms are determined in terms of the universal 
cusp anomalous dimension Fcusp = CFOis/rc + 0{a1), which plays a central role 
in the renormalization of Wilson lines with light-like segments [25]. The one-loop 
expressions for the non-logarithmic terms in the anomalous dimensions can be 
deduced from the explicit results for the Wilson coefficients derived in [9, 14]. 
They are 

%(a.) = -\^ + 0{al) , r,,(a.) = J^ + o{al) . (2.68) 

4 TT 2 TT 

It may seem surprising that after hard and hard-coUinear scales have been 
integrated out the operators of the low-energy theory still know about the large 
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Figure 2.3: SCET graphs contributing to the anomalous dimension of a 
soft-coUinear current. Full lines denote soft fields, dashed lines 
coUinear fields, and dotted lines soft-coUinear fields. 

scales V ■ pc- ~ E and Ps+ ■ Pc- ~ EAqcb, as is evident from the appearance 
of the logarithms in (2.67). The reason is that in interactions involving both 
soft and coUinear particles there is a large Lorentz boost 7 ~ Ps ■ Pc/ ^/vl vl ~ 
E'/Aqcd connecting the rest frames of soft and coUinear hadrons, which is fixed by 
external kinematics and enters the effective theory as a parameter. This is similar 
to applications of HQET to 6 — ^ c transitions, where the fields depend on the 
external velocities of the hadrons contaiping the heavy quarks, and 7 = w;, ■ = 
0(1) is an external parameter that appears in matrix elements and anomalous 
dimensions of velocity-changing current operators [2, 26]. In the case of SCET, 
the boost parameter 7 can, in some cases, reintroduce a "long-distance type" 
dependence on the large energy through low-energy operator matrix elements (see 
for example Section 6). 

We will now explain how the results for the anomalous dimensions can be 
obtained from a calculation of UV poles of SCET loop diagrams [16]. The rel- 
evant diagrams needed at one-loop order are shown in Fig. 2.3. They must be 
supplemented by wave-function renormalization of the quark fields. The gluons 
connected to the current are part of the Wilson lines Wc and Ss- We regularize 
IR singularities by keeping the external lines off-shell. The results for the sum of 
all UV poles must be independent of the IR regulators. For the heavy-to-coUinear 
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current the pole terms obtained from the three diagrams are (here and below we 
omit the —iO in the arguments of logarithms) 



(2.69) 



\ e iv ■ pc- /i / e ■ pc^ 2e 

while for the soft-to-collinear current we obtain 

2 2 3 \ / 2 2 -pI 3 \ 

e 2ps+-pc-li^J e 2ps+-Pc- e 

We quote the contributions to the operator renormalization constants in units 
of Cpo.s/''^'^ (in the MS subtraction scheme in 4 — 2e dimensions). The three paren- 
theses in the first line of the above equations correspond to the soft, coUinear, and 
soft-collinear contributions, where the first two terms include the corresponding 
contributions from wave-function renormalization. The 1/e poles of the soft and 
coUinear graphs depend on the IR regulators, but this dependence is precisely 
canceled by the soft-collinear contribution. By construction, the sum of the soft, 
coUinear, and soft-collinear contributions is IR finite and only contains UV poles, 
whose coefficients depend on the ratios v ■ Pc-/ ^J'' and ps+ - Pc~l This follows 
since IR divergences in both the full and the effective theory (which are equiva- 
lent at low energy) are regularized by the off-shellness of the external quark lines. 
The one-loop contributions to the anomalous dimensions 7^/1 and 7^5 are given by 
—CfolsI2'j\: times the coefficients of the 1/e poles in the above expressions. They are 
in agreement with the results (2.67) and (2.68) obtained from the scale dependence 
of Wilson coefficients. 

The calculations presented above make it evident that there is an intricate in- 
terplay between the soft, coUinear, and soft-collinear diagrams. In dimensional 
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regularization the dependence of the anomalous dimensions on the hard or hard- 
coUinear scale enters through the loop integral involving the soft-collinear exchange 
and thus seems to be related to very small momentum scales. However, care must 
be taken when assigning physical significance to the scales associated with individ- 
ual diagrams in SCET, because in the soft and collinear diagrams a cancellation of 
IR and UV divergences takes place. The logarithms appearing in their divergent 
parts should be interpreted as [cf. (2.70)] 

lnZ^ + ln^ = ln^+ln^, with = ^^^^^^ , (2.71) 

and thus arise from a cancellation of physics at the hard scale = 2ps+-Pc- and at 
the soft-collinear scale m^^. In the sum of all graphs, the soft-collinear contribution 
precisely cancels the IR piece of the soft and collinear parts, see (2.70). This 
interpretation is consistent with the fact that the anomalous dimensions measure 
the change of operator matrix elements under infinitesimal variations of the UV 
cutoff /i. They are therefore insensitive to the physics at low scales by construction, 
and the large logarithms in (2.67) are really of short-distance nature. 

We would like to close this Section with a (somewhat historic) remark that 
dates back to the time before soft-collinear fields were introduced to SCET/j. 
After all, for exclusive B decays one identifies the expansion parameter A with the 
ratio Aqcb/E, and therefore the soft-collinear scale m^c is parametrically smaller 
than Aqcd- This raises some concerns as to whether the soft-collinear fields have 
a right to exist. So consider for a few seconds (and not much longer) the soft- 
collinear sector to be absent. Then the following puzzle arises: We know from the 
explicit expressions for the Wilson coefficients of the currents that their anomalous 
dimensions must depend on a scalar product of the collinear momentum with a 
momentum characterizing the soft quark (i.e., v ■ or Pc~)- However, the 
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SCET Feynman rules imply that the first graph in Fig. 2.3 can only be a function 
of the soft momentum (i.e., v ■ Ps or p^), while the second one can only depend 
on the coUinear momentum (i. e. , p^), as is in fact confirmed by the explicit 
calculation. The apparent "factorization" of soft and collinear degrees of freedom 
in SCET// (in the hypothetical absence of soft-coUinear fields) would thus lead to 
the conclusion that the anomalous dimensions of the currents are independent of 
the products v ■ pc^ or ps+ ■ Pc-, in contradiction with the results for the Wilson 
coefficients. In fact, factorization would be a generic property of SCET//, since the 
interaction Lagrangian in (2.56) would no longer exist. This is obviously incorrect; 
for example it is well known that 5 — > tt transition form factors do not factorize, 
although they can be described using SCET//. 

Although the above argument strongly motivates the use of soft-collinear fields 
in SCET//, it does not necessitate it. Alternative prescriptions include the intro- 
duction of "analytic regulators" [27] where propagators are raised to a non-integer 
power, and explicit infra-red regulators in the Lagrangian [28]. However, these 
schemes are not favorable in our opinion because analytic regulators are very com- 
plicated and cannot be introduced on the Lagrangian level in a gauge-invariant 
fashion, and the proposed infra-red regulators in [28] break many important sym- 
metries such as gauge invariance and Lorentz invariance. The soft-collinear sector 
solves the problem in a much more simple and elegant way. 

2.4.3 Four-quark Operators 

In many cases in i?-decay physics, amplitudes calculated in SCET receive con- 
tributions from hadronic matrix elements of four-quark operators, which can be 
expressed in terms of the leading-order light-cone distribution amplitudes of a light 
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final-state meson and of tlie initial state B meson. Tlie relevant SCET operators 
can be taken as [14, 16] 

Qic){s,t) = [^W,]{sn)^T^T,[W^^^]{0) [qsSs]{tn)^T,T,[Slh]{0) 

POO rn-P 

= / due-''^' dae''"Q^c){^,a), (2.72) 
Jo Jo 

where the color label C = S* or O refers to the color singlet-singlet and color octet- 
octet structures Ti ® T2 = 1 1 or ® T4, respectively. Again, n ■ P is the total 
momentum carried by all coUinear particles, which is fixed by kinematics. (Strictly 
speaking, this is a momentum operator.) In light-cone gauge, u = n-ps corresponds 
to the plus component of the momentum of the spectator anti-quark in the B 
meson, while a = n-p^ denotes the minus component of the momentum of the quark 
inside a light final-state meson. It is conventional to introduce a dimensionless 
variable u = a/n ■ P & [0, 1] corresponding to the longitudinal momentum fraction 
carried by the quark. 

The momentum-space operators Q(c)(cu,cr) obey the rather complicated RG 
integro-differential equation 

d 



Q{C){^,un- P) = - duj' rfM'7(c)(cJ,u;', n-P, /i) (5(c) (u;',m' n-P). 
am/x Jq Jq 

(2.73) 

To obtain the anomalous dimensions at leading order we compute the 1/e poles of 
the diagrams shown in Fig. 2.4 in dimensional regularization and add the contribu- 
tions from wave-function renormalization. Note that only soft-collinear gluons can 
be exchanged between the soft and coUinear currents. For the color-singlet case 
Ti ® T2 = 1 ® 1 we find that the sum of the four diagrams with soft-collinear ex- 
changes (but not each diagram separately) is UV finite. The anomalous dimension 
is then a combination of the anomalous dimensions for the two non-local currents 
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Figure 2.4: SCET graphs contributing to the anomalous dimension of the 
four-quark operators Q(c)(a;,(T). Full lines denote soft fields, 
dashed lines coUinear fields, and dotted lines soft-coUinear fields. 



in (2.72). At one-loop order we obtain 



6{u-u')V{u,u)+6{u-u)H{u,u',fi) , (2.74) 



where (with u = 1 — u 
V{u,u) = 



u 



^ c{T{)\ e{u'^- u) + ^ f + c(ri) ) e{u - u') 



u \u — u 



1 - ciVi) 
+ ^^(^-^) 



(2.75) 



with c(l) = 0(75) = 1, c(7^) = is the Brodsky-Lepage kernel [29] for the evolu- 
tion of the leading-twist light-cone distribution amplitudes of a light meson, which 
we have reproduced here using the Feynman rules of SCET. The plus distribution 
is defined as 

[fiu,u')]+ = f{u,u) - 5{u-u') j dwf{w,u), (2.76) 

^0 
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which coincides with the conventional definition if the distribution acts on functions 
g{u) but not if it acts on functions g{u'). The function 



H{u, u', /i) = (^In - ^ 5{lu -u')-u 



(2.77) 



uj{uj — uj') uj'{uj' — uj) 

is the analogous kernel governing the evolution of the leading-order 5-meson light- 
cone distribution amplitude [30]. Here the plus distribution is symmetric in the 
two arguments and defined as 

du' [/(cu, uj')]+ g{uj') = / du' f{u, u') [g{u') - g{u)] . (2.78) 
Jo 

For the color-octet case Ti (g) T2 = T4 ® T4 things are more complicated. In 

that case the sum of the soft-coUinear exchange graphs shown in the last line in 

Fig. 2.4 does not vanish. However, it provides contributions such that, in the end, 

the dependence on the IR regulators drops out. Our final result for the anomalous 

dimension in the octet-octet case is 

1 «s r 

7(0)(c^,c^',u,u',ri ■ P,/i) = -^-^ 6{uJ-uj')V{u,u')+6{u-u')H{uj,uj',fi) 

2 71 \ n-vuj{n-Fy 4 / 

2.4.4 Decoupling transformation 

The leading-order interactions between soft-collinear fields and soft or collinear 
fields in the SCET Lagrangian can be removed by a redefinition of the soft and 
collinear fields [15]. In analogy with the decoupling of soft gluons in SCETj [11], 
we define new fields 

qs{x) = W,,{x+) qf\x) , h{x) = iy,,(x+) h^^\x) , = ^,,(x_) i^^\x) , 

A^(x) = WscM '^(x) WUx+) , A(f(x) = Ssc{x.) 4°)'^(x) Sl{x.) . 

(2.80) 
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The quantities Wsc and Ssc have already been introduced in (2.60), where we 
studied their gauge transformation properties. The new fields with "(0)" super- 
scripts are invariant under soft-coUinear gauge transformations, but there is no 
longer a soft-collinear background field in their transformation law. Essentially, 
the leading power effect of the above field redefinition is that the substitutions 
A'^ A'^ + ^n^^ n ■ Asc{x+) and Ai^ A^ + ^n'^n ■ Asc{x-) from the multipole 
expansion are reversed. As a result, when the new fields are introduced in the 
SCET Lagrangian the terms Cg, Cc, and Cgc retain their original form, while the 
leading-order interaction Lagrangian vanishes! Residual interactions between 
soft-collinear and soft or collinear fields start at 0(A^/^) [15]. After the field redef- 
inition it is convenient to introduce the gauge-invariant building blocks [14, 15] 

5ft (x) (x) = WUx+) Slix) q,{x) = Q,(x) , 
Si^^\x) = Wl{x+) Slix) h{x) = :K{x) , (2.81) 

lyf t(^) ^(o)(^) ^ slix^) W^x) e(x) = X(x) , 

which are obviously identical to (2.60) and invariant under all three types of gauge 
transformations. 

The fact that interactions of soft-collinear fields with other fields can be de- 
coupled from the strong-interaction Lagrangian does not necessarily imply that 
these fields can be ignored at leading order in power counting. The question 
is whether the decoupling transformation (2.80) leaves external operators such as 
weak-interaction currents invariant. The analysis of the previous sections indicates 
that in some cases the soft-collinear exchange graphs contribute to the calculation 
of the anomalous dimensions. Let us then study what happens when the decou- 
pling transformation is applied to the various types of operators. 

Under the transformation (2.80), the heavy-to-collinear and soft-to-collinear 
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currents in (2.63) and (2.65) transform into (setting x_l = for simplicity) 

w-c] (x+) r [si h] (x_) ^ x(x+) ^J,(o) r w^o) n^-) , 

(2.82) 

[e lyj (x+) r [si q,] (x_) ^ x(x+) sJ,(o) r w^M Q.(x„) . 

We observe that the soft-colhnear fields do not decouple from these currents but 
rather form a light-like Wilson loop with a cusp at x = [16]. The anomalous 
dimension of the combination {Sl^ Wgc) is the universal cusp anomalous dimension 
times a logarithm of the soft-coUinear scale, see the first expressions in the second 
lines in (2.69) and (2.70). After adding the contributions from the soft and collinear 
sectors, the dependence on the IR regulators drops out. However, the coefficient of 
the logarithm oiv-pc- in the heavy-collinear current andp^+'Pc- in the soft-colhnear 
current is unchanged, since both the soft and the collinear part are independent of 
these large scales. This cancellation also explains why the anomalous dimensions 
of the heavy-to-collinear and soft-to-collinear currents involve — Fcusp and — iFcusp, 
respectively: 



r 

cusp 



2 



r 

cusp 



In — — + In — + In 

{-Pi){-Pi) /i^ /i 

2v ■ Pc- , —2v ■ Ps , —pf. 

, {-2v ■ Ps) {-pi) ^ /i^ 



Tcusp In ■ 

■^Ps+ Pc- 

In 



2^""^ i2vp,.y 

(2.83) 



Similar arguments were used by Korchemsky in his analysis of the off-shell Sudakov 
form factor [31]. 

The effect of the field redefinition (2.80) on the four-quark operators is different. 
The color singlet-singlet operator is invariant, namely (setting x = for simplicity) 

Q(5)(s, t) Xisn) T^ X(0) Q,(tn) T2 ?{(0) , (2.84) 

since the additional soft-colhnear Wilson lines come in pairs W^^ Wgc = 1 and 
Sl^Ssc = 1- The color octet-octet operator is however not invariant, since in that 
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case the cancellation cannot happen. Instead, one obtains the soft-coUinear object 
[16] 

Tr [Ssc Ta Si Wsc Tb Wl] (0) . (2.85) 

The presence of this functional in the octet-octet case explains why soft-coUinear 
modes give a non-zero contribution to the anomalous dimension of the operator 
Q[0)- However, since this operator does not mix into the singlet-singlet operator 
(5(s), this effect does not propagate into physical decay amplitudes (as hadronic 
matrix elements of color-octet currents vanish). The decoupling of soft-coUinear 
fields from the color singlet-singlet operator implies that, to all orders in perturba- 
tion theory, the anomalous dimension of the four-quark operator Q{s) is the sum of 
the anomalous dimensions of the two currents X(sfi) Fi X(0) and Qs(tn) r2 3^(0), 
in accordance with the one-loop result obtained in the previous section. This ob- 
servation has important implications for applications of SCET to proofs of QCD 
factorization theorems. 

2.5 Reparameterization Invariance 

In section 2.4.2 we have briefly touched on reparameterization invariance, on which 
we expand here. Operators in SCET must be invariant under redeflnitions of the 
light-cone basis vectors n and n that leave the scaling properties of flelds and 
momenta unchanged [14, 22, 32]. This property is referred to as reparameterization 
invariance, and has proved to be a quite powerful tool. For example, it can be used 
to derive constraints on the Wilson coefficients of SCET operators, often relating 
the coefficients of some operators to those of others. Reparameterization invariance 
is a consequence of the invariance of QCD under Lorentz transformations, which is 
not explicit (but still present) in SCET because of the introduction of the light-cone 
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vectors n and n. 

Commonly one distinguishes between three classes of infinitesimal transforma- 
tions, corresponding to two different transverse boosts and a longitudinal boost: 

Type I: + e'^, invariant (with ~ A) 

Type II: ^ n'' + e''j_ , invariant (with ~ 1) (2.86) 

Type III: n^" ^n^/a, ^ an'' (with a ~ 1) 

The scaling properties for the parameters are given in parenthesis. As an exam- 
ple, note that the coUinear Wilson line Wc is invariant under type I and type III 
transformations. To derive the transformation property under type II, note that 
by definition [n ■ DcWc] = 0. The variation implies that with 6 n ■ Dc = e± ■ Dc± 
one can derive the change 6Wc from 6[n ■ DcWd = [32]. It follows that the 
correct type II transformation behaviour is Wc —>■ Wc — {n ■ Dc)^^ e± ■ Dc±Wc- 
The transformation laws for soft Wilson lines are found using arguments along the 
same line. 

It is then straightforward to show that, for example, the current operators in 
(2.64) are separately invariant under type I and type II transformations to leading 
order in A. In other words, reparameterization invariance links these operators with 
operators that appear at subleading order in A. The only non-trivial point in this 
analysis concerns the transverse coUinear derivative -D^j_, which has no n- vanishing 
variations at leading order in A under both type I and type II transformations, 

D'c^ * - ' Di:^-'-fn.Dc+0{X^) , Di:^ D':^-"^ e^.Dc^+0{X^) . (2.87) 

Note that in both cases the object t/lW'^ {ilp ciW) = ijiJ^cL is left invariant. For 
type I transformations this follows from in ■ Dc W^ = 0, whereas for type II trans- 
formations it follows since = 0. This is an important result, because without 
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the extra factor of i/i, the second operator in (2.64) would not be invariant under 
a type II reparameterization. 

The type III transformations are the most powerful ones and have non-trivial 
consequences. For example, consider the current operator in (2.65). Type III 
transformations dictates the dependence of the Wilson coefficient on the plus- 
component of the soft momentum to all orders in perturbation theory. This 
provides valuable information about the convergence of convolution integrals of 
hard-scattering kernels (Wilson coefficients) with the 5-meson light-cone distri- 
bution amplitudes, which will be an important ingredient to factorization proofs. 
With a slight abuse of notation, the soft-to-collinear current can be written in the 
following form, in which we used translation invariance to extract the dependence 
on the large energy E: 



4(0)rV>,(0) ^ / dtDit,E,fj.) [^W,]iO)T[Slq,]itn)+Oi\) . (2.88) 



Performing a variable substitution t at and a subsequent Type III transfor- 
mation (which also rescales the energy) one finds that the current operator is in- 
variant only if its Wilson coefficient obeys the homogeneity relation Di{t,E,fi) = 
a Di{at,aE, fi). Taking into account the canonical dimension of the coefficient, it 
follows that 



where the coefficient functions are dimensionless. Since the dependence of the 
Wilson coefficients on the renormalization scale /i is logarithmic, we conclude that 
to all orders of perturbation theory D{t, E, fi) ~ 1/t modulo logarithms. Therefore 
the above argument determines the behavior of the momentum-space coefficients 
on the soft momentum pi to all orders in perturbation theory. 




(2.89) 



CHAPTER 3 
STRUCTURE FUNCTIONS 



3.1 The 5-meson Light-Cone Distribution Amphtudes 

In section 2.4.3 on four-quark operators we learned a valuable lesson that touches 
the core of QCD-Factorization: Recall that color singlet- singlet operators decouple 
from the soft-collinear messenger sector, and therefore the coUinear and soft sectors 
do not interact with each other. Exclusive decay amplitudes are then expressible in 
terms of convolution integrals and matrix elements of operators of the form (2.84). 
The product Qs{tn) r2 !K(0) has a non- vanishing overlap with the B meson, while 
X{sn) Ti X(0) overlaps with the final light meson. The use of gauge-invariant 
building blocks provides the necessary (finite length) Wilson line to make the bilo- 
cal operators gauge invariant, since e.g. X(sn) Fi X(0) = [fPFc](sn) Fi [1FJ^](0) = 
^{sn) [sn, 0] Fi .^(0). The light-like separation of the fields explains why the corre- 
sponding non-perturbative wave functions after projection onto meson states are 
called "light-cone distribution amplitudes" (LCDAs). For the B meson, two such 
functions, (p^ and (f)^, arise for the two-particle Fock state configurations. They 
are defined in the HQET trace formalism as [33] 

^ {0\Qs{tn)r:K{0)\B{v)) (3.1) 



J- ^^75 



zF(/i) 



where F{fi) denotes the asymptotic value of \JMb Jb in the heavy-quark limit. 
{Mb is the 5- meson mass and fs its decay constant.) The variable u can be 
interpreted as the plus component of the spectator momentum. Often times it 
suffices to only consider these two-particle LCDAs, since currents with an extra soft 
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gluon are power suppressed. However, they become also important if the quantity 
of interest is itself power suppressed, which is the case for e.g. B light meson 
form factors at large recoil. It will then be necessary to also define wave functions 
containing an additional soft gluon or a derivative. The latter can, however, be 
eliminated using the equations of motion. Let us define the three-particle LCDAs 
of the B meson in analogy with (3.1) as [34] 



|Q,(tn) A>:^isn) T ?{(0)|i?) = - !f ^^^^ 



tr 



Since n ■ ^l^ = 0, the full light-cone decomposition of the gluon field As reads 
= •4's,± + i/^v ■ As and there are in principle two more wave functions one can 
define, namely 

{0\Q^{tn)vAsisn)T:KiO)\B) = -'-^^^^^^ 

11 + 1^ ( ~ i ~ ^"^'^^ 

tr T^^AA^{t,s) + ^A^{t,s) 

The three-particle LCDAs yli(co'i, 1^2), ^4(1^1, 1^2) are defined as the Fourier 
transforms of y4i(t, s), . . . , A^it, s), in close analogy to (3.1). 

Equation of Motion constraints. The above LCDAs are related through the 
equation of motions. In particular, the equation of motion for the light quark field 
ilp Qs = 0, which in terms of collinear fields implies Qs{x) = — ^^(x) Qs(x), 
allows us to express (p^ entirely in terms of the leading LCDA (p^ and the three- 
particle wave function Ai. It is straightforward to use (3.1) with F = 7^^ and act 
with a partial derivative with respect to the position z of the light quark field. 
Since ^ = (f^/2) n ■ d + {1^/2) n ■ d + ^± is the partial derivative in all directions, it is 
crucial to keep z slightly off the light cone, and set it to the final value z = tn only 
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after the derivative has been performed. When z is not hght-hke, the Wilson hne 
[z, 0] = PexTp[i j^dz ■ As{z)] is no longer given in terms of the product Ss{z)Sl{^), 
but rather as the product of path-ordered exponentials 



[^,0] = Ss{z)Pe^^ 



1 

da z ■ AAaz) 



Sm . (3.4) 



Expanding the Wilson line and taking a partial derivative will then lead to the 
combination^ ^s{z) — daJf^sip^z) on the left-hand side of (3.1), while the wave 
functions on the right-hand side carry additional dependence on 2^ 7^ 0, which 
needs to be taken into account when performing the derivative. The result is [34] 



Jo 

duJiduJ2S{uJ — uji — UJ2) Ai{uJi,LJ2) (3.5) 

POO 

+ I da duiduj2S{uj — iUi — aLU2) Ai{uJi,uj2) ■ 



Finally we state the explicit expression that eliminates 4>^{uj) by solving the above 
equation. Taking a derivative with respect to u leads to the desired result 



(P^{u;) = - I duj' + ! duo, duj2 -—^ A,{uj,, UJ2) 

<*oo 



(3.6) 



dull I duj2- — — Ai{uJi,uj2) 
LUUJ1UJ2 



The remaining LCDAs ^2, . . . , ^4 are not relevant to the present work. 

As mentioned earlier, many applications require only one of the above LCDAs 
at leading power, (f)^{uj,fi). Physically, the sensitivity to this universal, non- 
perturbative function arises in processes which involve hard interactions with the 
soft spectator quark in the i?-meson, for example in B ^ ^l^u, B ^ K*'~f, etc. 



^This combination coincides with a gluon field in Fock-Schwinger gauge 
4fsg(,z) = /q daaGf"'{az)n^-fu. 
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A generic factorizable decay amplitude may be written as 



A 



r—T{Eu;,fi)^^{u;,fi), 



(3.7) 



where it is assumed that the /i dependence cancels between the hard-scattering 
kernel T and the LCDA. In writing the amplitude in this way, a first step of 
scale separation was achieved, since the kernel depends on the physics associated 



whereas the LCDA is a universal, non-perturbative function which "lives" on low 
scales /i ~ Aqcd- Since there is no one scale /i at which neither of the two quantities 
contain large logarithms it is crucial to resum those large (Sudakov) logarithms 
to gain full control over the separation of physics at different scales. We thus 
have to derive and solve the renormalization-group evolution for the LCDA or, 
equivalently, the hard-scattering kernel [24, 30]. 

3.1.1 Renormalizat ion- Group Evolution 

Different operators that share the same quantum numbers can mix under renor- 
malization. It is therefore appropriate to write the relation between bare operator 
O'^^'^^uj) = J dte'"^^ Qs(tn) r J{(0) and renormalized one as 



In the case at hand the operator is, up to a Dirac trace, the product of F{fi) and 
the LCDA (j)^{u,fi). The renormalization-group equation is an integro-differential 
equation in which the anomalous dimension is convoluted with the LCDA [30]: 



with large energy scales, i.e. does not contain large logarithms for /i ~ ^E'Aqcd 





oo 




(3.9) 
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Figure 3.1: One- loop diagrams for the calculation of the anomalous dimen- 
sion of the 5-meson light-cone distribution amplitude. 



The anomalous dimension 7+ is related to the renormalization factor Z+ through 

'j+{u,u',fi) = - I duj ^' ^7^(tD, uj', jj,) - 7ir(aJ (5(cj - t^') . (3.10) 

J a m /i 

Here is the universal anomalous dimension of local heavy-light currents in 
HQET, which determines the scale dependence of F[fi). We may separate the 
on- and off-diagonal terms of 7+ and express it as 

7+(u;,ci;',/i) = Tcnspio^s)^^- + l{as) S{uj - u') + u T{uj,uj' ,as) (3.11) 

UJ 

to all orders in perturbation theory. The cusp anomalous dimension rcusp(«s) 
appears as the coefficient of the In/i term and has a geometric origin [25]: Since 
an effective heavy-quark field h{0) can be expressed as the product of a free field 
h^^\0) and a Wilson line 8^(0) extending from (—00) to along the w-direction, the 
matrix element in (3.1) contains Ss{z) 5^(0) 5*^(0) which can be combined to form 
a single Wilson line with a cusp at the origin, as shown in Fig. 3.2(a). (Recall that 
soft Wilson lines Sg extend along the n-direction.) The appearance of the single 
In fi term distinguishes the anomalous dimension of the 5-meson LCDA from the 
familiar Brodsky-Lepage kernel [29] . On the one-loop level the In fi term appears 
in the calculation of the first diagram in Fig. 3.1, where the gluon from the Wilson 
line Ss{z) 5*^(0) connects to the heavy-quark Wilson line 51,(0). 

By evaluating the diagrams in Fig. 3.1, we find the one-loop expressions [30] 
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(denoted by the superscript (-^)) to be Fculp = 4, 7*^^) = —2, and 



cusp 



(3.12) 



uj'{uli' — uj) uj{uj — Uj') 
where we have used that 7^^^^ = —3 is the one-loop coefficient of the anomalous 
dimension of heavy-to-light currents. The subscript + denotes the standard "plus 
distribution" which is defined as 



duj' 



(3.13) 



for any (test-) function f{uj'). This ensures that J du!'r{LLj,uj') = 0. 

The ffist key toward solving the RG equation (3.9) concerns the off-diagonal 
term uj T {u , u' , fi) in the anomalous dimension (3.11). We observe that 



00 

/ 



duj'ujT{uj, u', as){uj'y = u^3^{a, as) 



(3.14) 



on dimensional grounds. The dimensionless function 3^ can only depend on the (in 
general complex) exponent a. We can use a power- law ansatz 



1^0 



(3.15) 



with an arbitrary mass parameter /iq. The function / solves the RG equation 
(3.9), if the exponent a{fi) and the normalization U{a{fi),fi) obey the differential 
equations 

d 



dlnfi 
dU{a{fi),fi) 
dlnfi 



(3.16) 



-7(as) - 3'(a(/i), as) - \n — Tcnsp{as) ■ 



The ffist equation can be immediately integrated and yields a(/i) = f] + g{fi, fio) 
with initial value 77 = a(/io) and 



asi/J.) 



9(^,1^0) 



as(A»o) 



da 



rcusp(ct) • 



(3.17) 
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With this solution at hand, the second equation integrates to 

/dcx 
7(a) + 5'(?7 + 5(o(a),a) , (3.18) 

where g^{a) = g{fi, /i^), go{a) = g{fia, fJ'o), and fi^ is defined such that as{fia) = a. 
Note that g{fio, /iq) = and U{ri, /iq) = in this construction. 

The second key to the solution concerns the initial condition 4>^{uj, /io)- Defin- 
ing the Fourier transform fo(t) of the LCDA at scale /io with respect to ln(u;//io) 
allows us to express the u dependence in the desired power-law form 

cj)^(co,ij.o) = ^ J dtMt)[y^ ■ (3.19) 

— oo 

We therefore obtain an exact analytic expression for the solution of the RG equa- 
tion (3.9) as the single integral [30] 

oo 

0^(cj,/i) = ^ j dt (po{t) f{LJ,fi,fio,r] = it) . (3.20) 

— oo 

3.1.2 Asymptotic behaviour 

The solution (3.20) enables us to extract the asymptotic behaviour of the LCDA 
(j)^(uj,fi) as — * and u; — > oo by deforming the integration contour in the 
complex t plane. We hence need to study the analytic structure of the integrand 
ipo{t) f{uj,n,HQ,it) ~ ^J^^+9i^J■,^J■o) _ jf ^ jg very small we can deform the contour 
into the lower half plane, and then the position of the nearest pole to the real axis 
determines the u dependence of (f)^{uj,^). Similarly the nearest pole in the upper 
half plane dominates for very large u. 

Let us study the analytic structure of /(tu, /i, /iq, it) at leading order in RG- 
improved perturbation theory. Using the one-loop expressions (3.12) and the def- 
inition (3.14) we find 
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= + a) + - a) + 27^.] . (3.21) 

and 7^; denote the logarithmic derivative of the Euler-Gamma function and the 
Euler- const ant, respectively. Plugging this result into eq. (3.18) with rj = it one 
obtains (using the short-hand notation g = g{jji, Hq)) 

r{i - it) r{l + it + g) ^ ' 

The function f{uj,fi,fiQ,it) has poles along the imaginary axis, and the closest 
to the real axis are located at t = i and t = —i{l — g). Using the one- loop 
expression (3.17) we observe that the function g vanishes at /i = /io by definition 
and grows monotonously as fi increases. Therefore the position of the pole in the 
lower complex plane approaches the real axis under renormalization evolution, as 
illustrated in Fig. 3.2(b). These poles "compete" with the singularities arising 
from (fo(t) for the nearest position to the real axis. Let us assume that, for a given 
model of (f)^{u,fiQ), the LCDA grows like lu^ for small cu and falls off like cu^^ for 
large u. The corresponding poles of the function (poit) are then located at t = —i6 
and t = iC,- We therefore obtain the asymptotic behaviour of the renormalized 
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(a) -00 (b) 

Figure 3.2: (a) Left: The cusp in the Wilson line ^,(^) ^1(0)^^(0). (b) 
Right: Poles of the function /(..., it) in the complex t plane. 
The upper pole remains stationary under renormalization flow, 
whereas the position of the lower pole moves toward the real axis 
for increasing /x. 
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LCDA as 




,min{l, 5+g) . 



min(l,e)+g. f^^ ^ 



for UJ 



0, 



oo . 



(3.23) 



The two immediate observations are that, regardless of how small the value of S is, 
evolution effects will drive the small uj behaviour toward linear growth, and that 
the renormalized LCDA at a scale fi > fiQ will fall off slower than 1/uj irrespective 
of how fast it vanishes at /i = /iq. 

The emergence of a radiative tail after (even infinitesimally small) evolution 
seems, at first sight, a very strange property of the LCDA, because it implies that 
the normalization integral of (f)^{uj,fi) is UV divergent. This can be understood 
as the corresponding /oca/ operator Q(0) F J{(0) requires an additional subtraction 
when renormalized [33]. However, this is not an obstacle for practical applications, 
since only (j)^{uj , fi) / uj modulo logarithms appears. An integral over this function 
remains UV finite as long as g{fi, fi^) < 1, at which point the pole aX t = — g) 
reaches the real axis and the formalism presented above breaks down. 

It is evident from (3.23) that evolution effects mix different moments of the 
LCDA. For example, the first inverse moment of (j)^{uj,fi) defines a parameter 



which is connected to a fractional inverse moment of order l — g{fi, fi') at a different 
scale /i'. This makes it impossible to calculate the scale dependence of A^^u) in 
perturbation theory without knowledge of functional form of the LCDA. 

3.2 The Shape Function 

The shape function is a non-perturbative structure function that encodes the Fermi 
motion of the heavy quark inside the B meson. While the LCDA discussed in the 



oo 




(3.24) 
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last section has a rough interpretation as a probabihty distribution for the plus 
component of the spectator-quark momentum (the evolution effects do not quite 
fit into this interpretation), the shape function describes the plus component of 
the residual momentum of the heavy quark. It enters the calculation of inclusive 
B light particles decays such a.s B ^ X^l^u and B ^ X^'-f, and is defined as 
the forward matrix element of the bilocal heavy-to-heavy current 

2Mb ^ ' 2 \ 2 ) ^ ' 

The soft function S{uj) coincides with the shape function f{kj^) introduced in [35]. 
Only one function arises in the HQET trace formalism, because v is the only 
vector available by external kinematics. The shape function S{uj) has support for 
G ] — oo. A], which can be understood as follows: 
The 5-meson momentum can be decomposed into Pb = rribV + k, where rrib is 
the 6-quark mass, and k is the dynamical "residual" momentum. For the sake of the 
argument, the variable u can be thought of as the plus component n-k. Let us work 
in the rest frame and choose the coordinate system such that the three-momentum 
k points in the +z direction. One may then parameterize k'^ = uv'^ + \k\n'^, so 
that ui = n ■ k holds. The calculation of 

Ml = = {mi, + ujf + 2\k\{mb + uj) (3.26) 

leads to the endpoints of the allowed uj interval for the cases = and |A;| oo, 
which can be read off as u; = [Mb — rrib) and uj — > —rrib, respectively. In the 
heavy-quark limit oo, one therefore finds the support interval mentioned 

above, where {Mb —rn.b) ™^Z^ A is used. The quantity A ~ O(Aqcd) is an HQET 
parameter and depends on the particular definition of the heavy-quark mass. 
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t ^ ■ " 



Figure 3.3: Radiative corrections to the shape function. The bilocal HQET 
operator is denoted by the black square. A mirror copy of the 
first graph is not shown. 



3.2.1 Renormalizat ion- Group Evolution 



Renormalized shape function 

According to (3.25), the shape function is defined in terms of a hadronic matrix 
element in HQET and thus cannot be calculated perturbatively. However, the 
renormalization properties of this function can be studied using perturbation the- 
ory. To this end, we evaluate the matrix element (3.25) in HQET using external 
heavy-quark states with residual momentum k. For the time being, v ■ k kept 
non-zero to regularize infra-red singularities. The relevant one-loop graphs are de- 
picted in Fig. 3.3. Adding the tree contribution, we obtain for the matrix element 
of the bare shape-function operator 0{u)) = hV 5{u) — in-D)h (expressed in terms 
of renormalized fields) [36] 

^barc(a;) = Zy,5{uJ-n-k)-jj^V{l + e) 



(47r) 



dll 



-l-2t 



5{ijj — n ■ k) — 5{ijj — n ■ k + I) 



-|- 6{n ■ k — uj) {n ■ k — uj) ^{n ■ k — uj + 5) 



-l-e 



(3.27) 



where 5 



-2v ■ k, and 



(47r)2- 



r(2e)r(i - e) s- 



-2e 



(3.2^ 
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is the off-shell wave function renormalization constant of a heavy quark in HQET. 
The next step is to extract the ultra-violet poles from this result, which determine 
the anomalous dimension of the shape function. The renormalized shape function 
is related to the bare shape function through 

<S'rcn(t^) = / duj' Zs{uJ,Uj') Ski^T-ciuj') , With 

(3.29) 

r7 ( f\ c" / /\ ^-^^GL^ / f\ 

Zs{UJ,Ll! ) = d[uj - UJ ) -\ Zs{UJ,Ll! ) + . . . . 

An 

The result for Zs following from (3.27) must be interpreted as a distribution on 
test functions F{uj') with support on the interval — oo < u' < A. We find 

zs{uj,uj') = - + -ln^-^- ]6{uj-Lj') ' 



e A — io e J e \ lu' — lu 

(^T- (3.30) 
e2 ej e \uj' - u J ^ 

Evidently, the renormalization factor Zs depends on the parameter A setting the 
upper limit on the integration over u' in (3.29), which combines with the plus 
distribution (see (3.13)) to form a star distribution in the variable {u' — u). The 
star distributions are generalized plus distributions defined as [36, 37] 

dxFix) (-\^^ = frfx^^^^^^— ^ + F(0)ln-, 

XJ^. In X U 



<0 \"^/* ^0 



(3.31) 

'<0 \ ^ J * Jo ^ U 2 U 

where F{x) is a smooth test function. For later purposes, we note the useful 
rescaling identities 

1 \ / 1 ^ ^^^^^ 



Ax / , \ X j . \ X j . \ X j . 2 

(3.32) 
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We can now determine the renormalized shape function from (3.29). The result 
must once again be interpreted as a distribution, this time on test functions F{uj) 
integrated over a finite interval — Ahad < uj < A. (In practice, the value of Ahad is 
set by kinematics or by virtue of some experimental cut.) The result is 



5'parton(t^) = 6{uJ - U ■ k) <! 



CpOts 



'- + L ( 

24^ HAhad + n-Zc 



(3.33) 



TT 



n ■ k — UJ 
+5{n ■ k — u) In 



9(n ■ k — Lu) n ■ k — id , n ■ k — lu + 5 
^ ^ ' In h In 



2 Ahad + n ■ k , 6{n-k — uj) , 5 

/i n ■ k — UJ ^ b fi 



+ 



+ S{n ■ k — u) \n 



While it was useful to keep the heavy quark off-shell in the calculation of the ultra- 
violet renormalization factor, the limit 6 = —2v ■ A; — » can be taken in the result 
for the renormalized shape functions without leading to infra-red singularities. 
This gives 



parton 



(lu) = 5{iu — n ■ k) 



CpOLs 



TT 24 



(3.34) 



TT 



In 



n ■ k — UJ 



,n ■ k — UJ /i ) \n ■ k — UJ 

where the star distributions must now be understood as distributions in the variable 
{n ■ k — uj). We stress that these results for the renormalized shape function are 
obtained in the parton model and can in no way provide a realistic prediction for the 
functional form of S{uj). Only the dependence on the ultra-violet renormalization 
scale fi can be trusted. 



Evolution of the shape function 

The next goal is to solve the integro-differential evolution equation 
d 



dhifj, 



S{u,fi) = - / duj''ys{^^,uj',fi) S{uj',fi) . 



(3.35) 
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In the upcoming calculation of the differential decay rate in inclusive B Xu I v 



decays we will require the shape function at an intermediate scale /ij ~ lym^jAqco- 
If the shape function was known at some lower scale /xq, it would be necessary to 
evolve it according to (3.35). 

At one-loop order, the anomalous dimension for the shape function is twice the 
coefficient of the 1/e pole in the renormalization factor Zg. From (3.30), we obtain 

(3.36) 



7s(u;,u; ,/i) = 



TT 



21n^-lV(^-^') 



k — UJ I \ U)' — UJ 



The evolution equation (3.35) can be solved analytically using the aforemen- 
tioned general method presented in the previous section for the case of the 5-meson 
LCDA [30]. It is convenient to change variables from ujtouj = K — ujE [0, oo[ and 
define S{uj) = S{K — Cj). The renormalization-group equation then reads 

S{i2j,fi) = - dci)'75(a),a)',yu) ^(u)',/!) , (3.37) 



dln/i Jo 
where the anomalous dimension can be written in the general form 

7s(a),a;',/i) = 2 rcusp(Q!s) In ^ + 7(0^) 6{oj - uj') + 2Q{oj,oj',as) . (3.38) 

L UJ 

This is obviously analogous to the RGE (3.9) governing the LCDA. The loga- 
rithmic term containing the cusp anomalous dimension can again be interpreted 
geometrically. Since the heavy-quark field h{x) in HQET can be represented as 
a Wilson line along the v direction, the field !K{x) entering the SCET formalism 
contains the product of a light-like Wilson line (along n) and a time-like Wilson 
line (along v), which form a cusp at point x. The shape function contains two 
such cusps, each of which produces a contribution to the anomalous dimension 
proportional to Fcusp In [25]. The one- loop coefficients of the remaining terms in 
(3.38) are 

7o = -2Cp, go{u;,u;') = -To(^-^^^^^) . (3.39) 

\ OJ — UJ I , 
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The general solution of (3.37) can be obtained (in close analogy to our previous 
discussion) using the fact that on dimensional grounds 



(3.40) 



where the function 3^ only depends on the exponent a and the coupling constant. 
We set 5'(0, as) = by definition, thereby determining the split between the terms 
with 7 and Q in (3.38). The integral on the left-hand side is convergent as long as 
Rea > — 1. At one-loop order we find from (3.39) 



a. 



3'{a,as) = To— '4j{l + a)+'jE 



(3.41) 



where ip{z) is the logarithmic derivative of the Euler F function. Relation (3.40) 
implies that the ansatz 

r+23(M,/io) 



/(^,/W,/io,r) 



expUs{T,fi,ij,o) 



(3.42) 



with 



9{fJ',fJ'o) 



Us{t,h,Ho) 



da .^^-^P^") 



"s(/io) 



Pia) ' 



(3.43) 



o's{^J.o) 



da 

pW) 



g{n, Ha) + 7(«) + 3'{t + 2g{na, /^o), «) 



provides a solution to the evolution equation (3.37) with initial condition 
f{Lj, /io, yUO) t) = i'^/fJ'o)'^ at some scale /xq- Here /i^ is defined such that as{fia) = 
and r can be an arbitrary complex parameter. Note that g{fi, fio) > if /i > fiQ. 
We now assume that the shape function S'(a),/io) is given at the low scale /xq and 
define its Fourier transform with respect to ln((i//io) through 



S{LJ,fIo) 



2n 



dtSoit) (— 



(3.44) 



74 



The exact result for the shape function at a different scale /i is then given by 

r"00 



'S'P, yu) = — / dt So{t) f{u, fi, fio, it) 



(3.45) 



With the help of this formula, it is straightforward to derive explicit expressions 
for the evolution of the shape function from the hadronic scale /io up to the in- 
termediate scale /ij at any order in renormalization-group improved perturbation 
theory. Setting r2 = as{fJ'o)/as{f^i) > 1, we obtain for the evolution function at 
leading order 

T{l + it) 



r{l + it + ^ Inrs) 



(3.46) 



where 



— - — - (r2 — 1 — In r2) + \v? r2 



In r2 

r2 



(3.47) 



To , 7o , , ^ 

--rlE lnr2 - — lnr2 + O 



(r2 - 1) asifi) 



This result is valid as long as (Fo/Z^o) liir2 < 1, which is the case for all reason- 
able parameter values. Missing for a resummation at next-to-leading order are the 
0{as) contributions to Vs, which vanish for fi fiQ. For all practical purposes, 
given the intrinsic uncertainties in our knowledge of the shape function, it will 
be sufficient to use the equations given above. As mentioned earlier, we typically 
have Hi ~ mc, and so the running between /ij and /io should be performed in a 
theory with nj = 3 light quark flavors. The relevant expansion coefficients are 
then Fo = f , Fi = ^ - f vr^, 70 = -|, and /Sq = 9, /3i = 64. 

The leading-order result presented above can be simplifled further. When (3.46) 
is inserted into (3.45), the integration over t can be performed analytically. Setting 



75 



r] = (Tq/Po) \nr2 > 0, the relevant integral is 



where 

5o(t)= / ^^(c^',/io) - (3.49) 

is the Fourier transform of the shape function as defined in (3.44). The integrand 
of the t-integral has poles on the positive imaginary axis located at t = in with 
n > 1 an integer. For uj < uj' the integration contour can be closed in the lower 
half-plane avoiding all poles, hence yielding zero. For uj > iIj' we use the theorem 
of residues to obtain 

1=1 duj'R{uj,uj')S{uj',fio), (3.50) 



^0 

where 



(3.51) 



j=0 

Note that R{u,(l!') 6{uj — u') in the limit r/ — 0, corresponding to /ij — > /io, as 
it should be. 

Our final result for the shape function at the intermediate hard-collinear scale, 
valid at leading order in renormalization-group improved perturbation theory, can 
now be written in the simple form [36, 38, 39] (valid for /ij > /io, so that ?7 > 0) 

r(^) io 

with Vs as given in (3.47). 

From the above equation one can derive scaling relations for the asymptotic 
behavior of the shape function for ci) ^ and a) ^ oo (corresponding to ^ A 
and UJ — > — oo). If the function S{uj,fio) at the low scale fiQ vanishes proportional 
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to Cj'^ near the endpoint, the shape function at a higher scale /ij > /i vanishes faster, 
proportional to Similarly, if S{u,fio) falls off like u)~^ for u) — ^ oo, the shape 

function renormalized at a higher scale vanishes like c2;~™™^^'^^+''. Irrespective of 
the initial behavior of the shape function, evolution effects generate a radiative tail 
that falls off slower than This fact implies that the normalization integral 

of S{lj, fi) as well as all positive moments are ultra-violet divergent. The field- 
theoretic reason is that the bilocal shape-function operator contains ultra-violet 
singularities as 2;_ — 0, which are not subtracted in the renormalization of the 
shape function. The situation is analogous to the case of the 5-meson light-cone 
distribution amplitude discussed in [30, 33], see Section 3.1.2. These divergences 
are never an obstacle in practice. Convolution integrals with the shape function 
are always cut off at some finite value of lj by virtue of phase-space or some 
experimental cut. 

3.2.2 Properties of the shape function 

Information about the shape function can be extracted from a study of its moments 
using a local operator product expansion. Naively, one would define the moments 
Mn = j^^duJUJ^S{u). Prior to our work in [36], the general understanding was 
that, at tree level, Mq = 1 is fixed by the normalization of the shape function. 
Ml = by virtue of the HQET equation of motion iv ■ D h = (which implic- 
itly uses the pole-mass definition of rrib), and M2 = — Ai/3 is given by the matrix 
element of the kinetic-energy operator. However, it was not clear how to system- 
atically include radiative corrections to these relations. Apart from the obvious 
dependence on HQET parameters such as A and Ai, the moments Mn must also 
depend on the renormalization scale fi. However, as was studied in the last Sec- 
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tion, quantum corrections result in the appearance of a radiative tail, rendering 
all moments with > ultra-violet divergent. Since in any physical process the 
integration over the shape function is restricted to finite intervals, it suffices to 
define the moments with an ultra-violet cutoff, such that 



MAr(Auv,/i)= / duJLU S{iU,fl) 
-Auv 



(3.53) 



In the following Subsections we will expand the finite quantities M7v(Auv,/u) in a 
local HQET operator product expansion, assuming that Auv is large compared to 
Aqcd- The only relevant operators up to dimension 5 are of the form h [in ■ D)™ h 
with m = 0,1,2, whose matrix elements are given as 1,0,— Ai/3, respectively 
[35, 40, 41, 42]. 



Moments in the pole scheme 

We can write an expansion of the form 

M^(Auv,/i) = A{^v <( ^(Auv,/i) + ^(Auv,/i) 



3A2 



O 



Aqcd Y 



Auv; Jj' 

(3.54) 

where the coefficients K^^"* can be determined from a perturbative matching cal- 
culation using on-shell external 6-quark states with residual momentum k. By 
evaluating the moments of the renormalized shape function in (3.34), one finds at 
one-loop order [36] 



Mr°'^(Auv,/i) = (n-A;) 



N ) Cf^s |^j^2 Auv + n-k ^ Auv + n-k ^ vr^ ^ 



vr 



/i 



24 J 



CpOis 1 I ]^ _^ 2 In + n ■ k _ 2 



A 



uv 



n ■ k 



1 



(3.55) 



We then expand this result in powers of n ■ k/Ayy. Keeping the first three terms 
in the expansion yields 



78 



parton / » \ 

lAuv,/^J 
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TT \ /i /i 24^ 
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Mr*°"(Auv,/i)=n-A; 



CpOis 



TT 



C^Q^g ^UV 
TT 

Auv I -2 In 



A2 



/i 



Auv 



Mr*°'^(Auv,/x) = (r^-A:)2 



1 - 



TT 



In- 



/i 

2 Auv 

/i 



/i 24 

Auv At 

^2 1 



-2 In 



TT 



Auv 
/i 24 
Auv 



A^v In ^ + n ■ A; Auv ( -2 In 



(3.56) 



On the other hand, we need to calculate the one-loop expressions for the local 
operators h {in ■ D)^ h. The relevant diagrams are the same as in Fig. 3.3, where 
now the black square represents the local operators. The result is non-trivial 
when keeping v ■ k non-zero to regularize infra-red singularities. However, for the 
matching calculation we need the limit v ■ k —* 0, which can be taken without 
problems in the sum of all diagrams (but not for each individual diagram). In that 
case, the one-loop corrections vanish and the matrix elements reduce simply to their 
tree-level values. It follows that in (3.56) we must identify (n-k)"' — ^ {h {in-D)"'h). 
Substituting the results for the HQET matrix elements given earlier, we obtain for 
the Wilson coefficients of the first three moments 



K, 
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Cpas 



TT 



In 
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/i 24 
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(1) _ CpOis 
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TT 
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21n- 
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uv 



In 



A^ 

Auv 



K. 
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TT 





_ Cpas A Auv 




TT \ fi 




Cpots , Auv 
= —2 m 




TT fl 


Auv 




^^ 


fi 24 



2 

(3.57) 



We observe that this result reduces to the naive moment relations mentioned above 
in the tree-level approximation. The perturbative quantum corrections can be 
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trusted as long as the ratio Auv//^ is of 0(1). 
The shape-function scheme 

We have stressed before that the value of the moments depend on the definition of 
the heavy-quark mass. So far, our calculations have assumed the definition of the 
heavy-quark mass as a pole mass, which is implied by the HQET equation 

of motion iv ■ D h = 0. Results such as (3.57) are valid in this particular scheme. 
A more general choice is to allow for a residual mass term 6m in HQET, such that 
iv ■ D h = 6m h with 6m = O(Aqcd) [43]. It is well known that the pole mass is an 
ill-defined concept, which suffers from infra-red renormalon ambiguities [44, 45]. 
The parameter Apoie = Mb — '^-1°^'^, which determines the support of the shape 
function in the pole-mass scheme, inherits the same ambiguities. It is therefore 
advantageous to eliminate the pole mass in favor of some short-distance mass. 
For the analysis of inclusive 5-meson decays, a proper choice is to use a so-called 
low-scale subtracted heavy-quark mass mf,(/i/) [46], which is obtained from the 
pole mass by removing a long-distance contribution proportional to a subtraction 
scale /i/, writing m^"^'^ = mb{fif) + fij g{as{^)-, = m^^fif) + 6m. As long as 
mb{fif) is defined in a physical way, the resulting perturbative expressions after 
elimination of the pole mass are well-behaved and not plagued by renormalon 
ambiguities. Replacing the pole mass by the physical mass shifts the values of 
n ■ k and u by an amount 6m, since n ■ (m^°'°f + k) = mi,{fif) + {n ■ k + 6m), and 
because the covariant derivative in the definition of the shape function in (3.25) 
must be replaced hj in ■ D — 6m [43]. At the same time, Apoic = A(/ij) — 6m, 
where A(/ij) = Mb — mi,{fif) is a physical parameter. Note that this leaves the 
parameter u = A — u and hence the shape function S{u, fi) invariant! This follows 
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since Apoic — cUpoie = — (cjpoic + Sm), where Upoic denotes the value in the 

pole-mass scheme used so far. We now choose 6m such that the first moment Mi 
vanishes, thereby defining a low-scale subtracted heavy-quark mass to all orders in 
perturbation theory (with /i/ = Auv), called the "shape-function mass" mf^ [36]. 
The shape-function mass can be related to any other short-distance mass using 
perturbation theory. To this end, one uses the fact that (3.56) implies a relation 
to the pole mass 

2 i-Xi] 



pole SF/ \ , Cpasifi) 

TT 



2 In 



ln^ + ... 

/i 



(3.5^ 



Given the above expression, it is easy to relate also to the potential-subtracted 
mass introduced in [47] and to the kinetic mass defined in [48, 49]. 



PS 



kin / 



3tt 



(3.59) 



After the introduction of the shape-function mass, the coefficients Kn^ in (3.57) 
vanish by definition. At one-loop order, the remaining coefficients of the zeroth and 
second moment remain unchanged since 6m = 0{as)- Proceeding in an analogous 
way, we can use the second moment to define a physical kinetic-energy parameter, 
commonly called /i^. This quantity can be used to replace the HQET parameter 
Ai, which like the pole mass suffers from infra-red renormalon ambiguities [50]. At 
one-loop order, we obtain [36] 

/i2(Auv,^) _Mf^^(Auv,/x) 



-Ai) 



TT /i 



CpasifJ') Auv 1 
TT \ ^ 2 



+ ... . 

(3.60) 



This definition is similar to the running parameter n^. defined in the kinetic scheme 
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[48, 49]. At one-loop order, the two parameters are related by 



2tt 



(3.61) 



Given a value for the kinetic energy in the shape- function scheme for some choice 
of scales, we can solve (3.58) and (3.60) to obtain values for mf^ and /i^ at any 
scale, using the fact that and Ai are scale independent. 

3.2.3 Moments of the scheme-independent function S{uj,ii) 

The variable a; = A — u; and with it the function S{i2}, fi) = S{uj) are independent 
under redefinition of the heavy-quark mass. It will be useful to rewrite the moment 
relations derived above in terms of these quantities, defining a new set of moments 

Miv(/i/,/i)= j dtou^ S{u,ix) . (3.62) 







This yields 



HU,,,,) = l-^f^(,n'tl + i,,t^ + Q (3.63) 

TT \ fi fi 24: J 

TT V At 2 J 3yUj 

^Si^ = A(/.„;.), M2(/i/,/i) ^ /i^(/i/,/i) ^^( 

where the parameters A(/ij,/i) = Mb — mf^{fLf,fL) and fiKfif,^) should be con- 
sidered as known physical quantities. The moments with > 1 give simply a 
restatement of the shape-function scheme definition. Another interesting aspect is 
that the expression for Mo{fif,fi) can be used to extract the precise form of the 
asymptotic tail of the shape function, since /i/ is considered to be much larger than 
Aqcd in our calculation. 
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Asymptotic behavior of the shape function 



Taking the derivative of the zeroth moment Mq in (3.62) with respect to /i/, one 
obtains 



dmf{fif,fi) 



From (3.63) we find at one-loop order 
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(3.65) 
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The precise definitions of A and /i^ are not specified at this order. (Note that the 
shape function cannot depend on the value of the cutoff fif.) Relation (3.65) is 
a model- independent result as long as a) ^ Aqcd- We stress the remarkable fact 
that this radiative tail of the shape function is negative, in contrast with the naive 
expectation based on a probabilistic interpretation of the shape function as a mo- 
mentum distribution function. The point is that the definition of the renormalized 
shape function requires scheme-dependent ultra-violet subtractions. From (3.65) 
it follows that the shape function must have a zero, which for sufficiently large fi 
is located at a value ujq ^ A + iij 



CHAPTER 4 
INCLUSIVE SEMILEPTONIC DECAYS 



Due to experimental cuts in the measurement of the inclusive B Xul~u 
decays one is generally faced with a situation in which the hadronic final state 
Xu is constraint to have large energy Eh ~ Mb, but only moderate invariant 
mass sh = rn\- ~ Aqcd Mb- This kinematic region, called the "shape-function re- 
gion" , is the dominant phase space in which the final state cannot contain charmed 
hadrons. The calculation of the inclusive differential decay rate uses the optical 
theorem and assumes quark-hadron duality. This assumption is justified when in- 
tegrating over large portions of the phase space, corresponding to the summation 
over many final hadronic states. To calculate the total decay rate, i.e. including 
all of phase space, all kinematic quantities are integrated over a domain of order 
Mb, and one can perform an operator product expansion (OPE) to systematically 
compute power corrections. However, as mentioned above, when it is necessary to 
calculate the differential decay rate in the shape-function region, the calculation 
becomes more difficult because of the presence of three separated mass scales: the 
hard scale Mb, the hard-coUinear scale ^/Aqcd Mb, and the soft scale Aqcd- To 
properly disentangle the physics associated with these scales requires a sophisti- 
cated effective field-theory machinery. 

A systematic treatment consists of matching QCD onto SCET/ in a first step, 
in which hard quantum fluctuations are integrated out. The degrees of freedom 
in SCET/ are referred to as soft and hard-collinear to indicate that the invariant 
final state momenta have fluctuations of order the intermediate scale -^Aqcd Mb, 
much larger than Aqcd- It is therefore possible to treat them perturbatively. The 
expansion parameter of SCET/ is hence defined as A = Aqcd/ Eh- 
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In a second step, SCET/ is matched onto HQET, and hard-coUinear modes 
are integrated out. The resulting expressions for inclusive differential decay rates 
have the factorized form dT ~ H J S [51]. The function H contains the hard 
corrections, the jet function J, which describes the properties of the final-state 
hadronic jet, contains the hard-coUinear effects, and the shape function S accounts 
for the internal soft dynamics in the B meson [35, 52]. The ® symbol implies a 
convolution over a light-cone momentum variable u associated with the residual 
momentum of the b quark inside the B meson. 

4.1 Factorization theorem 

Using the optical theorem, the hadronic physics relevant to the inclusive semilep- 
tonic decay B Xy_l^u can be related to a hadronic tensor W^'^ defined via 
the discontinuity of the forward 5-meson matrix element of a correlator of two 
flavor- changing weak currents = n7^(l — 75)6 [53, 54, 55]. We define 

W^" = -Im (-^(^)l^n-^(^O) ^ T^'^ = t [ d^xe''-'' T{J^^{0),r{x)}. (4.1) 
TT 2Mb J 

Here v is the 5-meson velocity and q the momentum carried by the lepton pair. 

After the field redefinition b{x) = e"*"^''^'^ which is always the first step 

in the construction of an effective heavy field, the phase factor in (4.1) becomes 

^i{q-mtv)-x ^ e~*^'^, where p = m^v — q corresponds to the momentum of the jet 

of light partons into which the 6-quark decays. We will assume that these partons 

can be described by hard-coUinear fields. This is justified in the shape-function 

region, in which case the current correlator can be expanded in non-local light-cone 

operators [35, 40, 52]. 

The hadronic tensor in (4.1) factorizes at leading power [36], as can be shown 

in the following way: The jet momentum (and likewise the momentum of the 
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hadronic final state) scales like ~ E{X, 1, ^/X). (For the jet momentum, p_L = 
by choice of the coordinate system.) The jet invariant mass, ~ EAqcd, defines a 
hybrid, intermediate short-distance scale. The appropriate effective field theory for 
integrating out the short-distance fluctuations associated with the hard scale p_ is 
SCETj, see Section 2.3 but with hard-coUinear and soft degrees of freedom instead 
of coUinear and ultrasoft ones. Below a matching scale fih ~ rrii,, the semileptonic 
current can be expanded as 

u{x)ril - l^)h'{x) = I ds a(s) XU^ + sn) Tf J{(x_) + . . . , (4.2) 

i=i 

where the dots denote higher-order terms in the SCET expansion, which can be ne- 
glected at leading power in Aqcd/"16. The hard-collinear light-quark field Xhc{,x) = 
Sl{x^)Wl^{x) C,hcix) and the soft heavy-quark field 3-C(x_) = Sl{x^) h{x^) are 
SCET building blocks that are invariant under a set of homogeneous soft and 
hard-collinear gauge transformations [13, 14, 15]. These are the ingredients of 
SCET/ after the field redefinitions (2.38) and (2.40). Because of the reduced Dirac 
basis between coUinear spinors in (4.2), we may choose any three independent 
Dirac structures that do not vanish between the spinors. A convenient choice is 

15^ = 7^^(1-75), r^ = ^^(l + 75), r^ = ^(l + 75). (4.3) 
The current correlator in (4.1) then becomes 

T'"' = t I d^xe-'P-"" J2 I dsdtC*{t) Ci{s) x 

i,j=i (4.4) 

T { J{(0) Xhcitn), Xhc{x + sn) J{(x_)} + . . . . 
In a second step, the hard-collinear fluctuations associated with the light-quark 
jet can be integrated out by matching SCET onto HQET at an intermediate scale 
A^i ~ i^mbAQCD- At leading order the SCET Lagrangian (when written in terms of 
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the gauge-invariant fields such as Xhc) does not contain interactions between hard- 
colhnear and soft fields, due to the decoupling transformation (2.38). Since the 
external 5-meson states only contain soft constituents, we can take the vacuum 
matrix element over the hard-collinear fields, defining a jet function 

{n\T{XhcAi^),Xhc,i{^ + sn)}\n)=5kiJ{x + {s-t)n) + ... , (4.5) 

where k, I are color indices, and we have used translational invariance to determine 
the dependence on the coordinate vectors. Shifting the integration variable from x 
to z = x + {s — t)n, with = and introducing the Fourier-transformed Wilson 
coefficient functions Ci{n ■ p) = J ds e^'^^'^ Ci{s), we then obtain 

T^- = z ^ C*{n ■ p) Ci{n ■ p) / (fze-'^'' :H(0) f^^ J{z) J{(^_) + . . . . (4.6) 
In the next step, we rewrite the bifocal heavy-quark operator as [35] 

'K{Q)T'K{z_) = {hSs){0)Te'--^+ {Slh){0) = h{0)re'--^+ h{0) 

= J (kue--2''''-'h{0)T5{cu-in- D)h{0), (4.7) 

where F may be an arbitrary Dirac structure, and we have used the property 
in ■ D Ss = Ssin ■ d oi the soft Wilson line Ss- When this expression is used in 
(4.6), the resulting formula for the correlator involves the Fourier transform of the 
jet function, 

j d'ze-''^-'J{z) =^_J{p^), (4.8) 

however with p^ replaced by the combination pf^ = p^ + Intuitively, this 

happens as only the plus component of the soft residual momentum of the b quark 
adds to the hard-collinear momentum p. Using that p^^^ = P-, we now obtain 

3 

T'^-' = iJ2Hijin-p) duoJ{pl)hf^^i,_T]5{u-in-D)h + ..., (4.9) 
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where Hij{n ■ p) = C*{n ■ p) Ci{n ■ p) are called the hard functions. 

We need the discontinuity of the jet function, J(p^) = ^\ni[ij(jp')], in order 
to compute the hadronic tensor. The S-meson matrix element of the soft operator 
is evaluated using the HQET trace formalism, which allows us to write [2] 



2Mb ' ' 2 \ 2 ^ 

at leading power in the heavy-quark expansion. Clearly, the soft function S{uj) is 
the shape function, which was studied in Section 3.2. This gives the factorization 
formula 

W'"' = H.A^ ■ P) (f r t i±i^ J dco J{pI) Sico) + .... (4.11) 

In the final expressions (4.9) and (4.11) the dependence on the three scales 
n-p nib, vt ~ ■'^feAqcD and u ~ Aqcd has been factorized into the hard, jet, and 
shape functions, respectively. The factorization formula (4.11) was derived at tree 
level in [35, 52], and was conjectured to hold to all orders in perturbation theory 
in [51]. The derivation presented above [36] is equivalent to an all-order proof of 
this formula first presented in [11] (see also [8]). The limits of integration in the 
convolution integral are determined by the facts that the jet function has support 
for p5 — O5 ^^'i shape function has support for — 00 < u; < A. The argument 
pIj of the jet function can be rewritten in the "partonic variable" n ■ p and the 
"hadronic variable" n ■ Ph ~ O(Aqcd): 

pIj = p^ + n ■ puj = n ■ p {n ■ Ph — (K — oj)) = n ■ p {n ■ P^ — Cj) , (4-12) 

where Ph = Mb v — q = p + Av is the 4-momentum of the hadronic final state, and 
the variable a) = A—u > 0. Finally, n-Pn = Eh — \Ph\ = SHf^Efj+O^AQf^j^/mh) is 
related to the hadronic invariant mass and energy of the final state. The usefulness 
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of this variable has also been emphasized in [38, 56]. We shall see below that 
expressing the convolution integral in terms of the new variable Cj eliminates any 
spurious dependence of the decay spectra on the 6-quark pole mass. 

Using the fact that the Wilson coefficients Cj are real and hence Hij is sym- 
metric in its indices, we find 

i,j=i ^ 

+ 2H,2 V ■ p. v^v'' + 2{H,2 + if23) {p^-v"" + P^'.v^) + 2(2/fi3 + H^z) , (4.13) 

V ■ p- 

which may be compared with the general Lorentz decomposition of the hadronic 
tensor given in [37]: 

W^"" = Wi [p^v" + p'^v^' -g^'vp- ie^"""^p^Vf3) - W2 g"" 

+ W3 v^v'' + {p^v'' + p'^v") + W^p'^p'^ (4.14) 

We see that the structure function W2 is not generated at leading order in the 
SCET expansion. Since only the Wilson coefficient Ci is non-zero at tree-level, the 
structure function Wi receives leading-power contributions at tree level, whereas 
W4 and 14^5 receive leading-power contributions at 0{as{mb)). The function Ws 
receives leading-power contributions only at 0(a^(mb)), which is beyond the ac- 
curacy of a next-to-leading order calculation. 

4.2 Matching calculations 

We derive expressions for the perturbative functions Hij{n ■ p) and J {pi,), that 
enter the factorization formula in (4.11) at next-to-leading order in a^. To this 
end, we match expressions for the hadronic tensor obtained in full QCD, SCET, 
and HQET, using for simplicity on-shell external 6-quark states. 
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4.2.1 Hard functions 

Perturbative expressions for the hadronic functions Wi in the decomposition (4.14) 
have been obtained in [37] by evaluating one-loop Feynman graphs for the current 
correlator T^'^ using on-shell external quark states with residual momentum k 
(satisfying v ■ k = 0) in full QCD. The leading terms in the region of hard-coUinear 
jet momenta are 
1 



2^ 



2 ^1 = ^(Pk 



' Sln^y- 101n?/ + - + AL2{l-y)' 



47r 



47r V l-y 3 

+ ... , 



1 u/ 2^ C'^fts 2 ylny 



2 Air 1-y \l-y 



(4.15) 

whereas and do not receive leading-power contributions at this order. Here 
as = CisifJ'), y = n ■ p/rrib, and pi = p^ + n ■ pn ■ k. The star distributions are 
generalized plus distributions, which also play an important role in the renormal- 
ization of the shape function. They have been introduced in (3.31) and obey the 
rescaling identities (3.32). In order to find the hard functions Hij, we need the 
discontinuity of the current correlator (4.4) between on-shell heavy-quark states in 
SCET, keeping i,j fixed. The corresponding tree diagram yields 

D'^''^ = K 6{pI) , with K = Ub{v)T'^^^r'^ui,{v), (4.16) 

where ui){v) are on-shell HQET spinors normalized to unity, and the quantity 
K corresponds to the Dirac trace in (4.11). The interpretation of this result in 
terms of hard, jet, and soft functions is that, at tree level, J^^\p'i) = S{p'^) and 
'S'parton(^) = ^(^ — n ■ k) in the free-quark decay picture. It follows that the 
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Figure 4.1: One- loop diagrams contributing to the current correlator in 
SCET. The effective current operators are denoted by crossed 
circles, and hard-collinear propagators are drawn as dashed lines. 
Mirror graphs obtained by exchanging the two currents are not 
shown. 



convolution integral J duj J(p'^) S{u!) in (4.11) produces S{p1), and comparison with 
(4.15) shows that = 1, while all other hard functions vanish at tree level. 

At one-loop order we need to evaluate the SCET diagrams shown in Fig. 4.1. 
The first three graphs contain hard-collinear gluon exchanges, while the last two 
diagrams contain soft exchanges. For the sum of all hard-collinear exchange graphs, 
we find [36] 
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(4.17) 
(4.18) 



The sum of the soft contributions is given by 
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(4.19) 



where L = ln{n ■ p/fi). The 1/e poles in the sum of the hard-collinear and soft 
contributions are subtracted by a multiplicative renormalization factor Zj applied 
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to the bare current correlator in (4.4), where 

is the (momentum-space) current renormahzation constant in SCET [9]. Since the 
wave-function renormahzation factor of on-sheU heavy quarks is equal to 1, the 
sum of (4.17) and (4.19) after subtraction of the pole terms is matched with the 
results in (4.15), so that the hard functions equal 

Hn{n-p) = l + ^(^-AL' + 10L-A\ny-^-AL,il-y)-^-12 

Hi2{n ■ p) = — — + 1 , 

47r 1-y \l-y J 

rr I- ^ CpOs V ( I - 2y \ 

Hi3{n-p) = — — Iny - 1 . 

A7rl-y\l-y J 



(4.21) 



4.2.2 Jet function 



The SCET loop graphs in Fig. 4.1 determine the one-loop contributions to the 
product of the jet function and the shape function in (4.11). We may write this 
product symbolically as J^^^ ® S^^^ + J^^^ ® S^^\ where the ® symbol means a 
convolution in uj. Although the realistic shape function is a hadronic quantity, 
we may use its partonic version to determine the perturbative jet function. Its 
one-loop contribution J^^^ must therefore be extracted from the results (4.17) and 
(4.19). To this end, we need the renormalized shape function at one-loop order 
in the parton model, which has been done in Section 3.2.1. It follows that the 
renormalized jet function is given by the distribution [36] 



2\ jr/ 2\ . ^F^Cs 



(4.22) 
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It will often be useful to separate the dependence on n ■ p and n ■ Pu in this result 
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by means of the substitution = yp^, where = mb{n ■ Ph — u}) according to 
(4.12). Using the identities (3.32), we find 



2\ 9/ "2 \ c-/-2\ , C^pOis 



(4.23) 
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4.3 Sudakov resummation 

Equations (4.21) and (4.22) determine the short-distance objects Hij and J in the 
factorization formula (4.11) at one- loop order in perturbation theory. However, 
there is no common choice of the renormalization scale fi that would eliminate 
all large logarithms from these results. Likewise, the shape function, being a 
hadronic matrix element, is naturally renormalized at some low scale, whereas 
the short-distance objects contain physics at higher scales. The problem of large 
logarithms arising from the presence of disparate mass scales can be dealt with 
using renormalization-group equations. Proceeding in three steps, our strategy 
will be as follows: 

1. We match QCD onto SCET and extract matching conditions for the hard 
functions Hij at a high scale fih ~ f^ib- At that scale, no large logarithms 
appear and so the hard functions can be reliably computed using perturbation 
theory. We then evolve them down to an intermediate hard-coUinear scale 
l^i ~ A/^^foAqcD by solving the renormalization-group equation 

— ^ Hij{n ■ /i) = 27j(n ■ p, /i) Hij{n ■ p, /i) , (4.24) 
a In /i 

where 7j = is the anomalous dimension of the semileptonic heavy-to- 
coUinear current in SCET, see the first line in (2.67). 
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2. Starting from a model for the shape function S{uj,fio) at some low scale 
fio = few X Aqcd large enough to trust perturbation theory. Such a model 
could be provided by a QCD-inspired approach such as QCD sum rules or 
lattice QCD, or it could be tuned to experimental data. We then solve the 
integro-differential evolution equation (3.35) to obtain the shape function at 
the intermediate scale /Xj. 

3. Lastly we combine the results for the hard functions and for the shape func- 
tion with the jet function J in (4.22) at the scale /ij, where J is free of large 
logarithms and so has a reliable perturbative expansion. The dependence 
on the matching scales fih and /ij cancels in the final result (to the order at 
which we are working). 

We begin with the evolution of the hard functions. The anomalous dimension 7j 
for the SCET current is twice the coefficient of the 1/e pole in the renormalization 
factor Zj in (4.20). More generally [9, 24], 

jj{n ■p,fi) = -r,,,p(a,) In ^ + 7' (a.) = ^ f - In ^ - ^"j + . . . , (4.25) 

where T^usp = Cf^s/'^ + • . • is the universal cusp anomalous dimension governing 
the ultra-violet singularities of Wilson lines with light-like segments [25] . The exact 
solution to the evolution equation (4.24) can be written as [36] 

Hij{n ■ p, /ii) = Hij{n ■ p, fit) exp Unin ■ p, fih, fii) , (4.26) 

where 



Unin-p, ij^h, iJ.i) = 2 

as{^J.h 



da 



TL ' p f da^ 
r™sp(a)|ln^- j ^|+7'(«) 



. (4.27) 



(3{a 

) 

Setting ri = as(/ij)/as(yU/i) > 1, and expanding the evolution function to 0{a 
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we obtain 



- -g^ In ri 



(4.28) 
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(4.29) 



(As a reminder, the QCD /^-function is given in (2.18), and we expanded all anoma- 
lous dimensions as in (2.24).) The terms proportional to l/as{fih) resum the lead- 
ing, double logarithmic terms to all orders in perturbation theory. The remaining 
0(1) terms in Vh contribute at leading, single-logarithmic order. At next-to- 
leading order, the corrections proportional to the coupling agifJ^h) are included. In 
our case, the only piece missing for a complete resummation at next-to-leading 
order is the 0{as) contribution to Vh, 
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(4.30) 



which is independent of the kinematic variable n • p and vanishes for fii ^ fih. To 
compute these terms would require to calculate the cusp anomalous dimension to 
three loops (knowledge of r2) and the anomalous dimension 7' to two loops (7^). 
While the former has recently been computed in [57], the latter is still missing 
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up to date. This implies a universal, process-independent small uncertainty in 
the normalization of inclusive 5-decay spectra in the shape-function region. We 
stress, however, that this uncertainty cancels in all ratios of decay distributions, 
even between B — > l^v and B X^'y spectra. 

It is appropriate to perform the running between fi^ and fii in a theory with 
rif = 4 light quark flavors, since the intermediate scale /ij will be of order in 
our applications below. The relevant expansion coefficients far as they are 

known, 



16 „ 2576 16 ^ „ 96488 5152 , 176 , 160^ 

"3 



3 27 3 81 27 15 9 ^' 

25 ^ _ 154 21943 

7^ = -f • (4.31) 



4.4 Differential decay rates and spectra 

As mentioned earlier, the hadronic tensor is most naturally expressed in terms of 
the variables n ■ Ph and n-p in the shape-function region. It is thus useful to derive 
expressions for the decay rates in terms of these variables. Our theoretical results 
are valid as long as n ■ Ph can be considered as being of order a hadronic scale 
(say, a few x Aqcd), whereas n-p is integrated over a domain of order nib ^ Aqcd- 
It is this integration which provides a sampling over sufficiently many hadronic 
final states needed to ensure quark-hadron duality [58] . The distribution in n ■ p 
will be described in terms of a "partonic" scaling variable y = n ■ p/nih, while the 
distribution in the orthogonal light-cone component is described in terms of the 
dimensionfull hadronic variable P+ = n ■ Ph = Eh — \Ph\- At leading order in 
^QCT)/^b, we obtain from [37] the triple differential decay rate 
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dx dy dPj, 



{l + x-y) — + x 2 



+ ... , 

(4.32) 



where x = 1 — x, and x = 2Ei/mb is a scaling variable proportional to the en- 
ergy of the charged lepton measured in the 5-meson rest frame. The quantity 
Ftree = ^1^1 V„fep(m^°^'^)^/(1927r^) denotes the leading power, tree-level expression 
for the total B — > I'v decay rate. The hadronic function Wi are written in the 
factorized form, which are found to be [36] 
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(4.33) 



where the dots represent power corrections in Aqcd/"^;)- The result for the rescaled 
jet function J{p'^,y, fJ'i) is the expression in (4.23). Furthermore we need the Su- 
dakov exponent Vh in (4.29). This function is independent of the kinematic vari- 
ables y and p^. In addition, we need 
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(4.34) 
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For simplicity, we have identified the high-energy matching scale /i/^ with the heavy- 
quark mass rub. In the variables x, and P+, the phase space is remarkably simple: 



< P+ < Mb -2Ei = rribX + a, 



rrib 



<x<y<l 



(4.35) 



If El is integrated over a domain of order nib ^ Aqcd? i- e. x is integrated over a 
domain of order unity, one can replace the second condition byO<x<i/<lat 
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Figure 4.2: Hadronic phase space for the hght-cone variables P_ and P+ 
(left), and theory phase space for rub = 4.8 GeV (right). The 
scatter points indicate the distribution of events as predicted by 
the model of [37] . In each plot the solid line separates the regions 
where sh < M'jj (dark gray) and sh > M"^ (light gray), whereas 
the dashed line corresponds to P+ = M^/Mb- The dotted line 
in the first plot shows the contour where = {Mb — Mo)^. 



leading power in Aqcb/ i^b- If, on the other hand, the lepton energy is restricted to 
be close to its kinematic limit, Ei ^ Mb/2, then x = O^AqcD/^b), and at leading 
order the rate (4.32) can be simplified to 



mdydP, 24r„„,^(l-,)^ + ..., (4.36) 

with < ?/ < 1 and < P+ < Mb - 2Ei. 

In general, the hadronic tensor can be described in terms of the quantities 
P± = Eh =F \Ph\, whose true phase-space is M^r < P+ < P- < Mb, corresponding 
to a triangular region in the (P_,P+) plane. The variable P_ is related to our 
parton variables by P_ = n-p + A = rriby + A. In our theoretical description based 
on quark-hadron duality P+ starts from 0, while the small region with P_ < A is 
left unpopulated. This is illustrated in Fig. 4.2. Contours of constant hadronic or 
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leptonic invariant mass in the {P-,P^) plane are easy to visualize, since 

s^ = P| = P+P_, = {Mb - P+){Mb - P-) (4.37) 

are given by very simple expressions. The solid and dotted lines in the left-hand 
plot in Fig. 4.2 show the contours where sh = M^, and = {Mb — MdY, respec- 
tively, which can be used to separate B l^v events from semileptonic decays 
with charm hadrons in the final state. The dashed horizontal line shows the max- 
imum allowed value of P+ when a cut Ei > {M^ — M^) / {2Mb) is applied to the 
charged-lepton energy, which implies P+ < M'^/Mb- We will see later that this 
cut, which is obviously another way of eliminating the charm background, allows 
for a systematic treatment of the theoretical prediction, and is therefore of great 
interest. In the right-hand plot, we indicate the density of events in theory phase 
space obtained using the model of [37].^ It is apparent that the vast majority of 
events is located in the shape-function region of small P+ and large P_. 

In the remainder of this section, we present analytic results to the order of 
accuracy we have worked so far, for a variety of spectra in P — X„/~z/ decays. 
They are obtained by integrating over the scaling variable y before integrating 
over the hadronic variable P+, changing variables from P+ to = mi,{P^ — Co). 
The integral over the shape-function variable u is left until the end, so that our 
formulae are model independent. We will always present fractional decay rates 
normalized to the total inclusive rate 



T{B ^ Xul v) = Ttot = Ttree 



CFas{mk) /25 . 

i H ZTT 

Atx V 2 



+ . . . , (4.38) 



^While not rigorously implementing shape-function effects beyond tree level, the 
model of [37] has the advantage that it interpolates between the shape-function 
region and the remainder of phase space, where a local operator product expansion 
can be employed. On the contrary, our more rigorous discussion here is limited to 
the region of hard-coUinear jet momenta. However, the scatter plot shown in the 
figure provides a reasonably realistic impression of the population in phase space. 
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where the dots represent higher-order perturbative corrections as well as power 
corrections of order (AQCD/''^b)^ and higher. This procedure offers the advantage 
of eliminating the strong sensitivity to the heavy-quark (pole) mass. The integrals 
over the parton variable y encountered in our analysis can be reduced to a set of 
master integrals defined as 

b _ 



l + b 



h{b,z) = dyy 
Jo 

h{b,z) = j^dyyHny=^^[lnz-^^j, 

/3(M) = [dyyVy = ^ (in^ . - ^ + ^) , (4.39) 

Jo 1-y ^i + ^ + A i + b + jj' 

h{b,z)= \ dy y L2{1 - y) = -— L2{1 - z) — — — , 

Jo l + b l + b 

where b > —1 and z < 1 are a arbitrary real numbers. 
4.4.1 Charged-lepton energy spectrum 

Near the charged-lepton energy endpoint, where we can assume that Mb — 2Ei is 
of order a hadronic scale, the underlying event falls into the shape-function region 
of large n ■ p and small n ■ Ph- Starting from the triple differential rate in (4.36), 
we obtain for the normalized energy spectrum [36] 

Jo 1 47r 

2 mk{MB - 2Ei - u) 



Ttot dEi nib 
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+ (4/2(0) -3) In 



mbjMB - 2Ei - to) 



+ (7-7r2-3/2(a) + 2/3(a)) |. (4.40) 
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Here 



T(a) 
fn{a) 
Hia) 



6 Ji(l-a,l) -/i(2-a,l) 

4(1 - g, 1) - 4(2 - g, 1) 
/i(l-a,l)-4(2-a,l) ' 
llvr^ 49 



(4.41) 



+ (6 - c)h{a) - 4/3(a) - 2U{a) - 4/5(0) 



6 2 

At leading power in AgcD/'Tif, the heavy-quark mass in the denominator of the 
prefactor on the right-hand side of (4.40) can be replaced by m;, + uj = Mb — (^1 
which removes any sensitivity to the definition of m?,. This replacement can indeed 
be justified by studying power corrections to the shape function [59, 60]. 

All our results for decay rates will have a similar structure, but the definitions of 
the functions T, 4, and H will be different in each case. The tree-level result can be 
recovered by setting Vh = and a = 0, in which case T(0) = 1, and the spectrum 
is simply given in terms of an integral over the shape function [35] . Using the above 
result, it is straightforward to calculate the fraction Fe = T{Ei > Eo)/Ttot of all 
B — > Xul'v events with charged- lepton energy above a threshold Eq. Defining 
Ae = Mb- 2Eo, we find 



Fe{^e) 



Jo Mb-uj 



+ 



Atx 



2lu- '"^'^''r'" + (4/.(a)-7)lu 
+ (14 - ^2 _ ^ 2/3(0)) I . 



4. 

mb{AE - uj) 



(4.42) 



Note that the fraction Fe{Ae) is given in terms of a weighted integral over the 
shape function, with a weight factor of order Aqcd/''^6 that vanishes at the upper 
end of integration. As a result, only a small fraction of events is contained in the 
lepton endpoint region. 
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4.4.2 Hadronic P+ spectrum 

Applying a lower cut on the cliarged-lepton energy restricts the variable P+ to be 
less than A^;. However, a cut on P+ does not restrict the lepton energy to be in 
the endpoint region. Still, the fraction of events with P+ < A^; samples the same 
hadronic phase space as the lepton-endpoint cut, but it contains significantly more 
events. Such a cut therefore offers an excellent opportunity to determine the CKM 
matrix element \Vub\- 

To determine the fraction of events that survive, we integrate over x and y in 
the range < a; < y < 1 before integrating over P+. The result is 



Pp(A; 
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Air 
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mi)(Ap — uj) 



+ (7-7r2-3/2(a) + 2/3(a)) | 



(4.43) 



where now 



T(a) = 6/i(2 - a, 1) - 4/i(3 - a, 1) 
llvr^ 49 



H(a) 
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+ (6 - c)/2(a) - 4/3(a) - 2 /4(a) - Af^ia) - 4/5(0 



and 



/n(a) 
A/4(a) 



3/„(2 - g, 1) - 2J„(3 - g, 1) 
3/1(2 -g,l)- 2/1(3 -g,l) 
/4(3-g,l) 



(4.44) 



(4.45) 



3/1(2 - g, 1) -2/1(3 -g,l) ■ 
The contribution A/4 arises from the terms contained in the structure functions 

W4 and in (4.33). A comparison of the result for Pp(Ap) in (4.43) with the 

expression for Fe{Ae) in (4.42) yields that the cut on hadronic P+ contains a 

much larger fraction of all B Xul^u events. As a matter of fact, Fp{Ap) is 



102 



directly given in terms of an integral over the shape function, without a weight 
function of order Kc^cD/mb. (At tree level, Fp(Ap) = duj S{Cj).) Since the 
shape function peaks around ci) ~ A ^ 0.5 GeV, we expect a high efficiency for 
values of Ap in the vicinity of the "optimal cut" Ap = M^/Mb, which eliminates 
the charm background. 

4.4.3 Hadronic invariant mass spectrum 

A cut on the hadronic invariant mass in the final state constitutes the ideal separa- 
tor between B I'v and B ^ X^l'v events, since any final state containing a 
charm hadron has invariant mass above Mjj. Let us discuss the cut ^/sh < Md by 
examining the phase-space picture of Fig. 4.2. The available phase space for a cut 
on P+ < Ml) /Mb is fully contained. In addition there is a triangle-shaped region of 
larger P+, which culminates in a cusp where P+ = P_ = M^. Near the cusp, both 
light-cone momentum components are of the same order, and hence this portion of 
phase space should not be treated using our theoretical description based on the 
coUinear expansion. A priori, it is not evident that we can compute the fractional 
rate Fm{sq) = T{sh < so)/rtot in a controlled heavy-quark expansion. 

To see what happens, it is instructive to first ignore radiative corrections. At 
tree level, it is straightforward to obtain 

As „ 

aaV_ a; 



Fm{so) = jdCjS{Cj)+ JdLuSiiu)\^-^j (^2--^J, (4.46) 

A, 

where A^ = sq/Mb- The calculation of this event fraction requires knowledge of the 
shape function over a wider range in u than in the case of the event fraction with 
a cut on P+. The first integral is the same as for Fp{As) in (4.43) and corresponds 
to the region in phase space where P+ < Sq/Mb- The second integral corresponds 
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to the phase space above the dashed hne in Fig. 4.2. The region near the cusp 
corresponds to the upper integration region in the second integraL Note that, due 
to the rapid fall-off of the integrand, the tip of the triangle region only gives a 
power-suppressed contribution to the decay rate. When radiative corrections are 
included, the result for the integrated hadronic invariant mass spectrum becomes 
rather complicated. In general, we may split up the contributions into 



with F^r(so) = Fp(A, 



(4.47) 



The box contribution is given by the expression for the rate fraction Fp{Ap) in 
(4.43) evaluated with Ap = = sq/Mb- For the remaining contribution from 
the triangular region, we obtain [36]. 
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where 
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[7-7r2-3/2(a, z) + 2fsia,z)] 



contains the same functions T, H and /„ as defined in (4.44) and (4.45), but with 
all master integrals replaced by /„(6, 1) — * Inib, z). In addition, we need 
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where 



ki{z) = 4 6/i(2-a,2) -4Ji(3-a,2) = 4T(a, 2 



(4.51) 



6/2(2-0, z) -4/2(3-0, z) 



6/1(2- a, z) -4/1(3 -a, z) 



with the constant a given in (4.34). For orientation, a typical numerical value is 
a ~ 0.335 for /ij = 1.5 GeV. 

As mentioned above, the phase-space region near the cusp where u ~ ^/so or 
P ~ gives a power-suppressed contribution to the decay rate. More precisely, 
we find that the contribution is given by 



Fm(so) 3 e 



Cpasifii) 6 /A 



3— a 
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+ 



(4.52) 



TT (3 — a)^ \\/^y 3 — 

where the dots represent higher-order power corrections. For sq ~ mfeAqcD; the 
above result thus scales like (AQCD/'o^fi)^^""''^^, whereas /^a/(so) is of 0(1). For 
consistency, we should therefore omit the term in (4.52), which can be done by 
replacing all occurrences of y/so in upper integration limits in (4.48) and (4.50) 
with 00, which we will use this prescription in our numerical analysis in Section 4.6. 



4.4.4 Combined cuts on hadronic and leptonic invariant 
mass 

Bauer et al. have proposed to reduce the sensitivity to shape-function effects in the 
extraction of \ Vub\ by combining a cut on hadronic invariant mass with a cut > 
on the invariant mass squared of the lepton pair [61]. The first plot in Fig. 4.3 shows 
that this eliminates a large portion of the events with large P_. We indicate the 
remaining phase space for several choices of cuts, for example (sq, ^o) — i^D^ 0) 
(solid line), {M^, 6 GeV^) (dashed line), and ((1.7 GeV)^, 8 GeV^) (dotted line). 
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Figure 4.3: Phase-space constraints (left) and weight functions (right) 
for combined cuts on the hadronic and leptonic invariant 
mass: (sq, go) = i^D^ 0) (sohd), {Ml, 6 GeV^) (dashed), and 
((1.7GeV)2, SGeV^) (dotted). 



For the corresponding event fraction at tree level, we obtain 




A^/j/o 



(4.53) 

where = 1 — / {mi,M b) ■, and Ag = sq/Mb as above. For a fixed hadronic- 
mass cut Sq, the effect of the additional cut on is to broaden the support of the 
first integral, while at the same time reducing its weight due to the prefactor. To 
illustrate this point, we show in the second plot in Fig. 4.3 the weight functions 
under the integral with the shape function for the three different choices of (sq, go) 
mentioned above. The sensitivity to the precise form of the shape function is 
reduced because the weight functions become progressively more shallow as the 
value of go is raised. However, this reduction comes at the price of a significant 
reduction of the rate, raising questions about the validity of the assumption of 
quark-hadron duality. We will see in Section 4.6 that the relative uncertainty 
due to shape-function effects is not strongly (although somewhat) reduced when 
imposing an additional cut on g^. 
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4.5 Model-independent relations between spectra 

As we have seen in the previous section, all spectra and event fractions are ex- 
pressed in terms of weighted integrals of perturbative functions over the (universal, 
i.e. process independent) shape function. Therefore they all require the knowledge 
of the function form of the shape function, which cannot be calculated using ana- 
lytic techniques. One way would be to adopt a specific model or extract the shape 
function from experiment. Alternatively, it is possible to derive model- independent 
relations between different decay distributions in which the shape function has 
been eliminated [35]. The most promising strategy is to relate event fractions in 
semileptonic B l^v decays to a weighted integral over the B Xs'^ photon 

spectrum, which at present provides the most direct access to the shape function. 
While it is straightforward to derive such relations at tree level, radiative correc- 
tions introduce non-trivial complications [62, 63, 64, 65]. Since our formalism [36] 
has yet to be applied to the B Xg'j photon spectrum, we will instead derive a 
relation between the charged-lepton energy spectrum and a weighted integral over 
the P+ spectrum, which is in many respects very similar to the photon spectrum. 
We wish to construct a perturbative weight function w{A, P_|_) such that at leading 
power in AqcoMb 



dEi—= dP+w{A,P+)-—, A = Mb-2Eo. (4.54) 



This relation is independent of the shape function and hence insensitive to hadronic 
physics. The construction of the weight function is straightforward order by order 
in perturbation theory. Using the results of the previous section, we find 
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(4.55) 



where 



, , , 3952 - 5416a + 2988a2 - 838a3 + 120a^ - 7a5 
hi(a) = 2 — 2 



+ c 



(6-a)(4-a)2(3-a)(2-a)2 
20 - 8a + a^ 



(6-a)(4-a)(2-a) 



L / N , 20-8a + a2 

/12(a) = -4- -, (4.56) 

^ ^ 6-a 4-a 2-a ' ^ ' 



5056 - 6744a + 3556a2 - 942a3 + 127a^ - 7a 



(6-a)(4-a)2(3-a)(2-a)2 
The lesson we learned from this prototype relation (4.54) should be applied in 
the future to the B —>■ Xs'~f photon spectrum. Using similar methods, it will be 
possible to construct a shape-function independent relation of the form 

Fp(A) = i- l^dE, ^ w.{A, E,) , (4.57) 



where at tree level the weight function is simply Ws(A, E^) = 1. It is also possible 
(although far more complicated) to relate the B Xg'-f photon spectrum to the 
hadronic invariant mass spectrum Fm{so) in (4.47). However, because of the larger 
integration domain over the shape function, such a relation would require input 
of the photon spectrum beyond the region where it is currently experimentally 
accessible. 

The alert reader might wonder about the appearance of the renormalization 
scale fii in the weight function, since w{A,P^) is formally independent of the 
scale fii (because there is nothing to cancel a potential /i, dependence in (4.54)). 
Expanding the resummed result for the weight function to first order in ag, we 
obtain the simple expression 



I _ 2(A - P+) 

ll-loop Mr -P. 



l + ^(~-ln^-^^ 17 



47r \ 3 rrib 36 



(4.58) 
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Figure 4.4: Weight function w(A,P+) entering the rate relation (4.54) for 
A = M^/Mb and three different choices of the intermedi- 
ate scale, namely fii = 1.5 GeV (solid), 2.0 GeV (dashed), and 
1.0 GeV (dotted). The weight function is formally independent 

of /i,. 

in which the dependence on /ij has canceled. However, since this formula contains 
a large logarithm and the scale to be used in as is undetermined, it should not be 
used for phenomenological applications. A visualization of the resummed weight 
function in (4.55) for different choices of fii is given in Fig. 4.4. As is apparent, 
the function is pretty uneventful, since apart from the rational prefactor, the P+ 
dependence is given through a single logarithm. 

4.6 Numerical results 

We are now ready to study the implications of our analysis for phenomenology. 
We start by deriving the numerical values for the shape-function mass and kinetic 
energy including errors. We then present a model for the shape function which 
satisfies all theoretical constraints, as derived in Chapter 3. Finally, we present 
numerical results for the various decay rates and spectra investigated in Section 4.4. 
We use the two-loop running coupling constant in the MS scheme, normalized 
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such that as{Mz) = 0.119. Our standard choice of the intermediate matching 
scale is fii = 1.5 GeV. This corresponds to setting fif = mfeAhad with a typical 
hadronic scale Ahad ~ 0.5 GeV. The values of the strong coupling evaluated at 
these scales are as{mh) ~ 0.222 and as{^ii) — 0.375. The corresponding values of 
the perturbative parameters a and c defined in (4.34) are a ^ 0.335 and c ~ 0.614. 
Finally, the leading-order Sudakov factor in (4.33) takes the value e^"^"^'"^^'^ ^ 1.21. 

4.6.1 Shape-function mass and kinetic energy 

A value for the shape-function mass can be obtained by combining the relations 
(3.58) or (3.59) with existing predictions for the 6-quark mass in the relevant renor- 
malization schemes. The potential-subtracted mass at the scale /i/ = 2 GeV has 
been determined from moments of the bb cross section and the mass of the T(15') 
state [66]. Using the first relation in (3.59), we find [36] 771^^(2 GeV, 2 GeV) = 
m^^(2 GeV) = (4.59±0.08) GeV. From a similar analysis the kinetic mass has been 
determined at the scale /i/ = 2 GeV to be mf''{l GeV) = (4.57 ± 0.06) GeV [49]. 
From the second relation in (3.59) it then follows that m^^(l GeV, 1 GeV) = (4.65± 
0.06) GeV. Computing the scale dependence of the shape-function mass (using the 
fact that the pole mass is RG invariant), we obtain at the intermediate scale the 
values m^^(/ii,/ii) = (4.61 ±0.08) GeV and m^^(/ii, /Xi) = (4.65±0.06) GeV, respec- 
tively. Alternatively, we may use relation (3.58) in conjunction with an experimen- 
tal determination of the 6-quark pole mass from moments of inclusive B ^ X^l^v 
and B X^'j decay spectra. Using the average value Apoie = (0.375 ± 0.065) GeV 
obtained from [67, 68, 69, 70], we find mf^{fii,fii) = (4.67 ± 0.07) GeV. It is quite 
remarkable that these different determinations of the shape-function mass, which 
use rather different physics input, give highly consistent results. Combining them. 
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we quote our default value for the shape-function mass at the intermediate scale 
fii = 1.5 GeV as 

mfifii, fii) = (4.65 ± 0.07) GeV . (4.59) 

The corresponding A parameter is A(/ij,/ij) = (0.63 ± 0.07) GeV. 

A value of the kinetic-energy parameter in the shape-function scheme can be 
obtained from (3.60) or (3.61). Using the first relation and the experimental value 
-Ai = (0.25 ± 0.06) GeV^ [68, 69, 70] yields /i^(/ii,/ii) = (0.271 ± 0.064) GeV^. 
Alternatively, we may use the result for the kinetic-energy parameter obtained in 
the kinetic scheme, GeV)]kin = (0.45 ±0.10) GeV^ [49], to get from (3.61) the 
value fiKfii, fii) = (0.254±0.107) GeV^. Again, the two determinations are in very 
good agreement with each other. Combining them, we obtain 

nKlii, /ii) = (0.27 ± 0.07) GeV^ . (4.60) 
4.6.2 Model shape functions 

In our analysis of decay rates below, we will adopt a model for the shape function 
S{i2},fii) at the intermediate scale. For the purpose of illustration, we use a two- 
component ansatz for the shape function that is a generalization of the model 
employed in [37, 41]. The form we propose is [36] 

^ ^ (|)'-:.P (-4) , ,) , 

(4.61) 



where A and b are model parameters, and A differs from the pole-scheme parameter 
^poic by an amount of 0(as(At)). In the limit ^^(/i) — > this function reduces to 
the familiar model used in [37, 41]. The radiative tail ensures the correct leading 
asymptotic behavior of the shape function as displayed in (3.65). This in turn gives 
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the correct power-like dependence of shape-function moments on the integration 
cutoff. In our model, this tail is glued onto a "primordial", exponential function 
such that the combined result is continuous. The normalization factor is given 
by 
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(4.62) 



TT 

which is determined such that the integral over the shape function from ct; = 
to fif + A{fif,fi) coincides with the first expression in (3.63) up to second-order 
corrections. By evaluating the first moment of the model shape function, we find 
that the model parameter A is related to the HQET parameter A in the pole 
scheme and the shape-function scheme as 

A = Apoic H 1= , A = A(/ii, ^i) + ^i[—= - 1 ] . (4.63) 

Finally, the model parameter b can be adjusted to reproduce a given value for the 
second moment of the shape function. 

Table 4.1 collects the parameters of the model shape functions at the intermedi- 
ate scale fii = 1.5 GeV corresponding to different values of A(/ij, /ij) and /i^(/ii, /ij). 
The left-hand (right-hand) plot in Fig. 4.5 shows three models for the shape func- 
tion obtained by varying the parameters A and /i^ in a correlated (anti-correlated) 
way. In both cases, the solid, dashed, and dotted curves refer to different values 
of A, as indicated in the table. In Fig. 4.6 we illustrate the renormalization-group 
evolution of the shape function as studied in Chapter 3. The sharply peaked solid 
line shows our model function evaluated with A = 0.495 GeV and b = 3.0, which 
we use as an ansatz for the function S{uj,^q) at the low scale fio = 1 GeV. (For 
comparison, the dotted gray curve shows the default choice for the shape function 
adopted in [37, 41], which exhibits a very similar shape except for the missing 
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Table 4.1: Parameters and moments of the model shape functions at the 
intermediate scale /ij. The running quantities mf^, A, and /i^ are 
defined in the shape-function scheme and evaluated at Hf = = 
fii = 1.5 GeV. 



Model 


Lines 


m^* [GeV] 


A [GeV] 


fil [GeV2] 


A [GeV] 


b 


SI 


Dotted 


4.72 


0.56 


0.20 


0.611 


2.84 


S2 








0.27 


0.617 


2.32 


S3 








0.34 


0.626 


1.92 


S4 


Solid 


4.65 


0.63 


0.20 


0.680 


3.57 


S5 








0.27 


0.685 


2.93 


S6 








0.34 


0.692 


2.45 


S7 


Dashed 


4.58 


0.70 


0.20 


0.751 


4.40 


S8 








0.27 


0.753 


3.61 


S9 








0.34 


0.759 


3.03 



radiative tail.) The broad solid curve gives the shape function at the intermediate 
scale yUj = 1.5 GeV as obtained from the evolution equation (3.52). The barely vis- 
ible dashed-dotted curve shows our default model for the function S{u},^i), which 
coincides with the solid line in the left-hand plot. The beautiful agreement of the 
two curves gives us confidence in the consistency of our models adopted for the 
shape function at the intermediate scale. 

4.6.3 Predictions for decay spectra and event fractions 

We are now ready to present our results for the decay spectra and partially inte- 
grated event fractions in B ^ X^l^v decays. In order to illustrate the sensitivity 
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Figure 4.5: Various models for the shape function at the intermediate scale 
/ij = 1.5 GeV, corresponding to different parameter settings in 
Table 4.1. Left: Functions SI, S5, S9 with "correlated" parame- 
ter variations. Right: Functions S3, S5, S7with "anti-correlated" 
parameter variations. 
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Figure 4.6: Renormalization-group evolution of a model shape function from 
a low scale /iq (sharply peaked solid curve) to the intermediate 
scale /ij (broad solid curve). See the text for an explanation of 
the other curves. 



to shape-function effects we use all nine shape functions SI through S9 in Table 4.1, 
thus varying the parameters A(yUj,/ij) and /i^(/ij,/ij) independently. This is quite 
conservative because we neglect any possible correlation between them. For each 
physical quantity we draw three bands corresponding to the three different values 
of A. The width of each band reflects the sensitivity to the variation of /i^. 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Ap [GeV] A, = So/Mb [GeV] 



Figure 4.7: Fraction of 5 — > X„/^P events with hadronic light-cone mo- 
mentum P+ < Ap (left), and fraction of events with hadronic 
invariant mass sh < sq (right). In each plot, the three bands 
correspond to the values A = 0.63 GeV (solid curves), 0.70 GeV 
(dashed curves), and 0.56 GeV (dotted curves). Their width re- 
flects the sensitivity to the value of /i^ varied in the range between 
0.20 and 0.34 GeV^. The arrow indicates the point at which the 
charm background starts. 

The following predictions for spectra and rate fractions refer to the leading term 
in the heavy-quark expansion. We note that our calculations would break down 
if the cuts on kinematic variables were taken to be too strict, because then the 
spectra would become dominated by hadronic resonance effects. Parametrically, 
this happens when the quantities Ap, A^, or A^ become of order Aq^p/Mp ~ 
50MeV. On the other hand, taking too large values of Ap, A^, or Ag such that 
they are not of O(Aqcd) anymore invalidates the collinear expansion. 

In Fig. 4.7 we show predictions for the fractions of all B X^l^v events 
with hadronic light-cone momentum P+ < Ap, and with hadronic invariant mass 
squared sh < sq. Recall that, for Ap = A^ = sq/Mb, the hadronic invariant 
mass fraction Fm differs from the fraction Fp by the contribution of the events 
in the triangular region above the dashed line in Fig. 4.2. Comparing the two 
plots, we observe that this additional contribution is predicted to be very small. 
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2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 0.2 0.4 0.6 0.8 



[GeV] Ae = Mb-2Eo [GeV] 

Figure 4.8: Charged-lepton energy spectrum in the region near the kinematic 
endpoint (left), and fraction of events with charged-lepton energy 
El > Eq (right). The meaning of the bands and the arrow is the 
same as in Fig. 4.7. 

(Note that for large values of we even find a negative contribution to the rate 
from the triangle region for some choices of the shape function. This feature is 
unphysical and should be fixed by the inclusion of power corrections to our leading- 
order predictions.) The arrows on the horizontal axes indicate the points Ap^g = 
Ml) /Mb, beyond which final states containing charm hadrons are kinematically 
allowed. With this choice of the cut, both rate fractions capture about 80% of 
all events. While it is well known that a hadronic invariant mass cut ^/sH < 
provides a very efficient discrimination against charm background [71, 72, 73], here 
we observe that the same is true for a cut on the P+ variable. 

Our results for the lepton energy spectrum Si{Ei) = (l/Ftot) (dT/dEi), and for 
the event fraction with a cut Ei > Eq, are displayed in Fig. 4.8. The right-hand plot 
shows that with A^; = M^^/Mb only about 10-15% of all events are retained, and 
the theoretical calculation is very sensitive to shape-function effects. Such a cut is 
therefore much less efficient than the cuts on sh or P+. As a result, an extraction 
of \Vuh\ from the charged-lepton endpoint region is theoretically disfavored. 



116 

The shape- function sensitivity is rather small for values of Ap and near the 
charm threshold. This is to some extent a consequence of our improved knowledge 
of the shape-function parameters. On the other hand we observe an interesting "fo- 
cus mechanism" in that the three bands in the P+ and the hadronic invariant mass 
(but not for the charged lepton energy) event fractions start to converge near the 
charm threshold. This is due to a subtle interplay between the jet function and the 
shape function, which is explained in more detail in [36]. The important observa- 
tion is that the leading logarithm in the jet function has the opposite sign as found 
in a straightforward partonic calculation. The convergence is somewhat surprising 
since the model shape functions at the intermediate scale /ij are rather different 
for values around u ~ Aqcd (see Fig. 4.5), and the argument that they share the 
same norm is not applicable yet. One way of thinking about this mechanism is 
to notice that the broadening of the shape function under renormalization-group 
evolution from a low scale up to the intermediate scale (Fig. 4.6) is a perturbative 
effect, which should not lead to an increased shape- function sensitivity. Because 
the convolution of the shape function with the jet and hard functions is indepen- 
dent of the scale /ij, the broadening of the shape function must be compensated 
by perturbative logarithms in the jet function. 

This focus effect did not take place in earlier studies such as [37, 71, 72, 73, 74], 
where parton-model spectra were convoluted with a primordial shape function. As 
mentioned earlier, in the parton model the leading Sudakov logarithm comes with 
the opposite (negative) sign, hence causing an anti-focus effect of the radiative 
corrections. This also explains why our prediction for the hadronic invariant mass 
fraction Fm exhibits a smaller shape-function sensitivity than what has been found 
in most previous analyses. 
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w [GeV] Ap [GeV] 

Figure 4.9: Left: Four examples of shape functions with identical normal- 
ization and first two moments, corresponding to A(/ij,/ij) = 
0.63 GeV and = 0.27 GeV^, but different functional 

form. Right: Corresponding results for the event fraction 
Fp{Ap). 

We shall stress again that the knowledge of the first few moments does not 
determine the functional form of the shape function. In order to demonstrate 
this, we have constructed several functions of different functional forms, but with 
identical norm, first, and second moment. To correctly account for the asymptotic 
behaviour u ^ Aqcd, we have added the second term in (4.61) to each of the 
functions. They are shown on the left-hand side of Fig. 4.9. All of them are thus 
allowed by the required moment analysis and could be considered as variations of 
the shape function. Some of the curves are admittedly extreme, given the fact 
that the b photon spectrum is identical to the shape function at tree level. 

Any information about the shape function can be used to eliminate some of these 
curves. At present, however, we stress that using the model shape functions SI 
through S9 sample even the variation of the functional form, as the event fraction 
predictions on the right-hand side of Fig. 4.9 fall within the bands studied earlier. 

We summarize our phenomenological results in Table 4.2, in which we compare 
the shape-function sensitivity (which is the only uncertainty considered here) of 
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Table 4.2: Comparison of different theoretical methods using inclusive B- 
decay rates to extract the CKM matrix element \Vub\- The error 
on the efficiency represents the sensitivity to the shape function 
only. All results refer to the leading term in the heavy-quark 
expansion. 



Method 


Cut 


Efficiency 


Hadronic invariant mass 


SH < Ml 


(81.4^|-2)% 




sh < (1.7 GeV)^ 


(78 2+H)% 

V ' "'^ —5.27 ' ^ 




sh < (1.55GeV)2 


(72.7i^i)% 


Hadronic P+ 


P <Ml = 0.66 GeV 

imb 






P+ < 0.55 GeV 




Charged-lepton energy 


El > = 2.31 GeV 


(12.5^:1)% 




El > 2.2 GeV 


{22.2tli)% 


Combined {sh, (f) cuts 


SH<Ml, q^>0 


{7A.6ttl)% 


[tree level only] 


Sh < Ml, q^>6 GeV^ 






Sh < (1.7GeV)2, > 8 GeV^ 


(33.4^^1)% 



different kinematic cuts. The first two blocks of entry address the high-efficiency 
methods of cutting on the hadronic invariant mass sh and the variable P+. The 
efficiency for the optimal cut sh Ml is remarkably precise; however, due to 
detector resolution effects it is necessary to relax that cut somewhat. We state 
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results for upper limits on the invariant mass of 1.7 GeV and 1.55 GeV, for which 
the relative uncertainty increases somewhat. Cutting on P+ introduces a relative 
uncertainty of order 10%, which is still acceptable. Again, experimental constraints 
might force us to move away from the optimal cut. However, it might be possible 
to stay closer to the optimal cut due to the apparent "buffer zone". (We will 
discuss the charm background in much more detail below.) 

The only low-efficiency method we discuss is a cut on the charged-lepton en- 
ergy. Unlike the previous two methods, such a measurement does not require the 
reconstruction of the neutrino in the decay B Xu and is therefore favored by 
experiment. However, due to the low efficiency one has to worry also about other 
theoretical uncertainties like weak annihilation effects [75], which are expected to 
contribute less than 3% to the total rate and can therefore be safely neglected in 
high-efficiency methods like the sh and P+ cuts. Our prediction for the lepton 
energy event fractions are significantly larger than have been reported in the past, 
and suggest that the extracted value of \Vub\ needs to be corrected to lower values. 

For completeness, we also give results for some combined cuts on hadronic 
and leptonic invariant mass, which have been briefly discussed in Section 4.4.4. 
Contrary to the other cases, these numbers refer to the tree-level approximation 
and so should be taken with caution. For reference, we quote again the (tree-level) 
result for the pure hadronic invariant mass cut, which differs significantly from 
the corresponding result including radiative corrections. While the additional cut 
on leptonic reduces the shape- function sensitivity, it comes along with a strong 
reduction of the efficiency. For instance, the combined cut ^/sH < 1.7 GeV and 

> 8 GeV^ employed in a recent analysis of the Belle Collaboration [76] has an 
efficiency of about 33% (at tree level and leading order in Aqcd/"^6), which is much 
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smaller than the efficiency of the pure hadronic invariant mass cut ^/sh < 1-7 GeV. 
However, the sensitivity to shape-function effects is only slightly better in the case 
of the combined cut. 



4.7 Charm background 

One of the main advantages of the P+ spectrum over the hadronic-mass spectrum 
is a better control of the charm background. In order to study this in more detail, 
we investigate the OPE prediction for the normalized = n ■ p/rrii, spectrum in 
B ^ Xcl^u decays, 

^ ^ = ^^|r^ k(3 - 2p+) + qP+{?>- 8p+ + ?>Pl) - q\2 - 3^+)! , (4.64) 

where Q < p+ < ^/g with g = (mc/mb)^, and 

/(^) = 1 - 8^ + 8^^ - - 12^2 In ^. (4.65) 

We only include tree-level contributions from dimension-3 operators [77]. 

In this approximation, inclusive charm events are located along a single line 
p+P- = m^, or in terms of the hadronic variables (P+ — A)(P„ — A) = {M^ — A)^, 
in Fig. 4.2. This line starts at the tip of the triangle = = Md and extends 
to the right while always staying above the solid line sh = M"^. Because 0{as) 
corrections redistribute these events into the light-gray segment above that line, 
the integral over the spectrum in (4.64) serves as an upper bound on the inclusive 
charm background. In the same approximation, the hadronic mass distribution is 
given by 

1^ ^ = TT^V^^ \z{3-2z)-g{A-3z)] , with z = ~^ . (4.66) 

Tc ds f{Q)e L J e 
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0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 

Ap [GeV] A, [GeV] 

Figure 4.10: Fraction of i? ^ Xul~v events with P+ < Ap (left) and sh < 



So (right). The bands are identical to the ones in Fig. 4.7. The 
curves represent the background from inclusive B — > X^l'i^ 
(solid) and exclusive B — > D^*H~iy decays (dashed), normalized 
to the total inclusive semileptonic charm rate. 

Here, s = sn/ml, e = A/mb, and phase space is such that {^/Q + s)^ < s < 
{g + e){l + e). 

The results for the event fractions Fp and Fm are summarized in Fig. 4.10. The 
bands show the inclusive B —>■ X^l'v event fractions for different shape-function 
models as discussed above, while the solid lines give predictions for the inclusive 
B Xc/~z/ background spectra as obtained from equations. (4.64) and (4.66). 
The latter are normalized to their total decay rate F^, which is about 60 times 
larger than F^. In the case of the fraction Fm-, the inclusive charm background 
starts right at the threshold = [nic + KY /Mb — M'^/Mb- In the case of the 
Fp event fraction, on the contrary, the inclusive charm background starts at a 
value Ap = m'^/mi) + A ^ 0.96 GeV, which is significantly larger than the value 
Ap = M^i/Mb ~ 0.66 GeV above which final states containing charm mesons 
are kinematically allowed. This is a consequence of the fact that there is a gap 
between the charm threshold sh = M"^ and the start of the inclusive h ^ c events. 
In reality this gap is filled by exclusive modes containing the hadronic final states 
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D, D*, or Dn, Dnn. The exclusive contributions to the hadronic P+ spectrum 
from B D^*H^iJ decays are 



dP+ 487r3 

{w^ - 1) {w + iy 



dP+ 487r3 P+ 

X l + ^^^i (4.67) 



w = v-v' 



w + 1 

where M^(,) /Mb < P+ < M^(.) andr(^,) = Mj~,{*)/Mb- The recoil variable 
is given by M^j,, + = 2w M^(,)P^. The corresponding contributions to the 
hadronic mass spectrum are given by 

^^=r^(,)-5(sH-M^„). (4.68) 
asH 

The exclusive contributions to the event fractions are given by the dashed lines 
in Fig. 4.10. To obtain these curves we have used an ansatz for the form fac- 
tors (w) that is consistent with experimental data on the recoil spectra and 
branching fractions. The fact that these exclusive modes are very well understood 
should help to model the background. The smooth onset of the D^*^ background 
is a direct consequence of the fact that the ideal P+ cut touches the charm region 
at only a single point in phase space (see Fig. 4.2). On the contrary, the region 
of phase space when applying the sh cut borders the background along the curve 
separating the light- and dark-shaded regions, which leads to a step increment in 
the event fraction Fm- As a consequence, one needs to move away from the ideal 
cut Sh = Ml) because of smearing effects due to experimental resolution. It is 
our hope that it will be possible to stay closer to the ideal cut when performing 
a P+ analysis. In this case, both the P+ and sh discrimination methods lead to 
comparable efficiencies and shape- function uncertainties, as emphasized earlier in 
the discussion after Table 4.2. 



123 



4.8 Theoretical accuracy of a \Vub\ measurement 

As a final remark to this long Chapter, let us comment on the applicability of 
the theoretical framework developed in this work. According to Fig. 4.2 most 
of the B — > Xul'v events are located in the shape-function region of large P_ 
and small to moderate P+, in which a systematic heavy-quark expansion using 
SCET power counting is valid. It allows us to calculate inclusive decay rates 
integrated over domains AP_ ~ Mb and AP+ ^ Mb-, where typically AP+ ~ 
^QCD- (In the examples above, AP+ = A^;, Ap, or As, respectively.) While the 
corresponding predictions for decay spectra and event fractions are sufficient to 
analyze experimental data over most of the phase space relevant to measurements 
of the CKM matrix element \Vub\i it would be of interest to extend the validity 
of the theoretical description outside the shape-function region. Let us in the 
following paragraph comment on such a prospect. 

In the case where AP„ ~ AP+ ~ Mb are both large, the decay spectra can be 
computed using a local operator product expansion. The resulting prediction for 
the normalized B — > X^l'v spectrum in the variable = p^/nib = (P+ — A)/mb 
reads [77] 



1 rfr„ 





463 






as 


H 


3^ 




/ 17A, ^ 


\18 


mi 



ln]5+ 

P+ 

3X2 

2ml 



Ai 

6ml 



(4.69) 



where < p+ < 1, and 



158 407)0 367p^ llSn^ lOOn^ llp^ 

Hp) = 1 1 1 

9 18 6 3 9 6 1 



7p6 



4 A6p , 16p^ 
gp- _| 

3 3^3 



\np — Ap^(3 — 2p) p . 



(4.70) 
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We include the contributions from dimension-3 operators at 0{as) and those from 
dimension-5 operators (whose matrix elements are proportional to the heavy-quark 
effective theory parameters Ai,2) at tree level, using [55, 37]. If the p+ spec- 
trum is integrated without a weight function, the tree-level power corrections from 
dimension-5 operators vanish. This is in accordance with the fact that subleading 
shape functions have zero norm at tree level [59]. 

An interesting question is whether it will be possible to match the Factorization 
approach for the shape-function region AP+ ~ Aqcd and the OPE approach for 
AP4. ~ Mb in some intermediate region of AP+ values that are numerically (but 
not parametrically) large compared with Aqcd- If the two predictions were to agree 
in an overlap region, this could be used to construct a theoretical description of 
inclusive B X„ l^u decay distributions that is valid over the entire phase space. 
While this is an exciting prospect, we note that performing a systematic operator 
product expansion in the overlap region is far from trivial. For a hierarchy of scales 
Aqcd ^ Ap ^ we obtain with A = Ap — A 

V A 3 A 36 J nib \3 A 9 J \ 

(4.71) 

For Ap ~ 1 GeV, the power correction in the second term leads to an enhancement 
of the fraction Fp by 5-10%. This is of similar magnitude as tree-level estimates 
of (zero-norm) subleading shape-function effects on the Ei and sh spectra [60, 78]. 

The leading-order term in Eq. (4.71) can also be obtained from Eq. (4.43) by 
taking the limit Ap ^ Aqcd- Interestingly, such an analysis uncovers that in 
this kinematic range there is an enhanced class of power corrections of the form 
(Aqcd/^)" with n > 2, which arise first at 0{as) (see also [79]). The leading 
corrections to the expression above are given by 



Fp(Ap) = 1 - 
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^Fp(Ap) = -fi(21u^-n4- (^-^2) 

37r V A Sy 3a' 
It follows that for sufficiently large A the event fraction can be expressed as a 
double expansion in A/mf, and Aqcd/A. Because of the hierarchy of scales AP^ -C 
AP+ Mb <C M^, again a two-step procedure is in order. (An example of such an 
approach can be found in [80].) Such a multi-scale OPE is non-trivial and left 
for future work. For now we merely note the estimate that subleading power 
contributions to the event fraction Fp result to about 10%. 

We are now ready to give an in-depth look at the prospect of determining \ Vub\ 
from a P+ measurement. A list of relevant theoretical uncertainties to the fraction 
of events with the optimal cut < M'^/Mb is [77] 

Fp = (79.6 ±10.8 ±6.2 ±8.0)%, (4.73) 

where the errors represent the sensitivity to the shape function (as stated in Ta- 
ble 4.2), an estimate of 0{al) contributions, and power corrections, respectively. 
The uncertainty due to our ignorance of the shape function is obtained in part by 
varying the parameters A = Mb — rrif, (4.59) and fi"^ (4.60) in the shape-function 
scheme, determining its first two moments {SFp = lg|%). Furthermore we vary 
the functional form of S, as depicted in Fig. 4.9 {SFp = 

Higher-order perturbative effects are estimated by studying the dependence on 
the matching scales /ij and fih, varied in the ranges 1.25 GeV < /.ij < 1.75 GeV 
{5Fp = ttl%) and m^/v^ <fih<V2mh {5Fp = tti%)- 

Every single error stated here can be reduced in the future within the sys- 
tematic framework presented in [36, 77]. A proper treatment of subleading power 
corrections in the shape-function region is now feasible. In addition, the leading 
shape-function uncertainty can be eliminated using model-independent relations 
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such as (4.57). Finally the perturbative uncertainty can be reduced by computing 
the 0{a'j.) corrections to equation (4.43). 

The CKM-matrix element | Vub \ can be extracted by comparing a measurement 
of the partial rate r„(P+ < Ap) with a theoretical prediction for the product of 
the event fraction Fp and the total inclusive B X„/"z/ rate. The resulting 
theoretical uncertainty on \Vub\ up to date is 

= (±7±4±5±4)%, (4.74) 

I Kb I 

where the last error comes from the uncertainty in the total rate (4.38) as can be 
found in [81, 82]. Because of the large efficiency of the P+ cut, weak annihilation 
effects [75] have an influence on {V^bl of less than 2% and can be safely neglected. 

QCD-Factorization has provided us with a powerful method to systematically 
improve the calculation of event distributions in the shape-function region. The 
hadronic P+ spectrum is the preferred application of this framework, and the 
possibility of a high-precision determination of \Vub\ through a measurement of the 
partial decay rate r„(P+ < Ap) is an exciting prospect. 



CHAPTER 5 
EXCLUSIVE RADIATIVE DECAYS 

The radiative, semileptonic decay B ^ '-/ lu provides a clean environment for 
the study of soft-colhnear interactions [83]. This process is particularly simple in 
that no hadrons appear in the final state. Yet, there is sensitivity to the light- 
cone structure of the B meson, probed by the coupling of the high-energy photon 
to the soft spectator quark inside the heavy meson. Therefore this decay mode 
serves as an opportunity to learn QCD factorization and Sudakov resummation for 
exclusive decays in a realistic environment [24]. Our main goal is to establish the 
QCD factorization formula [84] 

poo (h^il u] 

A{B- ^^l-iyi)(xMBfBQu T{l+,E^,m,,ii) (5.1) 

to all orders in perturbation theory and at leading power in AqcD/^ib. Here Qu = \ 
is the electric charge of the up-quark (in units of e), and T = 1 + 0{as) is a per- 
turbative hard-scattering kernel. The decay constant fs and the LCDA cf)^ have 
been introduced in Section 3.1. The physics underlying the factorization formula 
is that a high-energy photon coupling to the soft constituents of the 5-meson 
produces quantum fluctuations far off their mass shell, which can be integrated 
out in a low-energy effective theory. As usual we choose the photon direction to 
be the ^-direction, so that the photon momentum is E^n. The two transverse 
polarization states of the photon can be expressed in terms of the basis vectors 
£^ = -^(0, 1, =F^;0), which correspond to left- and right-circular polarization, re- 
spectively. 

The new complication, as compared to the factorization of the inclusive case of 
the last Chapter, is that the final theory is SCET//. Radiative decays are partic- 
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ularly simple in that the only external "collinear" particle is the electromagnetic 
photon, which carries a purely light-like momentum. Therefore the only hadronic 
constituents present are the HQET fields describing quarks and gluons of the B 
meson. While it is possible to match QCD directly onto this theory, a two-step 
matching procedure through the intermediate theory SCET/ has the advantage of 
integrating out off-shell propagators of order and E'^Aqcd in separate steps. 
This suggests a second stage of "perturbative" factorization [14, 85, 86], which we 
will establish below. It says that the hard-scattering kernel itself can be factorized 
as 

T(/+, E,, m„ f,) = H f ^, ■ J f ^) . (5.2) 

Factorization holds as long as the photon is energetic in the 5-meson rest 
frame, meaning that E^ is of the order of the 6-quark mass. In order to prove the 
factorization formula (5.1) one needs to show that [4]: 

1. The decay amplitude can be expanded in powers of transverse momenta (i.e. 
in powers of Aqcj:,/E^) and, at leading order, can be expressed in terms of a 
convolution with the 5-meson LCDA as shown in (5.1). 

2. After subtraction of infra-red contributions corresponding to the i?-meson 
decay constant and LCDA, the leading contributions to the amplitude come 
from hard internal lines, i.e., the hard-scattering kernel T is free of infra-red 
singularities to all orders in perturbation theory. 

3. The convolution integral of the hard-scattering kernel with the LCDA is 
convergent. 

4. Non-valence Fock states do not give rise to leading contributions. 
In the next Section we will address all of the above points. 
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5.1 Proof of factorization 

In the full theory, the hadronic part of the decay amplitude for B 7/// is 
given by the time-ordered product of a weak, flavor-changing current (u6)v-a and 
the electro-magnetic currents °™^b). However, attachments of the 

energetic photon field to the massive 6-quark line leads to power suppression [14], 
and so at leading power in SCET, the time-ordered product is matched onto trilocal 
operators of the form 

5^ ze /" d'xT { [uY{l - l,)b] (0), [g 4^^-) q] (x)} 

q=u,b 

Y [ dsdtCi{t,s,v-q,mb,fi)Qs{tn)4''^^\sn)^ri:K{0). (5.3) 

The Wilson coefficients receive contributions only from Feynman diagrams 
where the photon is emitted from the spectator quark line. The fields 4^'^™) and 
are located at the appropriate light cones, which refiects the SCET scaling prop- 
erties of the photon and spectator-quark momenta. ^I.'j™^ is the electromagnetic 
analog of the gauge-invariant collinear gluon field defined in (2.60). To first order 
in e we have 

/•O 

4(7^(0) = n„7M dweF^^'iwn). (5.4) 
J —00 

The Feynman rule for this object is simply e^*, where e is the photon polarization 
vector. Finally, from the fact that the leptonic weak current z^7^(l — 75) / is con- 
served (in the limit where the lepton mass is neglected) it follows that the relevant 
Dirac structures Fj in (5.3) can be taken as Fi = 7^(1 — 75) and F2 = n^(l + 75). 

We can use the trace formalism in HQET and the definition of the 5-meson 
LCDA (3.1). Performing the relevant traces, we find that at leading power in 
^QCB /frib the decay amplitude vanishes if the photon has right-circular polariza- 
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tion, while for a photon with left-circular polarization it is given by 

A{B^ ^il^i^i) = -^Vubeui{pi)f_{l - -f5)vu{pu) (5.5) 
X >/Mb J duj Ci(cj, n-q,v -q, nib, /i) (p+iuj, fi) + . . . , 

where the dots represent power-suppressed contributions. Only the SCET oper- 
ator with Dirac structure Fi = 7''(1 — 75) contributes to the decay amplitude. 
Ci denotes the Fourier transform of Ci, and we will continue the discussion in 
momentum space. The HQET parameter F{fi) is related to the physical 5-meson 
decay constant through fsV^B = Kp{miy, jj) F{jji) [1 + 0(AQCD/''Tib)], where at 
next-to- leading order in the MS scheme [87] 

i^.K,/^) = l + ^^f31n!^-2V (5.6) 
47r V /i y 

Combining these results, it follows that the terms shown in the second line of (5.5) 
equal those on the right-hand side of the factorization formula (5.1) if we identify 
the hard-scattering kernel as 

^T(/+,E^,mb,/i) = K~^{mb,fi) Ci(/+, 2£'^, E^, m^, /z) . (5.7) 

We thus completed the first of the four steps to prove factorization. 

We proceed to prove the convergence of the convolution integral in (5.1). The 
key ingredient here is to note that the invariance of SCET operators under repa- 
rameterization of the light-cone basis vectors n and n can be used to deduce the 
dependence of Wilson coefficient functions on the separation t between the com- 
ponent fields of non-local operators. This has been discussed in some detail in 
Section 2.5. In our case, invariance of the operators in (5.3) under the type III 
rescaling transformation n'^ n'^/a and n'^ — > an^ (with fixed v) implies that 

Ci(t, n ■ q,v ■ q, mi,, /i) = Ci{at, an ■ q,v ■ q, nib, (5.8) 
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to all orders in perturbation theory. In other words, the variables t and n-q can only 
appear in the combination n ■ q/t, but not individually. The coefficient functions 
can be factorized further by noting that v ■ q and rrib enter only through hard 
or hard-coUinear interactions with the heavy quark. The corresponding modes 
can be integrated out in a ffist matching step onto SCET/ and lead to a functions 
Hi{v-q, m^, /i), which depend on the Dirac structure of the weak current containing 
the heavy quark. The non-localities of the component fields in (5.3) result from 
the coupling of hard-coUinear fields to the soft spectator quark in the B meson. 
These effects live on scales of order mfeAqco and can be integrated out in a second 
step by matching onto SCET//, leading to the function J(t,n ■ q,fi). Therefore 



where = 2v -q/mh = 2E^/miy is a scaling variable of order 1. The corresponding 
result for the hard-scattering kernel obtained after Fourier transformation has the 
form shown in (5.2) if we identify 



where n ■ ql+ = 2E^l^. Since the dependence of the coefficient functions on the 
renormalization scale is logarithmic, it follows that to all orders in perturbation 
theory the Wilson coefficients in (5.9) scale like Ci ~ 1 modulo logarithms. (Cor- 
respondingly, the kernel scales like T ~ 1 modulo logarithms.) The convergence of 
the convolution integral in (5.5) in the infra-red region t oo, corresponding to 
the region ^ in the factorization formula (5.1), then follows to all orders in 
perturbation theory as long as the integral converges at tree level. Because the B 
meson has a spatial size of order 1/Aqcd due to confinement, the bilocal matrix 




(5.9) 




(5.10) 
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element must vanish faster than 1/t for t ^ 1/Aqcd, and so the integral over t is 
convergent. The absence of endpoint divergences in convolution integrals is con- 
nected to the question of whether the infra-red degrees of freedom have correctly 
been identified in the effective theory [86]. If we had found endpoint singularities as 
— > 0, we would have to include other, long-distance modes in the construction of 
the effective theory (such as soft-coUinear modes), thus invalidating factorization. 

Finally, let us demonstrate that more complicated projections onto the B meson 
involving higher Fock states or transverse parton momenta only enter at sublead- 
ing order in power counting. This can be seen from the rules for constructing 
SCET operators out of gauge-invariant building blocks (2.60), as explained in [14]. 
Projections sensitive to transverse momentum components contain extra deriva- 
tives and so are power suppressed. Projections corresponding to non-valence Fock 
states contain insertions of the soft gluon field ^l^. Since As scales like Aqcd, such 
insertions lead to power suppression unless this field is integrated over a domain of 
extension 1/Aqcd- Reparameterization invariance dictates the such an insertion 
must be accompanied by a factor of n in the numerator, or by a factor of n in 
the denominator. The only possibility to insert a factor of n in the numerator 
is through the combination j dufi^s±{un), which vanishes since the operator al- 
ready contains a factor of yi, and = 0. For a factor of n in the denominator, we 
might insert e. g. n/{un-q), which scales like Aqcd/"^6- One might object that an 
insertion of f ■ n invalidates the above argument. However, such a factor cannot 
appear as interactions with the heavy quark can be integrated out before the final 
matching onto SCET//. This also ensures that the factor must appear to the 
left of Pj. We have thus completed the proof of the factorization formula (5.1) to 
all orders in perturbation theory, and at leading power in AqcD/"^b- 
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Figure 5.1: One-loop diagrams in the full theory contributing at leading 
power to the B —>■ jlu decay amplitude. 



5.2 Calculation of the hard-scattering kernel 



We derive the Wilson coefficients Q in (5.3) by performing a matching calculation 
using on-shell external quark states. To this end, we assign incoming momenta 
rUbV to the heavy quark and / with = to the soft light quark. By construction, 
the Wilson coefficients are independent of the nature of the external states. Soft 
gluon emission from the external quark lines cancel in the matching. However, 
emissions from the internal quark propagator in the full theory lead to the build-up 
of the finite-length Wilson line S's(tn) 5*1(0). To one-loop precision, the coefficient 
functions can be written in the form [24] 



On 



Oil H Ci + 



(5.11) 



To obtain the NLO corrections q we evaluate the one-loop contributions to the 
decay amplitude in the full theory and in SCET. The relevant diagrams in full 
QCD are shown in Fig. 5.1. In addition there is a contribution from the wave- 
function renormalization for the heavy quark. Using anti-commuting 75 (and, as 
usual, working in 4 — 2e dimensions), the diagrams evaluate to the following bare 
expressions: 



QCD 



-2e 



2e 



In^ 



2E^ /_j 
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2(1 - 2 log x^) In 



2E.A 



7 



mt 



134 






Figure 5.2: One-loop diagrams in the effective theory whose contribution to 
the amphtude needs to be subtracted in the calculation of the 
Wilson coefficients. 



4 In'' + 



2 3^'^ 



Inx^ — 2/72(1 — x^) — 2 — TT^ 
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^QCD 



-2v ■ I 



-2e 



In' 



-2v ■ I 



-2e 



111 

12 



(5.12) 



The expression receives contributions from the weak current vertex correc- 

tion and wave-function renormalization for the heavy quark (first bracket), the 
electro-magnetic vertex correction and self-energy insertion (second term), and 
the box diagram (last term, which depends on f ■ / in addition to /+ [24, 83]). 

On the low-energy theory side we evaluate the one-loop order contributions, 
which can be written as Ci_ioop ® (OscET)treo + (^^66 ® (Oscet)i-Ioop- The relevant 
diagrams are shown in Fig. 5.2 and evaluate to 



SCET 



~2t 



1 



-2v ■ I 
/i 



e2 4 

-2e 



(5.13) 



+ ln' 



-2vl 



12 



^SCET 



0. 



The first term in the expression for ^4^'^'^"'" corresponds to the first two diagrams 
in the figure, while the second term is obtained from the last graph. Note that 
this is precisely the same contribution as obtained from the box diagram in the 
full theory. The difference of the expressions given in (5.12) and (5.13) determines 
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the NLO contributions to the Wilson coefficient functions. Subtracting the pole 
terms in the MS scheme, we find 

= _2 In^ ^ + (5 - 4 In x,) In ^ + In^ - 2 In^ x, 

_^ 2-3x^ lnx^-2L2(l-x^)-7-^, (5.14) 
1 — 4 

C2 = ^ _^ ■ 

These results for the Wilson coefficients agree with the expressions for the hard- 
scattering kernel given in [84] and [85] , apart from the contribution of the HQET 
decay constant F{fi). (The authors of [84] and [85] assumed an unconventional 
mft-dependent definition of the i?- meson LCD A.) After simplification we find the 
hard-scattering kernel at next-to-leading order [24] 



T{l+,E^,mb,n) = 1 + 



Air 



21n2^ + 21n^ + ln2 ^^■^^+ 



/i /i /i2 

"'^^''"'^ 2L2(l-x.)-5 



1-x^ " 4 



(5.15) 



The renormalization-scale dependence is, despite appearance, entirely given in the 
form 21n^(/i//+) — 21n(/i//+) -|- /i-independent terms, which can therefore be can- 
celed against the fi dependence of the LCDA (p^il^, jj) under the convolution inte- 
gral in (5.1). The remaining large logarithms \n{E^/lj^) need to be resummed to all 
orders in perturbation theory. This can be done by stating the result in the above 
form, factorizing the kernel into a hard function H (which contains ln(2£'^//i)) 
and a jet function J (which contains ln(2£'^ /+/yU^)), and solving their RGEs. 

5.3 Resummation of large logarithms 

The hard-scattering kernel T can be systematically factorized when performing the 
matching calculation in two separate steps QCD SCET/ SCET77 [14, 85, 86]. 
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The hard function H{2E^/fi,x.y) arises when integrating out off-shell fluctuations 
of order nib, while fluctuations around the intermediate scale mfeAgcD lead to 
the jet function. Since the scale mbAqcD arises only from a scalar product of 
a soft momentum with a coUinear momentum, one can simply set / = in our 
previous matching calculation, which ensures that only the hard fluctuations are 
integrated out [24]. After the function H has been identified in this way, the jet 
function is constructed using that J = T/H. Specifically, for / = the second and 
third diagrams in Fig. 5.1 involve scaleless integrals that vanish in dimensional 
regularization, while the box graph can readily be shown to vanish at leading 
power. The remaining contribution from the weak vertex correction and wave- 
function renormalization for the heavy quark yields 

1 5 21nx-, n 
- - 2 In^ 

(5.16) 



^QCD 



1=0 



nib 



-2e 



2e 



1 — 



In x^ 



2L2(l-x,)-6-- 



while the expression for v4^'"°|j_q is the same as that for A^^^ given in (5.12). 
After MS subtractions the above result determines the hard function Hi defined 
in (5.8). Using the first relation in (5.10), it then follows that 
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f2E^ 



Atx 



_21n2^ + 21n^ 



/i 



2Lo(l -X. 



4-^ 
12 



(5.17) 



The knowledge of the hard-scattering kernel (5.15) and the above result for the 
hard function yields the jet function, which is thus given by 



J 



7 



1 + 



In" 



1-^ 
6 



(5.18) 



\ ji^ J 47r \ ji^ 

In order to proceed we need RG equations obeyed by the various coefficient 
functions. The fact that the decay amplitude in (5.1) is scale independent links 
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the scale dependence of the hard-scattering kernel to the evolution of the LCDA 
(j)^{uj,fi) [33]. By analyzing the renormalization properties of this function, as we 
have done in Section 3.1.1, it follows that the hard-scattering kernel satisfies the 
integro-differential equation (for > 0) [30]^ 
d 



r^„snftts) In ^ + 7(as 
+ 



cusp V- a/ — T 



T(/+, / du; 1+ r(u;, /+, a^) T(u;, /i), 

(5.19) 



where Fcusp is the universal cusp anomalous dimension familiar from the theory of 
the renormalization of Wilson loops [25] and r(ci;, u;', a^) is given in equation (3.12). 
From the functional forms of the hard and jet functions given above, it follows that 
the hard component and the jet function obey the RG equations 
d 



dln/i 

J(/+,^) 



dln/i 



-Fcusp («s) 
Fcusp (a,) lii^^-y 



+ 7("s) -7'("s) 



(5.20) 
(5.21) 



+ / dujl+T{uj,l+,as) J{uj,fi) . 



The anomalous dimensions 7 and 7' do not have a simple geometric interpre- 
tation and must be determined by explicit calculation, unlike is the case for the 
cusp anomalous dimension [51]. From (5.17) and (5.18) we find 



a.. 



47r 



(5.22) 



We now discuss the general solution of the evolution equations (5.19) and (5.20). 
Exact solutions can be written down analogously to the discussion in Section 3.1.1. 
We define the dimensionless function 



3'(a, as 



du uj' V{uj^ uj' , as) ( — 



(5.23) 



^For simplicity of notation, we omit the arguments E^, rrih, and for the 
remainder of this section. Also, unless otherwise indicated, as = as{fi)- 
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This is the same function as previously defined in (3.14), and for convenience we 
restate the one-loop order result (3.21) as 



3^(^)(«, a.) = rcusp(a.) k(l + a) + ^(1 - a) + 2^^ 



(5.24) 



where ip{z) is the logarithmic derivative of the Euler F-function. We start by 
solving the first equation in (5.20) with the initial condition for H{fih) evaluated 
at a high scale /i/, ~ mi,, for which it does not contain large logarithms. We 
then evolve the function H{fi) down to an intermediate scale /i, ~ y^^TT^AqcD and 
multiply it by the result J{l^,fii) for the jet function, which at the intermediate 
scale is free of large logarithms and can be written in the general form J(/+, fii) = 
J[as(/ij), ln(2£'^ Z+Z/if)]. This determines the kernel T{l^,fii) at the intermediate 
scale. Finally, we solve (5.19) and compute the evolution down to a low-energy 
scale fj, ~ few x Aqcd- The exact solution is given by 



T(/+, /i) = H{fih) di^sifJ^i), Vr,] exp U (/+, fi, fj,i, fih, v) 



ri=0 



(5.25) 



where the notation tJ[as(/ii), V,,] means that one must replace each logarithm of 
the ratio 2E^ l+/fJ''i by a derivative with respect to an auxiliary parameter 77. The 
evolution function U is given by [24] 
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Given this exact result, it is straightforward to derive approximate expressions 
for the kernel at given orders in RG-improved perturbation theory, by using per- 
turbative expansions of the anomalous dimensions and /^-function to the required 
order. Unfortunately, controlling terms of 0{as) in the evolution function U would 
require knowledge of all anomalous dimensions at two-loop order, which at present 
is lacking. (It also requires the cusp anomalous dimension to three loops, which 
has been computed only recently [57].) We will, however, control the dependence 
on the variables /+ and to 0{as)- As usual, we write 

n=0 "■■ n=0 

and similarly for the anomalous dimensions 7 and 7'. The relevant expansion 
coefficients are listed in (4.31) and (4.31), and we also find 70 = —2Cf, 7o = 0. 
Defining the ratios ri = as{fii)/(ys{f^h) and r2 = as(/i)/as(/ii), we obtain our final 
result [24] 
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/+ [GeV] /+ [GeV] 

Figure 5.3: RG-improved predictions for the hard-scattering kernel at maxi- 
mum photon energy. Left: Results at /i = 1 GeV. The bands refer 
to different values of the intermediate matching scale: fif = A^mb 
(center), 2Ah'mh (top), O-SA/jm^ (bottom). Their width reflects 
the sensitivity to the high-energy matching scale /i^, varied be- 
tween 2ml 0.5ml. The dashed line shows the result obtained 
at one-loop order. Right: Dependence of the kernel on the renor- 
malization scale varied between 0.75 GeV and 2.0 GeV as in- 
dicated on the curves. 

- ^ Hnr,) - i^ ni + clnr^) _ ^^^^^^^^ ^ ^^^^^ (5 39) 
2fjo r(l-clnr2) 

corresponds to the function U evaluated with 77 = 0, /+ = /i, and 2E^ = 

Finally, let us study the effects from Sudakov resummation numerically. In the 
left-hand plot in Fig. 5.3 we compare the resummed hard-scattering kernel (5.28) 
with the expanded 0{as) approximation (5.15). The effect is maximized by setting 
the photon energy to its highest possible value of = mh/2. (Numerically, we 
work with the pole mass mt = 4.8 GeV.) The renormalization scale n is chosen to 
1 GeV, which is a reasonable scale for the 5-meson LCDA. The matching scales 
Hh and fii are varied around their "natural" choices fXh = '^E^ and /ij = ^y2E■yAh, 
where Ah = 0.5 GeV serves as a typical hadronic scale, are taken as default values. 
We use the two-loop running coupling normalized at as{mh) = 0.22 and set rij = 4 
for the number of light quark flavors. (For simplicity, we do not match onto a 
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three-flavor theory even for low renormalization scales.) We find that resummation 
effects decrease the magnitude of the radiative corrections, i.e., the resummed 
kernel is closer to the tree-level value (T = 1) than the one-loop result. The 
fact that after Sudakov resummation the radiative corrections are moderate in 
magnitude persists even for asymptotically large 6-quark masses. For instance, 
setting mf, = 50 GeV we find after resummation T(/_|_, fi) = 0.74 at /_|_ = /i = 1 GeV. 
(Fixed-order perturbation theory breaks down for such large values of the quark 
mass. From (5.15) we would obtain r(/_|_,/i) = 0.08 with these parameter values.) 
The figure also exhibits that our results are stable under variation of the two 
matching scales. Varying fi] and /i| by factors of 2 changes the result for the 
kernel by less than 10%. This suggests that the unknown NNLO corrections to 
the function Uq in (5.29) are perhaps not very important. 

The scale dependence of the resummed expression for the kernel is illustrated in 
the right-hand plot in Fig. 5.3, which shows the functional dependence of T{l^,fi) 
for maximal photon energy and several values of fi. The matching scales are set to 
their default values fih = f^b = 4.8 GeV and /ij = y^Ahmt, ^ 1.55 GeV. We observe 
a significant scale dependence of the kernel, especially as one lowers fi below the 
intermediate scale /ij. In other words, the second stage of running (for fi < fii) is 
numerically significant. 



CHAPTER 6 
EXCLUSIVE SEMILEPTONIC DECAYS 



6.1 Heavy-to- light form factors at large recoil 



Transitions of the B meson into light pseudoscalar or vector mesons are parame- 
terized in terms of scalar functions, called form factors, that depend only on the 
Lorentz invariant quantity g^, where q = Pb — Pai is the momentum transfer from 
the 5-meson onto the light final meson M. They are defined by the following 
Lorentz decompositions of bilinear quark current matrix elements: 
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where mp the mass of the pseudoscalar meson. The relevant form factors for B 
decays into vector mesons with polarization vector e and mass my are defined as 
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Figure 6.1: Gluon exchange contributions to heavy-to-light form factors. 

The flavor-changing current is denoted by a dashed line. The 
lines to the left belong to the B meson, and those to the right 
belong to the light meson M. 



In many applications form factors are needed near zero momentum transfer 
(g^ ^ 0), corresponding to a kinematic situation in which a flavor-changing weak 
current turns a heavy B meson at rest into a highly energetic light meson. Such 
heavy-to-light transitions at high recoil are suppressed in QCD by inverse powers 
of the heavy-quark mass rrih. (This is contrary to the case of heavy-to-heavy 
transitions such as i? D, which are unsuppressed in the heavy-quark limit 
[2].) The challenge in understanding the physics of these processes is to describe 
properly the transformation of the soft constituents of the B meson into the fast 
moving, collinear constituents of the energetic light meson in the final state. At 
lowest order in perturbation theory this transformation can be achieved by the 
exchange of a gluon, as depicted in Fig. 6.1. Naive power counting suggests that 
the nature of the gluons is hard-coUinear, as they turn the soft spectator quark 
into a collinear constituent of the final meson of energy 

E = » Aqcd (6.7) 

in the 5-meson rest frame. The meson momentum can be written as p'^j = En^ + 
0{m\,j/AE), which is nearly light-like. We assume that the partons inside the light 
meson carry a significant fraction of its total energy, and that the light constituents 
of the B meson carry soft momenta of order Aqcd- This justifies to describe the 
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process using SCET. It would then appear that the form factor is governed by 
a sufficiently hard gluon exchange, which can be dealt with using perturbative 
methods for hard exclusive QCD processes [29, 88]. The situation is, as we shall 
see below, not quite that simple [12, 14, 27, 86, 89, 90, 91]. Leading contributions 
to the form factors also arise when the exchanged gluon is of long-distance nature, 
leading to the "soft overlap" (also called Feynman mechanism). As an example, a 
straightforward calculation of the two diagrams in Fig. 6.1 yields 

where = (1, 0, 0, 0) is the 5-meson velocity, F denotes the Dirac structure 
of the flavor- changing current, the * product means that the two factors must be 
sandwiched between quark spinors, and the dots represent power-suppressed terms. 
Here we assigned an incoming momentum / to the light spectator anti-quark, and 
the outgoing parton momenta pi and p2 collinear to the light-meson momentum 
are parameterized as 

= xiEnf" + + . . . , pf^ = X2En^ - + . . . , (6.9) 

with the fractions Xi, X2, subject to the constraint X1 + X2 = 1. Assuming that / is 
soft and Xi of order unity one flnds that the amplitude (6.8) scales like 1/(-E'^Aqcd- 
However, assuming linear distributions for X2 — > and n ■ / — > in the endpoints 
of the meson LCDAs, it is evident from (6.8) that the corresponding convolution 
integrals for a hadronic amplitude diverge at these endpoints [92, 93]. In other 
words, the above analysis breaks down for such momentum configurations because 
the gluon exchange is no longer of short- distance nature. Naively, form factors 
would then be dominated by this soft overlap, since it is not suppressed by a 
perturbative coupling constant a^. However, it has sometimes been argued that 
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the summation of large Sudakov logarithms associated with the soft gluon exchange 
mechanism may lead to a strong suppression of the soft overlap term, essentially 
reinstating the perturbative nature of the form factors (see, e.g., [94, 95, 96]). 
The assumption of a short- distance nature of heavy-to-light form factors at large 
recoil is the basis of the pQCD approach to exclusive hadronic B decays [97], 
which is often considered a competitor to QCD factorization. Leaving aside some 
conceptual problems associated with this treatment [98], the issue of Sudakov 
logarithms is an intricate one, because contributions to the form factors can arise 
from several different energy scales. Besides the hard scale E ~ Mb and the 
hadronic scale Aqcd, interactions between soft and coUinear partons involve the 
intermediate hard-coUinear scale of order ^i^Aqco- To settle the question of 
Sudakov suppression of the soft overlap contribution is one of the main motivations 
for our work [89]. 

6.2 Spin-symmetry and factorization 

The discussion of heavy-to-light form factors in the heavy-quark limit E ^ Aqcd 
can be summarized by the factorization formula [90] 

/f-A^(E) = a{E, /ii) CM(/ii, E) + T,{E, /EA^, /i) ® ® 0m(/x) + . . . , 

(6.10) 

where the dots represent terms that are of subleading order in Aqcd / E. Ci and Tj 
are calculable short-distance coefficient functions, (pB and 0m denote the leading- 
order LCDAs of the B meson and the light meson M, and the products imply 
convolutions over light-cone momentum fractions. The functions Cm denote uni- 
versal form factors that only depend on the nature of the light final-state meson 
but not on the Lorentz structure of the currents whose matrix elements define the 
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various form factors. The first term in the factorization formula therefore imphes 
spin-symmetry relations between different form factors, which were first derived 
in [18] by considering the large-energy limit of QCD. The second term in (6.10), 
which arises from hard-collinear gluon exchange, breaks these symmetry relations 
[90]. 

The arguments E and y/^AqcD in the short-distance coefficients Cj and Tj 
are representative for any of the hard or hard-collinear scales in the problem, re- 
spectively. The form of the factorization formula shown in (6.10) assumes that the 
factorization scale /ii in the first term is chosen to lie between the hard scale E 
and the hard-collinear scale ^Z-EAqcd, whereas the scale fi in the second term is 
chosen to lie below the hard-collinear scale. The Wilson coefficients Ci receive con- 
tributions at tree level, whereas the hard-scattering kernels Tj start at first order 
in asi^/EAQci))- Naively, one would conclude that the spin-symmetry preserving 
term provides the leading contribution to the form factors in the heavy-quark limit. 
Only then the notion of an approximate spin symmetry would be justified. How- 
ever, the situation is more complicated because as written in (6.10) the universal 
functions Ca/(/^I7 E) still depend on the short-distance scale ^yE Aqcd', in fact, their 
£^-dependence remains unspecified. We will see how this dependence on the large 
energy enters the amplitudes when our technology of scale separation is applied. 
The matching of a QCD form factor onto matrix elements in SCET is done in two 
steps: QCD ^ SCET/ ^ SCET//. In close analogy to the previous applications 
discussed in this thesis, hard fluctuations with virtualities on the scale E ~ Mb 
are integrated out in a first step, while in the second step hard-collinear modes 
with virtualities of order ^^JEAQCD are removed. Both steps of this matching are 
understood for the hard-scattering term in the factorization formula, for which the 
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kernels Tj can be factorized as [86] 

QCD,/^)- (6-11) 

i 

Below we will focus on the universal functions C,m{E,^i) entering the first term 
in (6.10). (The notation /ii serves as a reminder that the factorization scale for 
this term is defined in the intermediate effective theory SCET/.) In [86, 91] these 
functions are defined in terms of matrix elements in the intermediate effective 
theory SCET/, which leaves open the possibility that they could be dominated by 
short-distance physics. Here we show [89] that the functions Cm renormalized at a 
hard-coUinear scale ^ihc ~ a/^Aqcd can be factorized further according to 

CM{lihc,E) = Y,Df\VEK^,lih,,li)®if'\li,E) , (6.12) 

k 

where the functions are defined in terms of hadronic matrix elements of 
SCET// operators. By solving the RG equation for these operators we show that 
the universal form factors C,m do not receive a significant perturbative suppression, 
neither by a power of as{-\/EK^^) nor by resummed Sudakov logarithms. This 
statement holds true even in the limit where Mb is taken to be much larger than 
the physical 5-meson mass. Perhaps somewhat surprisingly, we find that there 
exists a leading contribution to the soft functions Ca/, for which the hard-collinear 
scale ^E'Aqcd is without physical significance. Switching from SCET/ to SCET//, 
one describes the same physics using a different set of degrees of freedom. This 
observation leads us to the most important conclusion of this work, namely that 
the long-distance, soft overlap contribution to heavy-to-light form factors exists. 
However, we point out that even at low hadronic scales /i ~ Aqcd the functions 
^fc^^^ contain a dependence on the large recoil energy E which is of long-distance 
nature and cannot be factorized using RG techniques. As a result, it is impossible 
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to determine the asymptotic behavior of the QCD form factors f^~*^^ as E ^ oo 
using short-distance methods. 

In order to prove the factorization formula (6.10) one needs to show that all 
contributions to the form factors that do not obey spin-symmetry relations can 
be written in terms of convolution integrals involving the leading-order LCDAs 
of the B meson and the light meson. Specifically, this means showing that (i) no 
higher Fock states (or higher-twist two-particle distribution amplitudes) contribute 
at leading power, and (ii) the convolution integrals are convergent to all orders in 
perturbation theory. Point (i) can be dealt with using the power-counting rules of 
SCET along with reparameterization invariance [86, 91]. Here we use our formalism 
to complete step (ii) of the factorization proof. (The convergence of the convolution 
integral over the S-meson LCDA can also be shown using arguments along the lines 
of the discussion around equation (5.8).) 

6.3 Matching calculations 

In the matching of the intermediate effective theory SCET/ onto the final low- 
energy effective theory SCET//, hard-coUinear modes with virtuality of order 
QCD are integrated out, and their effects are included in short-distance coeffi- 
cient functions. Our primary goal in this section is to construct a basis of operators 
relevant to the matching of the universal functions (m in (6.10) onto the low-energy 
theory (at leading power in A = Aqcb/E), and to calculate their Wilson coeffi- 
cients at lowest order in perturbation theory. As always, matching calculations can 
be done using on-shell external quark and gluon states rather than the physical 
meson states, whose matrix elements define the form factors. The results for the 
Wilson coefficients are insensitive to infra-red physics. 
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6.3.1 Spin-symmetric contributions 

Heavy-to-light form factors are defined in terms of i? — M matrix elements of 
flavor- changing currents qVb. The spin-symmetric contributions to the form fac- 
tors are characterized by the fact that in the intermediate effective theory (i.e., 
after hard fluctuations with virtuality /i ~ ~ Mb are integrated out) they con- 
tain the Dirac structure F sandwiched between the two projection operators 
and |(1 + 1^). This implies that T can be decomposed into a linear combination 
of only three independent Dirac structures, which leads to symmetry relations be- 
tween various form factors [18]. As was shown in the equations (2.42) and (2.44), 
the relevant operators in SCET/ can be written as time-ordered products of the ef- 
fective Lagrangian with effective current operators J^j^{x) = [^hcWhc\{x) Tm h{xJ) 
defined by the matching relation 

qVh = Y,Klt{m^„E,^^)Jf^ + ... , (6.13) 

M 

where the dots denote power-suppressed terms. Here Tm = l,75!7±i/ is one of 
the three Dirac basis matrices that remain after the projections onto the two- 
component spinors C,hc and h. The variable E entering the coefficients K^^ is the 
total energy carried by coUinear particles (more precisely, 2E = n ■ pl°^), which 
in our case coincides with the energy of the final-state meson M. Since the La- 
grangian is a Lorentz scalar, it follows that for 5 — M transitions each of the three 
possibilities corresponds to a particular choice of the final-state meson (hence the 
label "M" on Tm), namely Fm = 1 for M = P a pseudoscalar meson, Tm = 15 
for M = V|| a longitudinally polarized vector meson, and F^y^ = for M = V± a. 
transversely polarized vector meson. The coefficients for the various currents 
that are relevant in the discussion of heavy-to-light form factors are summarized 
in Table 6.1, where we use the definitions (2.43) for the transverse tensors g'^ and 
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Table 6.1: Coefficients Klj arising in the leading-order matching of flavor- 
changing currents from QCD onto SCET/. The coefficients Ci are 
deflned in equation (2.44) and given beyond tree-level in [9]. We 
denote q = q/Ms and E = E/Mb, where q = Pb — Pm- 



Current 


M = P 


M -- 


= V\\ 


M = V± 


qTb 


{Tm = 1) 




= 75) 




q-f^^b 


{n'^Ci + v^'C5) 








qi^ibb 






+ V'Cs) 




qia^^qyb 


K - (1 - -E)n^] Cii 






-5f [C9 + (l-^)Ci2] 


qia'^'^^qub 




K - (1 - 




-i^l{CQ + ECi2) 



e'^. Once a deflnition for the heavy-to-light form factors is adopted, it is an easy 
exercise to read off from the table which of the coefficient functions Ci contribute 
to a given form factor. This determines the functions Ci{E,fii) in (6.10). (In 
general, these functions are linear combinations of the Ci in Table 6.1.) 

Next, we match the time-ordered product i J c/^xT {jjj''(0), £scETi(a;)} onto 
operators in SCET//, focusing flrst on operators which include soft and collinear 
fields only. The insertion of a subleading interaction from the SCET/ Lagrangian is 
required to transform the soft 5-meson spectator anti-quark into a hard-coUinear 
final-state parton. In the resulting SCET// interaction terms the soft and collinear 
fields must be multipole expanded, since their momenta have different scaling 
properties with the expansion parameter A [12], see Section 2.3.1. Soft fields are 
expanded about x+ = while collinear ones are expanded about x_ = 0. In 
general, collinear fields can live on different x+ positions while soft fields can live on 
different a;_ positions. The relevant operators can be written as matrix elements of 
color singlet- singlet four-quark operators multiplied by position-dependent Wilson 
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coefficients, i.e. [14, 16] 



/ 



ds dt Dk{s, t, E, /i) [(e Wc) {x+ + x±)... (lyj {x+ + x± + sn)] 

X [{q^Ss){x.+x^ + tn)...{Slh){x_+x^)], (6.14) 



where the dots represent different Dirac structures. There is no need to include 
color octet-octet operators, since they have vanishing projections onto physical 
hadron states and do not mix into the color singlet-singlet operators under renor- 
malization. The Fourier transforms of the coefficient functions defined as 



coincide with the momentum-space Wilson coefficient functions, which we will 
calculate below. (In the matching calculation, u is identified with the component 
n-l of the incoming spectator momentum, and U2 is identified with the longitudinal 
momentum fraction X2 carried by the coUinear anti-quark in the final-state meson.) 
When we rewrite the operators above in terms of the gauge-invariant building 
blocks defined in (2.60), all factors of the soft-collinear Wilson lines Ssc{0) and 
Wsc{0) cancel out, as explained after equation (2.84), since the coUinear fields are 
evaluated at x_ = and the soft fields at x+ = 0. Hence, we can rewrite the 
operators in the form (setting a; = for simplicity) 



At leading order in A, operators can also contain insertions of transverse deriva- 
tives or gauge fields between the coUinear or soft quark fields. Since 9^ ~ A and 
A'^j_ ~ A'^j_ ~ A, such transverse insertions must be accompanied by a factor of 
i^/in ■ ds ~ A^^. The inverse derivative operator acts on a light soft field and 
implies an integration over the position of that field on the n light-cone, the effect 




(6.15) 




(6.16) 
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of which can be absorbed into the Wilson coefficient functions. The appearance 
of ifi in the numerator is enforced by reparameterization invariance, as derived in 
equation (2.87). It then follows that in our case only single insertions of transverse 
objects are allowed [14, 24]. Finally, the multipole expansion of the soft fields 
ensures that the component fi ■ ps oi soft momenta does not enter Feynman dia- 
grams at leading power. Likewise, at leading power there are no operators that 
contain in ■ dc ~ (acting on coUinear fields) or n ■ Ac ~ A^, since these are 
always suppressed with respect to the corresponding transverse quantities. The 
new structures containing gluon fields are of the form 



and we define the corresponding Fourier-transformed coefficient functions in anal- 
ogy to (6.15) as Di:{u2,U3,uj, E, ^) and Dk{u!i,u!2,U2, E, In matching calcula- 
tions, U3 is identified with the longitudinal momentum fraction X3 carried by a 
final-state coUinear gluon, while ui and U2 are associated with the components 
n ■ Iq and n ■ Ig of the incoming soft anti-quark and gluon momenta. 

We are now in a position to construct a basis of four-quark operators relevant 
to the discussion of the universal functions (m in (6.10) [89]. These operators must 
transform like the current J^^ under Lorentz transformations. Also, the soft and 
coUinear parts of the four-quark operators must have non-zero projections onto 
the B meson and the final-state meson M. We set the transverse momenta of the 
mesons to zero, in which case there is no need to include operators with transverse 
derivatives acting on the products of all soft or coUinear fields. 




dr ds dt Dk{r, s, t, E, /i) [X(0) . . . A'^^{rn) . . . X(sn)] [Q,(tn) . . . J{(0)] , 



ds dt du Dkis, t, u, E, fi) [X(0) . . . X{sn)] [Q,(tn) . . . ^^^(un) . . . J{(0)] , 



(6.17) 



(6.18) 



153 

Case M = P: The resulting operators must transform as a scalar {Tm = 1). A 
basis of such operators is 

oP = g' [X(0) 1 75 X{sn)] [Q,{tn) ^ 75 ^0)] , 
of) = [X(0) 1 75 ^^±X(sn)] [Q,M 1 75 :k(o)] , 
of = [X(0) 1 75 4c±(rn) X(sn)] [Q,(tn) |75:K(0)] , 
of) = [X(0) 1 75 X(sn)] [Q,(tn) 4s±{un) ^ 75 :K(0)] . 
The soft and coUinear currents both transform like a pseudoscalar. 

Case M = Vj]: The resulting operators must transform as a pseudo-scalar 
(Tm = 75)- A basis of such operators is obtained by omitting the 75 between 
the two collinear spinor fields in (6.18), so that the collinear currents transform 
like a scalar. The Wilson coefficients for the cases M = P and M = V|| coincide 
up to an overall sign due to parity invar iance. 

Case M = V±: The resulting operators must transform as a transverse vector 
{Tm = 7±i/) • A basis of such operators is 

of-) = / [X(0) I 7±.75 X(sn)] [Q,{tn) ^ 75 ^(0)] , 

of-) = / [X(0) ^ {ze^^ - ^?i75) ^^1 X{sn)] [Q.(tn) ^ 75 n^)] , 

(6.19) 

of-) = / [X(0) I (.e;i-„ - ^^75) X(.n)] [Q.(tn) | 75 ?<(0)] , 

of-) = g' [X(0) I 7±.75 X(sn)] [Q,(tn) ^.^(Mn) ^ 75 ^K(O)] . 

The soft currents transform like a pseudoscalar, while the collinear currents trans- 
form like a transverse axial-vector. (The relative sign between the two terms in 
of-) can be determined by considering left and right-handed spinors X and using 
relation (6.22) below.) 
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It is not necessary to include operators with a derivative on the soft spectator 
anti-quark. All such operators would contain Qs {—i^i) ii-, which can be related to 
a linear combination of the operators o(*^^ and 04*^^ using the equation of motion 
for the light-quark field. 

Wilson coefRcients of four-quark operators 

We now calculate the momentum-space Wilson coefficients of the operators O^-^^ at 
tree level [89]. The relevant graphs must contain a hard-coUinear gluon exchange, 
which turns the soft 5-meson spectator anti-quark into a collinear parton that can 
be absorbed by the final-state hadron. These hard-coUinear gluons are integrated 
out when SCET/ is matched onto SCET//. However, at 0{as) we can obtain 
the Wilson coefficients in SCET// by directly matching QCD amplitudes onto the 
low-energy theory, without going through an intermediate effective theory. 

As we have seen, at leading power the universal form factors C,m receive contri- 
butions from ordinary four-quark operators as well as from four-quark operators 
containing an additional collinear or soft gluon field. We start with a discussion 
of the matching calculation for the operators 0^12-1 whose coefficients can be ob- 
tained by analyzing the diagrams shown in Fig. 6.1 in the kinematic region where 
the outgoing quarks are collinear and the incoming quarks are soft. The resulting 
expression for the amplitude in the full theory has already been given in (6.8). 
Note that, by assumption, n ■ / ~ Aqcd and 0:2 ~ 1, so that the result is well 
defined and it is consistent to neglect subleading terms in the gluon propagators. 

We wish to express the result for the QCD amplitude as a combination of 
SCET// matrix elements multiplied by coefficient functions. To this end, we must 
first eliminate the superficially leading term in (6.8) using the equations of motion. 
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Assigning momenta to the two outgoing collinear lines as shown in (6.9), it follows 
that 

r * = 1 7„r * - 2r * t + . . . , (6.20) 

where the dots denote power-suppressed terms. In the next step we use a Fierz 
transformation to recast the amplitude into a form that is convenient for our anal- 
ysis. Taking into account that between collinear spinors the Dirac basis contains 
only three independent matrices, we find for general matrices M and 

2{u^Muh) {v^Nv^) = {u^'^v^) {vqN'i^Muh) + {^'^l^v^) [v^N-iJ^Muh] 

+ {uii\l±.vd {v.N^i'I^Muh) . (6.21) 
To simplify the Dirac algebra we use the identities 

7±75 i = ie^±l±u i , 7±7l i = {dT " ^^Tl^) i ■ (6-22) 

Finally, we include a minus sign from fermion exchange under the Fierz transfor- 
mation and project the collinear and soft "currents" in the resulting expressions 
onto color-singlet states. 

Once in this form, the different contributions to the amplitude can be readily 
identified with matrix elements of the operators O'^^ 'i- For M = P, Vjj we obtain 
n(-P) _ n(^i) - 1 + ^2 p(P) _ n^^i^ - ^ 

(6.23) 

while for M = Vj_ we find 

(y^) _ Cf_ 1 ^(vi) _ Cf_ 1 

^ " AT 4E2^2^' ^2 - ^ 4^2^^ ^2- 

Wilson coefRcients of four-quark operators with an extra gluon 

The matching calculation for the operators Og^/^ proceeds along the same lines. 
However, in this case it is necessary to study diagrams with four external quarks 
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* > ( * * * ^ > > ^ * > * ■ • y- 

H 



Figure 6.2: Diagrams relevant to the matching calculation for operators con- 
taining an extra collinear gluon. 



and an external gluon, which we treat as a background field. Let us first discuss 
the case with three collinear particles in the final state. The relevant QCD graphs 
are shown in Fig. 6.2. Physically, they correspond to non-valence Fock states of 
the light final-state meson, but this interpretation is irrelevant for the matching 
calculation, which can be done with free quark and gluon states. 

By assumption, each of the three collinear particles carries large momentum 
components along the n direction. We denote the corresponding longitudinal mo- 
mentum fractions by Xi (quark), X2 (anti-quark) and X3 (gluon), where Xi -|- X2 + 
X3 = 1. The calculation of the diagrams in Fig. 6.2 exhibits that at leading power 
the amplitude depends only on the light-cone components n ■ I and n-pi = 2Exi of 
the external momenta, and that only transverse components of the external glu- 
ons fields must be kept. These observations are in accordance with the structure 
of the operator 0^^^\ One must add to the diagrams shown in the figure a con- 
tribution arising from the application of the equation of motion for the collinear 
quark fields that led to (6.20). One way of obtaining it is to include diagrams with 
gluon emission from the external collinear quark lines in the matching calculation. 
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Performing the Fierz transformation and projecting onto the color singlet-singlet 
operators 0^^\ we obtain the coefficient functions 



D 



-D 



(^i) 



D 



1 



^-1 

U2 



Us 
U2 



J_ _ 2Cf {u2 + UsY 
m N U2{ui+Us)\ ' 

2 1 (M2 + M3) 



(6.25) 



m U2{Ui+U2)_ 

Next, we calculate the contributions from diagrams with three external soft 
particles, which physically correspond to three-particle Fock states of the B meson. 
Again there are eight diagrams, analogous to those in Fig. 6.2. The initial soft gluon 
is attached to either an off-shell intermediate line or a collinear line. The diagrams 
with external gluons must again be supplemented by a contribution resulting from 
the application of the equation of motion VqlJ = 0{g) used in the analysis of the 
four-quark amplitude in (6.8). After Fierz transformation and projection onto the 
color singlet-singlet operators 0^^'^\ we find the coefficient functions 



D 



(P) 



-D 



(^|) 



1 



D 



8E^ul{ui+U2y 
1 



1 



2Cf \uJ2 I 



1 + 



1 U2\ UJ2 



1 1 



iV2 uj LUi iV2 m 



(6.26) 



8E^ul{uJi + uj2y 

The variables Ui and 002 correspond to the light-cone components n ■ Iq and n ■ Ig 
of the incoming momenta of the soft anti-quark and gluon, respectively. 



Endpoint singularities 

The Wilson coefficients of the SCET// four-quark operators become singular in 
the limit where some of the momentum components of the external particles tend 
to zero. When the B ^ M matrix elements of the operators O^.*^'' are evaluated, 
these singularities give rise to endpoint divergences of the resulting convolution 
integrals with LCDAs. For instance, the matrix elements of the operators 0(^^^ 
and 04^^'' involve the leading twist-2 projection onto the light meson M. The 
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corresponding LCDAs are expected to vanish linearly as ^2 0, whereas 
the corresponding coefficients contain terms that grow like I/m^, giving rise to 
logarithmically divergent convolution integrals. Similarly, the operators 02^^'^ and 
Og^^'' involve the leading-order projection onto the i?-meson LCDA which is 
expected to vanish linearly as u ^ [30]. Once again, logarithmic singularities 
arise because the corresponding coefficients contain terms that grow like l/u"^. 

While the logarithmic divergences in the convolution integrals could be avoided 
by introducing some infra-red regulators, they indicate that leading-order contri- 
butions to the amplitudes arise from momentum regions that cannot be described 
correctly in terms of collinear or soft fields. In Section 6.4, we will explain that 
in SCET// these configurations are accounted for by matrix elements of operators 
containing the soft-collinear messenger fields. 

6.3.2 Spin-symmetry breaking contributions 

The two-particle amplitude in (6.8) also includes contributions for which the Dirac 
structure is different from i/iT h. These give rise to symmetry-breaking contribu- 
tions to the form factors. We shall not derive a complete basis of all possible 
symmetry-breaking operators (there are many) but rather list the ones that enter 
at ffist order in a^. Since the spin-symmetry violating terms can be factorized in 
the form of the second term in (6.10), they are associated with a short-distance 
coupling constant as{\J EKqqj^). It is therefore appropriate in this case to include 
the coupling constant in the Wilson coefficient functions. 

Let r denote the Dirac structure of the flavor- changing currents g F 6, whose 
matrix elements define the form factors. To ffist order in a^, the spin-symmetry 
breaking terms can then be obtained from the matrix elements of two operators 



159 



given by 



( 1 ^ 



Qi = [X(o) 



75 



( 1 \ 



X(sn)\ |Q.(in)|7„ 



( 1 \ 



-75 



Q2 = [X(o) 



X{sn)] [Qs{tn) 



4 



I 1 \ 

rj{(o)] 



(6.27) 



-75 



V / 



75 

V / 

where each hne contributes for a different final-state meson M. The corresponding 
Wilson coefficients are 



Cf 7ra,, 



To = ^ 



(6.28) 



N 2MbEu2 ^ N E'^U2 uj 

Linear combinations of Ti and T2 determine (up to prefactors) the hard-scattering 
kernels Tj in (6.10). The matrix elements of the operators (5i,2 can be expressed 
in terms of the leading-order LCDAs of the B meson and the light meson M. 
Only the 5-meson LCDA called 0^"* contributes [33] because of the factor ^ next 
to Qs- This property holds true to all orders in perturbation theory [86]. The 
resulting convolution integrals are convergent. Evaluating these matrix elements 
we reproduce the spin-symmetry breaking terms obtained in [90]. 



6.4 Physics of endpoint singularities — a toy model 

Our strategy in the previous sections has been to perform matching calculations 
by expanding QCD amphtudes in powers of AqcdZ-E"; assuming that collinear 
and soft external momenta have the scaling assigned to them in SCET. We have 
then matched the results onto SCET// operators containing soft and collinear 
fields and read off the corresponding Wilson coefficient functions in momentum 
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Figure 6.3: Triangle subgraph whose spectral function can be used to study 
endpoint singularities on the coUinear side. 

space. A problematic aspect of this procedure has been the observation that, if 
the matrix elements of the effective-theory operators are expressed in terms of 
meson LCDAs, then the resulting convolution integrals do not converge. Endpoint 
singularities arise, which correspond to exceptional momentum configurations in 
which some of the partons inside the external hadrons carry very small momentum. 
The question naturally arises how one should interpret these singularities, and 
whether the results we found for the short- distance coefficient functions are in fact 
correct. 

The fact that endpoint configurations are not kinematically suppressed points 
to the relevance of new momentum modes. In the limit X2 — * 0, the scaling 
associated with the collinear anti-quark in the final state of Fig. 6.1 changes from 
(A^, 1, A) to (A^, A, . . . ). Likewise, in the limit n-l ^ 0, the scaling associated with 
the soft spectator anti-quark in the initial state changes from (A, A, A) to (A^, A, . . . ). 
The soft-collinear messenger fields in the SCET// Lagrangian have precisely the 
scaling properties corresponding to these exceptional configurations. (Since the 
modes in an effective theory are always on-shell, the transverse components of 
soft-collinear momenta scale like X^^^. We will see below that this is not really 
relevant.) 
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The purpose of this section is to analyze the interplay of the various modes 
present in SCET/j and to see how soft-coUinear messengers are connected with 
the phenomenon of endpoint singularities. We will do this with the help of a toy 
example. Consider the triangle subgraph enclosed by the box in the diagram shown 
on the left-hand side in Fig. 6.3, which is one of the two gluon-exchange graphs 
relevant to the B ^ M form factors (see Fig. 6.1). To understand the physics of 
endpoint singularities we focus on the side of the light meson M (an analogous 
discussion could be given for the i?-meson side). We study the discontinuity of the 
triangle diagram in the external coUinear momentum p^. The resulting spectral 
density q{p^) models the continuum of light final-state hadrons that can be pro- 
duced on the coUinear side. We will not bother to project out a particular light 
meson from this spectral density. For simplicity, we will also ignore any numera- 
tor structure in the triangle subgraph and instead study the corresponding scalar 
triangle, as discussed in [15, 89]. 

Let us define the discontinuity 



1 

2Tn 



I{p' + iO) - I{p' -iO) =D- e{p') (6.29) 

of the scalar triangle integral 

/ = zTT-'^/V'"'' f d^k — — (6.30) 

in c? = 4 — 2e space-time dimensions, where p is the outgoing coUinear momentum 
and / the incoming soft momentum. It will be convenient to define the invariants 
L? = —P — iO and = — (/ — p^ = 2/+ ■ p^ + . . . , which scale like ~ and 
Q"^ ~ A. (In physical units, ~ Aq^^^ and ~ -EAqcd with E > Aqcd-) We 
will assume that these quantities are non-zero. From [15], we can obtain explicit 
results for the discontinuity D and for the various regions of loop momentum k 
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that give a non-vanishing contribution. We find that 



D = \n^ + 0{e,\) 



(6.31) 



and that the momentum configurations that contribute to this result are those 
where the loop momentum is either collinear, meaning that k ~ (A^, 1, A), or soft- 
coUinear, meaning that k ~ (A^, A, A'^^^). The contributions from these two regions 



which add up to the correct answer. 

It is instructive to rewrite these results in a more transparent form. To this 
end we use Cutkosky rules to evaluate the discontinuities of the diagrams directly 
and perform all phase-space integrations except the integral over the light-cone 
component n-k oi the loop momentum, which we parameterize as n ■ /c = —X2n ■ P- 
As in previous sections, X2 denotes the fraction of longitudinal momentum carried 
by the anti-quark in the final-state. The exact result for the discontinuity of the 
scalar triangle is {x2 = 1 — X2) 



The collinear and soft-collinear contributions separately are divergent even though 
they correspond to tree diagrams (after the two propagators have been cut). We 
obtain 



are 




(6.32) 



D 




(6.33) 




(6.34) 
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The collinear contribution is infra-red singular for X2 ^ 0, which is an example of 
an endpoint singularity. In the present case, the singularity is regularized dimen- 
sionally by keeping e non-zero. The soft-coUinear contribution is infra-red finite 
but ultra-violet divergent, since X2 runs up to 00. Again, this divergence is reg- 
ularized dimensionally. Evaluating the remaining integrals one recovers the exact 
results in (6.32). 

Several comments are in order: 

i) The result for the spectral density in the full theory is finite. The endpoint 
singularity is regularized by keeping subleading terms in the hard-collinear prop- 
agator l/[—{k + /)^] ^ l/(x2Q^ + L"^)- When the subleading terms are dropped 
based on naive power counting (as we did in the analysis of the previous sections) 
the full-theory result reduces to the contribution obtained from the collinear re- 
gion. An endpoint singularity arises in this case, which however can be regularized 
dimensionally. In order for this to happen in the realistic case of external meson 
states, it would be necessary to perform the projections onto the meson LCDAs 
in d A dimensions. The factor {x2X2)~'^ in (6.34) would then correspond to a 
modification of the LCDAs, which renders the convolution integrals finite. 

ii) The collinear approximation fails for values X2 = 0(A). For such small 
momentum fractions the exact full-theory result is reproduced by the contribution 
obtained from the soft-coUinear region. In fact, the collinear and soft-coUinear 
contributions in (6.34) coincide with the first terms in the Taylor expansion of 
the full-theory result in (6.33) in the limits where X2 = 0(1) ^ L'^/Q^ and X2 = 
0(A) ^ 1. The fact that X2 runs up to 00 in the soft-coUinear contribution 
is not a problem. In dimensional regularization the integral receives significant 
contributions only from the region where X2 ~ L'^/Q'^ = 0(A). To see this, one 
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can introduce a cutoff 6 to separate the collinear and soft-coUinear contributions, 
chosen such that 1 ^ 6 ^ X. This cutoff is introduced as an infra-red regulator 



This yields Dq = - In 5 and Dsc = ln6 + \nf^ + 0{\/6). The sum of the two 
contributions once again reproduces the exact result. 

iii) Next, note that the soft-coUinear contribution precisely reproduces the 
endpoint behavior of the full-theory amplitude. It is irrelevant in this context that 
the soft-coUinear virtuality /c^^ ~ E'^X^ is parametrically smaller than the QCD 
scale Aqqj-). What matters is that the plus and minus components of the soft- 
collinear momentum, kgc ~ (A^, A, . . . ), are of the same order as the corresponding 
components of a collinear momentum in the endpoint region, where n-p ~ A rather 
than being 0(1). 

iv) Finally, we see that the endpoint divergences we encountered in the Sec- 
tion 6.3.1 were not regularized because we dropped the dimensional regulator when 
performing the projections onto meson LCDAs. This, however, does not affect the 
results for the Wilson coefficients of the SCET/t- operators, which remain valid. 
In the toy model, the collinear contribution is represented in the effective the- 
ory as the discontinuity of the integral (corresponding to the first diagram on the 
right-hand side in Fig. 6.4) 



which shows that the Wilson coefficient resulting from integrating out the hard- 
collinear propagator is simply —l/u2- 



in the collinear integral {J^ J^) and as an ultra-violet regulator in the soft- 
coUinear integral (Jq°° J^). The integrals can then be evaluated setting e — 0. 




(6.35) 
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Figure 6.4: CoUinear and soft-coUinear contributions to the scalar triangle 
graph in SCET//. Dashed (dotted) lines represent coUinear (soft- 
coUinear) propagators. 

To summarize this discussion, we stress that to reproduce the behavior of the 
full amplitude it is necessary to include in the effective theory contributions in- 
volving coUinear fields and those involving soft-coUinear messengers, as indicated 
in Fig. 6.4. The latter ones represent the exact behavior of the amplitude in the 
endpoint region. Whether or not endpoint configurations contribute at leading 
order in power counting is equivalent to the question of whether or not operators 
involving soft-collinear fields can arise at leading power. This shows the power 
of our formalism. Operators containing soft-collinear fields can be used to pa- 
rameterize in a systematic way the long-distance physics associated with endpoint 
configurations in the external mesons states. This will be discussed in more detail 
in the next section. It also follows that the scaling properties of operators contain- 
ing soft-collinear fields can be used to make model-independent statements about 
the convergence of convolution integrals that arise in QCD factorization theorems 
such as (6.10). We will exploit this connection in Section 6.7. 

6.5 Soft-collinear messengers and the soft overlap 

In Section 6.3 we have derived four-quark operators built out of coUinear and soft 
fields, which contribute at leading order to the universal soft functions Cm in (6.10). 
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Power counting shows that the products of these operators with their Wilson co- 
efficients scale like A^. (Here one uses that (ir ~ c/s ~ 1 and dt du ^ ^^^0 
Taking into account the counting of the external hadron states, \B) ~ A^'^^^ and 
(M| ~ A"^, it follows that the corresponding contributions to the universal form 
factors scale like A'^^^. In other words, heavy-to-light form factors are quantities 
that vanish at leading order in the large-energy expansion. The fact that the re- 
sulting convolution integrals were infra-red singular suggested that there should be 
other contributions of the same order in power counting, which cannot be described 
in terms of collinear or soft fields. 

An important observation made in the previous section was that the endpoint 
behavior of QCD amplitudes is described in the low-energy theory in terms of di- 
agrams involving soft-coUinear messenger fields. For instance, in the last graph in 
Fig. 6.4 the soft spectator anti-quark inside the B meson turns into a soft-collinear 
anti-quark by the emission of a soft gluon. The soft-collinear anti-quark is then 
absorbed by the final-state meson. In a similar way, endpoint singularities on the 
S-meson side would correspond to a situation where the initial state contains a 
soft-collinear spectator anti-quark, which absorbs a collinear gluon and turns into 
a collinear anti-quark. The SCET// Lagrangian contains the corresponding inter- 
action terms only at subleading order in A. However, because the universal form 
factors Ca/ themselves are power-suppressed quantities, these subleading interac- 
tions will nevertheless give rise to leading-power contributions. 

We will need the first two non-vanishing orders in the interactions that couple 
a soft-collinear quark to a soft or collinear quark. The results are most transparent 
when expressed in terms of the gauge-invariant building blocks introduced in (2.60). 
They are [15] 
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4f ^ = Qs 4s±Wl 9 , 4^') = a Ssc 4c± X , (6.36) 

and 

2 tn ■ o 2 

(6.37) 

where 

(r = -f-—^'ipsc±9 (6.38) 
2 in ■ Dsc 

contains the "small components" of the soft-coUinear quark field, which are in- 
tegrated out in the construction of the effective Lagrangian. The longitudinal 
components of the calligraphic gluon fields are defined in [15]. (Note also that 
n- As = Q and n ■ Ac = The soft-coUinear fields enter these results in combina- 
tions such as Wl^ 6 or Sl^ 9, which are gauge invariant. Soft and coUinear fields live 
at position x, while soft-collinear fields are evaluated at position x^ for C^g and x_ 
for Cg^. The measure d'^x associated with these interactions scales like A~^. The 
superscript on the Lagrangians indicates at which order in power counting (A^^^ 
or A) the corresponding terms contribute to the action. 

Next, we need the representation of the fiavor- changing SCET/ current J^j in 
(6.13) in terms of operators in the low-energy theory SCET//. As shown in [16], 
at leading power, and at the matching scale /i = fihc, the relation reads 

4? (a;) ^ JmH^) = + Ta/ {Si WscM + x±) (6.39) 

with a Wilson coefficient equal to unity. The anomalous dimensions of the currents 
are the same in the two theories. The product {S^ Wsc){0) arises since soft-collinear 
messenger fields cannot be decoupled from the current operator J'^j^ in SCET//, 
contrary to the case of the color singlet-singlet four-quark operators discussed in 
Section 6.3. 
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Figure 6.5: An artist's view of the soft-coUinear messenger contribution to 
the form factors. The shaded region contains soft-colhnear inter- 
actions. The arrows indicate the flow of the components n ■ Ps 
(left) and n ■ (right) of soft and coUinear momenta, which do 
not enter the soft-coUinear block. 

Using these results, we can write down a tri-local operator whose matrix ele- 
ment provides a long-distance contribution to the universal form factors. It is 

or = d^xd^yT{4J'\.) £«(,) jS\0) + 4;)(x) 4f (,) jS\o)} . 

(6.40) 

Even after the decoupling transformation (2.60) this operator contains arbitrarily 
complicated soft-collinear exchanges, as illustrated in the cartoon in Fig. 6.5. 

Note that the superficially leading term in the time-ordered product cancels 
[15]. This can be understood as follows: 

After the decoupling transformation the strong-interaction part of the SCET// 
Lagrangian no longer contains unsuppressed interactions between soft-collinear 
messengers and soft or collinear fields. In order to preserve a transparent power 
counting it is then convenient to define hadron states in the effective theory as 
eigenstates of one of the two leading-order Lagrangians Cs and Cc- For instance, 
we define a "SCET pion" to be a bound state of only collinear fields, and a "SCET 
B meson" to be a bound state of only soft fields. The SCET pion state coincides 
with the true pion, because the collinear Lagrangian is equivalent to the QCD 
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Lagrangian [12].^ The SCET B meson coincides with the asymptotic heavy-meson 
state as defined in heavy-quark effective theory. It is important in this context 
that time-ordered products of soft-colhnear fields with only collinear or only soft 
fields vanish to all orders as a consequence of analyticity [15]. Hence, soft-colhnear 
modes do not affect the spectrum of hadronic eigenstates of the collinear or soft 
Lagrangians. 

Once the SCET// hadron states are defined in this way, each term in the time- 
ordered product (6.40) can be factorized into a part containing all soft and collinear 
fields times a vacuum correlation function of the soft-collinear fields. These vacuum 
correlators must be invariant under rotations in the transverse plane. It follows 
that the correlator arising from the superficially leading term in the time-ordered 
product vanishes, since 

{n\ T {{aSMy-) (slWschkio) (wleux^)} \n) (6.4i) 

contains a single transverse derivative (see (6.38)). 

The leading terms in the time-ordered product in (6.40) scale like A^, since 
^m"* ~ A^/^. They thus contribute at the same order in power counting to the 
universal functions as the four-quark operators discussed in Section 6.3. When 
Og*^'' is added to the list of four-quark operators a complete description of the 
soft-overlap contribution is obtained. The Wilson coefficient of this new operator 
follows from the fact that there is no non-trivial matching coefficient in (6.39), 
that the current operators have the same anomalous dimensions in SCET/ and 

SCET//, and that the SCET// Lagrangian is not renormalized. Hence, to all 

^The endpoint region of the pion wave function is not described in terms of a 
Fock component containing a soft-collinear parton, but rather in terms of a time- 
ordered product of the SCET pion state with an insertion of the Lagrangian Cg^. 
This insertion is non-zero only in processes where also soft partons are involved. 
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orders in perturbation theory 



^5 if^hcfi) 



(6.42) 



which is in fact a universal function, independent of the labels "i" and "M". 
As before, Hhc denotes the hard-coUinear matching scale, at which the transition 
SCETj SCET// is made. The appearance of Ci{^hc) in the denominator of this 
relation is due to the fact that this coefficient was factored out in the definition of 



If we define by the B ^ M hadronic matrix elements of the operators 
0^^\ then the sum Cm = Xlfc -^i*^'' ^i^^"* describes the entire soft overlap con- 
tribution. Each term in this sum gives a contribution of order A^^^ to the form 
factors, in accordance with the scaling law obtained a long time ago in the con- 
text of QCD sum rules [99]. We have thus completed the derivation of the new 
factorization formula (6.12). Whereas the Wilson coefficients D^^^^ depend on the 
renormalization scale /i as well as on the hard-coUinear scale ^hc ~ ^-EAqcd, the 
characteristic scale of the hadronic matrix elements is the QCD scale Aqcd, 
not the hard-coUinear scale. While this is obvious for the matrix elements of the 
operator 0^^^\ it also holds true for the remaining matrix elements, for which the 
sensitivity to long-distance physics is signaled by the presence of endpoint singu- 
larities. It remains to discuss how our results are affected by single and double 
logarithmic corrections arising in higher orders of perturbation theory. 

6.6 Operator mixing and Sudakov logarithms 

Because the operators O^*^"* share the same global quantum numbers they can mix 
under renormalization. This mixing is governed by a 5 x 5 matrix of anomalous 




see (6.10) and (6.12). 
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dimension kernels, which are in general complicated functions of the light-cone 
variables Ui and cUj. The anomalous dimension matrix governing this mixing has 
the structure 

/ 

7ii 7i4 

722 723 
732 733 
741 744 



(6.43) 



V 751 752 753 754 755 / 

Because the operators 0[^^l consist of products of soft and coUinear currents (after 
decoupling of soft-collinear messengers), the entries in the upper left 4x4 sub- 
matrix can all be written as sums of soft and collinear anomalous dimensions 
for the corresponding current operators. We have taken into account that the 
operators 02^^^ and Og^^'* mix under renormalization (this mixing is determined by 
the mixing of twist-3 two-particle and three-particle LCDAs for light mesons), as 
do the operators 0[^^^ and 0[^^^ (this mixing has not yet been worked out but is 
allowed on general grounds). The operator O^'^^^ consisting of a triple time-ordered 
product requires the "local" operators as counter-terms; however, those operators 
do not mix back into 0^'^^\ 

Because of the structure of the anomalous dimension matrix it follows that 
the coefficient 755 is one of the eigenvalues, which governs the scale dependence 
of the Wilson coefficient in (6.42). This coefficient coincides with the well- 
known anomalous dimension of the current operators J^j^f and J'^^^ computed in 
Section 2.4.2. The corresponding "operator eigenvector" can be written as a linear 
combination (in the convolution sense) 

4 



(M) 



oi*" 



)(A/) 



(6.44) 



k=l 
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where the coefficients d\!^^^ are independent of the high-energy matching scale fihc- 
The other four eigenvectors are hnear combinations of the operators O^*^"* with 
k ^ 5. This observation has an important consequence: The decomposition of the 
form factors in terms of matrix elements of these eigenvectors implies that there 
exists a contribution from Oil^^\ The remaining contributions can be written 
as endpoint-finite convolution integrals that involve both two- and three-particle 
LCDAs of the initial and final mesons. Leaving such factorizable terms aside, it 
follows that the combination of operators contributing to the universal form factors 
must, up to an overall factor, coincide with the eigenvector Ol^^\ and hence 

df) = ^ (6.45) 

to all orders in perturbation theory. With a slight abuse of notation, let us now 
denote by the B ^ M hadronic matrix element of the eigenvector O^^^^ in 
SCETj/. Combining (6.10) and (6.42), we then find that the spin-symmetric imi- 
versal form-factor term can be rewritten as 

C^{E,^l)CM{,^J'l,E)\ = Ci{E, fi) CuifJ', E)\ (6.46) 



II 



This relation is not as dull as it seems; rather, it contains the remarkable message 
that for the soft overlap contribution to heavy-to-light form factors the interme- 
diate hard-collinear scale is without any physical significance. Switching from 
SCET/ to SCET// we merely describe the same physics using a different set of de- 
grees of freedom. In other words, there is no use of going through an intermediate 
effective theory. The RG evolution of the soft functions (m remains the same all 
the way from the high-energy scale E ~ Mb down to hadronic scales fi ~ Aqcd- 
The physics of the soft overlap term is thus rather different from the physics of the 
spin-symmetry breaking corrections in the factorization formula (6.10), for which 
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the hard-coUinear scale is of physical significance. Any spin-symmetry breaking 
contribution involves at least one hard-collinear gluon exchange, and the amplitude 
factorizes below the scale fihc- 

The result (6.46) allows us to systematically resum the short-distance loga- 
rithms arising in the evolution from high energies down to hadronic scales. The 
RG equation obeyed by the Wilson coefficient functions is [9, 16] 

■a{E,fi) = freu.p[as(/i)] In— + 7[a,(/i)]^ Ci{E,fi), (6.47) 



d m /i \ /i 

where the coefficient of the logarithmic term is determined by the cusp anomalous 
dimension [25]. Its solution is 

Q(E, /i) = Q(E, 11^) exp f/(/i,,, /i, E) , (6.48) 

where fih ~ '^E is the high-energy matching scale for the transition from QCD 
to SCET, at which the values of the Wilson coefficients can be reliably computed 
using fixed-order perturbation theory. The RG evolution function can be written 
as 



f/(/i,,/i,E)= / d«.^^"^P^"^ 



as(^lh) 



, 2E^ f da' 
In 



as{fJ,h) 



da^ (6.49) 



where l3{as) = dag/dln fi. The dependence on the high-energy matching scale 
Hh cancels against that of the Wilson coefficients Ci{E,fih) in (6.48). Note that 
after exponentiation the evolution function contains an energy and scale-dependent 
factor exp U{fih,f^,E) oc where 

aW^j da^. (6.50) 

In order for this scale dependence to be canceled, the low-energy hadronic matrix 
element must carry an energy dependence of the form [Kqct> / EY'^^\ 
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Evaluating the hadronic matrix elements at a low scale means that all short- 
distance dependence on the large energy E is extracted and resummed in the 
coefficient functions. However, in the present case the matrix elements still con- 
tain a long-distance dependence on the large scale E, which cannot be factorized 
[16]. The reason is that the large energy enters the effective theory as an exter- 
nal variable imprinted by the particular kinematics of soft-to-collinear transitions. 
Because of the large Lorentz boost 7 = E/rriM connecting the rest frames of the 
B meson and the light meson M, the low-energy effective theory knows about the 
large scale E even though hard quantum fluctuations have been integrated out. 
This is similar to applications of heavy-quark effective theory to 6 — * c transitions, 
where the fields depend on the external velocities of the hadrons containing the 
heavy quarks, and '-f = vi, ■ Vc is an external parameter that appears in the ma- 
trix elements of velocity-changing current operators [26, 2]. In the present case, 
it is perhaps reasonable to assume that the primordial energy dependence of the 
hadronic matrix elements at some low hadronic scale might be moderate, so that 
the dominant E dependence is of short-distance nature and can be extracted into 
the Wilson coefficient functions. However, there will always be some energy de- 
pendence left in the matrix elements; even if we assume that it is absent for some 
value of fi, it will unavoidably be reintroduced when we change the scale. As a 
result, it is impossible to determine the asymptotic behavior of the QCD form 
factors ff~*^\q^) using short- distance methods.^ 

Let us now proceed to study the numerical importance of short- distance Su- 

dakov logarithms. Given the exact results in (6.48) and (6.49), it is straightforward 

^The same phenomenon is known to occur in the case of the Sudakov form 
factor, for which the coefficient of the double logarithm is sensitive to infra-red 
physics [31, 100]. 



175 



to derive approximate expressions for the resummed Wilson coefficients at a given 
order in RG-improved perturbation theory by using perturbative expansions of the 
anomalous dimensions and /5-function. Unfortunately, controlling terms of 0{as) 
in the evolution function U would require knowledge of the cusp anomalous dimen- 
sion to three-loop order (and knowledge of 7 to two-loop order), which at present 
is lacking. We can, however, control the dependence on the recoil energy E to 
0{as). Following [24], we define the ratio r = as{^i) / asil^h) and obtain 



^ ijL In r 

2E\^^ 



^^h) 



where 

47r 



(6.51) 



4/32 \_as{iih 



f 1 Inr I + In^r 

V ^ ; 2/?o 

(r — 1 — Inr 



(6.52) 



Inr + 0(a,). 

^Po 



To /3o. 

The only piece missing for a complete resummation at next-to-next-to-leading 
order is the 0{a^ contribution to f/o, which is independent of E. The relevant 
expansion coefficients are Fq = y, Fi = — yvr^ — ^rif, 70 = — y, and 
j3o = 11 — ^Uf, Pi = 102 — yn/. We set n/ = 4 in our numerical work. 

In Fig. 6.6 we show the dependence of the Wilson coefficients Ci{E,fi) on the 
large energy E. We choose fih = 2£' for the high-energy matching scale and 
use the tree-level initial conditions Ci{E,fih) = 1, which is consistent at next-to- 
leading order. We fix /i at a low hadronic scale in order to maximize the effect 
of Sudakov logarithms. The maximum recoil energy in i? ^ vr transitions is such 
that 2£'max — 5.3 GeV. Obviously, for such values the perturbative resummation 
effects are very moderate. In the energy range 1 GeV < 2E < Mb the Wilson 
coefficients differ from unity by no more than about 20%. The extrapolation to 
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Figure 6.6: Energy dependence of the Wilson coefficients Ci{E,fi) at next- 
to-leading order in RG-improved perturbation theory. The three 
curves correspond to as(/i) = 1 (solid), as{fi) = 0.75 (dashed- 
dotted), and as(/i) = 0.5 (dashed), which are representative of 
typical hadronic scales. 

larger energy values shows that Sudakov suppression remains a moderate effect 
even for very large recoil energy. 



6.7 Factorization of spin-symmetry breaking effects 

Using the close connection between messenger exchange and endpoint singulari- 
ties discussed in Section 6.4, the formalism of soft-collinear fields can be used to 
demonstrate the convergence of convolution integrals in QCD factorization theo- 
rems. If messenger exchange is unsuppressed, the convolution integrals diverge at 
the endpoints, spoiling factorization. By the same reasoning, convolution integrals 
are finite if soft-collinear messenger contributions are absent at leading power. 

Let us apply this method to show, to all orders in perturbation theory, that the 
convolution integrals entering the spin-symmetry breaking term in the factorization 
formula (6.10) are free of endpoint singularities. As mentioned in the Introduction, 
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this is an essential ingredient still missing from the proof of this formula. With our 
technology the proof is rather simple and consists of only the following two steps: 

1. Soft-collinear messenger fields can be decoupled from the four-quark op- 
erators mediating spin-symmetry breaking effects, which are of the type shown 
in (6.27). The reason is that in the color singlet- singlet case the operators are 
invariant under the field redefinition (2.60) [16]. It follows that their matrix ele- 
ments factorize into separate matrix elements of soft and collinear currents, which 
can be written in terms of the leading-order LCDAs of the external mesons. If 
the messengers did not decouple, they would introduce unsuppressed interactions 
between the soft and collinear parts of the four-quark operators, thereby spoiling 
factorization. 

2. Time-ordered products containing interactions of soft-collinear messengers 
with soft or collinear fields do not contribute to the spin-symmetry violating term in 
(6.10). This follows from SCET power counting. The insertions of the Lagrangians 
Cqe and Cq^ in (6.40) cost a factor A'^^^, meaning that they can only come together 
with a leading-order current operator of the type X . . . 3-C ~ A^^^. (Note that there 
are no 0(A^/^) corrections to the current that could make the vacuum correlator 
(6.41) involving C^^J"^^ and C^^^^'^^ non-zero [16].) However, such a current will 
always be of the form (6.39) because of the projection properties of the heavy- 
quark and coUinear-quark spinors. It thus respects the spin-symmetry relations. 

These two observations imply that at leading order in Aqcd / E soft-collinear ex- 
changes do not contribute to the spin-symmetry breaking term in the factorization 
formula, and hence the corresponding convolution integrals are convergent. 



CHAPTER 7 
CONCLUSION 



QCD-Factorization theorems are statements about leading power properties of 
amplitudes to all orders in the strong coupling constant. The use of effective field 
theories enables one to prove such theorems rigorously, and thus take a significant 
step further from past fixed-order proofs. In this thesis we demonstrated this tech- 
nique on some decay modes that are particularly clean in the sense of simplicity 
of QCD-factorization formulae. The first part of this work was devoted to a de- 
tailed introduction to Soft-CoUinear Effective Theory (SCET), which describes the 
infra-red physics of strong interactions between light, but highly energetic partons 
with soft heavy or light degrees of freedom. The crucial advantage of this theory 
is that it allows for a systematic expansion in inverse powers of the heavy-quark 
mass already at the Lagrangian level, therefore allowing for all-order proofs. Fur- 
thermore, SCET enables us to perform perturbative calculations for systems in 
which local operator expansions break down due to long-distance effects. This is a 
significant improvement of our ability to compute the strong-interaction effects of 
weak decays of heavy mesons. Furthermore we investigated two non-perturbative 
structure functions: the i?-meson light-cone distribution function, which enters the 
calculation of exclusive decay amplitudes in the high-recoil region, and the shape 
function, which encodes the Fermi-motion of the heavy quark inside the i?-meson 
for the calculation of inclusive decay modes. In particular, we have derived impor- 
tant relations that link the (finite-interval) moments of the shape function to the 
values of local HQET operators. 

To gain full control over the separation of physics of different energy scales, 
it is necessary to perform a resummation of large logarithms. We have demon- 
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strated this technique for inclusive B Xul^v and exclusive ^ and 
B — s> Ml~u (where M is a light pseudoscalar or vector meson), using a sophisti- 
cated multi-step matching procedure and solving renormalization-group equations 
(RGEs). In the inclusive decay mode we matched QCD SCET/ ^ HQET. The 
intermediate theory SCET/ contains soft and hard-coUinear degrees of freedom, 
the latter being momentum modes with large energy E and moderate invariant 
momentum square of order E'Aqcd- In a second step, these modes are removed, 
resulting in a description within heavy-quark effective theory (HQET). The triple 
differential decay rate was then expressed at leading power in Aqcd / E and at next- 
to-leading order in renormalization-group improved perturbation theory in terms 
of Wilson coefficient functions H and J, convoluted with the matrix element of 
a single HQET operator, the shape function 5*. This formulation is valid in the 
shape-function region of phase space, where the Xu system is energetic and light, 
as anticipated. Model-independent results have been derived for various event 
fractions and spectra, including a cut on the charged-lepton energy, the hadronic 
variable = Eh — \Ph\, and the hadronic invariant mass. Finally we have com- 
puted numerical predictions using models for the shape function that are consistent 
with all known constraints. As a result, event fractions needed for the extraction 
of the CKM matrix element \Vub\ are found to be much larger than previously 
anticipated. This in turn implies that the numerical value of \Vub\ determinations 
from inclusive B decays might become smaller. To make a conclusive statement 
it will be necessary to include subleading power corrections (and possibly higher- 
order corrections) to the predictions derived here. For future determinations we 
proposed to use the P+ variable to discriminate against the charm background. A 
detailed study of the current theoretical uncertainties and the nature of the charm 
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background was also performed. 

QCD-Factorization was proved to all orders in the coupling constant and to 
leading power for the exclusive ^ 7 /~z/ decay amplitude. This mode provides a 
clean and realistic environment in which to study the aforementioned methodology. 
We performed a detailed matching calculation QCD SCET//. Large logarithms 
were resummed by first factorizing the hard-scattering kernel into a hard and a 
jet function. This is accomplished by matching in two steps QCD SCET/ 
SCETij. While the first matching step gives rise to the hard function, the jet 
function arises in the second step. Sudakov resummation was then achieved by 
solving the RGEs of the hard function and the combined hard-scattering kernel. 
We found that the resummation effects enhance the amplitude (as opposed to a 
"Sudakov suppression" ) . 

Finally we turned our attention to amplitudes that do not factorize. Form 
factors encode the exclusive decay amplitudes for B Pl^v and B Vl^v 
{P = light pseudoscalar meson, V = light vector meson). There are three differ- 
ent form factors for 5 — vr and seven for B p transitions. In the large recoil 
limit they obey approximate spin-sjTiimetry relations (corrections of which factor- 
ize), leaving only three independent universal form factors. In our formalism, the 
proof of this statement was rather simple; however, our analysis revealed also that 
there exists a soft contribution to the form factors which cannot be factorized. 
In SCETn this can be understood as the soft-coUinear sector does not decouple. 
As a result, there exists a long-distance operator O5, which mixes into the short- 
distance operators under renormalization. Since the short-distance operators do 
not mix back into O5, the anomalous dimension of O5, which is identical to the one 
of heavy-to-collinear current operators, is also an Eigenvalue of the entire anoma- 
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lous dimension matrix. The matrix element of the corresponding "Eigenvector" , a 
hnear combination of all operators, is endpoint-finite and contributes to the form 
factor without any suppression. In fact, the knowledge of the anomalous dimen- 
sion Eigenvalue allows to control the short-distance dependence on the large recoil 
energy. We found that there is no sizable Sudakov suppression to the universal 
form factors. In summary, the spin-symmetric soft overlap exists, and dominates 
over the factorizable hard-scattering contributions that break spin-symmetry. 

In light of these few examples of B decay processes, the power of the QCD- 
Factorization approach to strong-interaction effects with more than one relevant 
physical energy scale is apparent. While we are unable to solve QCD analytically, 
an approximation to leading power in the expansion parameter Aqco/'^f) can be 
of great help in our ability to predict Standard-Model processes. In some cases 
(e.g. inclusive decays) we are even provided with a framework in which corrections 
to the leading-power factorization formulae can be calculated systematically. In the 
future, we will be able to use such frameworks to improve the precision of theoreti- 
cal input quantities necessary for the extraction of Standard Model parameters and 
for searches of New Physics effects. We shall mention here the important examples 
of flavor-changing neutral decay modes such as -B ^ 7 and B ^ Xg , which 
are particularly sensitive to possible New Physics, but suffer from sizable QCD- 
uncertainties which can spoil valuable information. It will certainly be worthwhile 
to apply the general methodology presented in this Thesis and in the correspond- 
ing references to these processes, as well as to improve the accuracy of our results 
in both the power expansion and the coupling expansion to higher orders. On the 
exclusive side we have dealt with a few exemplary decay channels. It will be ex- 
citing to study how the lessons we learned are applicable to decay modes of bigger 
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phenomenological interest, such as the radiative B -ft'*7, or the non-leptonic 
B ^ TXTi and B Kn. As far as searches for New Physics effects are concerned, 
precision in Flavor Physics can provide us with valuable information. In both the 
experimental and theoretical community this fact is reflected in the ongoing dis- 
cussion on the potential of a "Super 5-factory" , that might complement the direct 
searches of New Physics at next-generation particle colliders. 



APPENDIX A 
LIST OF ABBREVIATIONS AND ACRONYMS 



c - coUinear 

CKM - Cabibbo, Kobayashi, Maskawa 

CLCG - Collinear Light-Cone Gauge 

h - hard 

he - hard-colhnear 

HQET - Heavy-Quark Effective Theory 

NLO - Next-to-Leading Order 

NNLO - Next-to-Next-to-Leading Order 

pQCD - perturbative QCD 

QCD - Quantum Chromo Dynamics 

RG - Renormahzation Group 

RGE - Renormahzation-Group Equation 

s - soft 

so - soft-colhnear 

SCET - Soft-CoUinear Effective Theory 

SLCG - Soft Light-Cone Gauge 

us - ultrasoft 
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